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Tests at Allied Chemical show 
Nitrogen Tetroxide ready for in- 
stant use even after 9 years of 
storage. Combustion tests with 
important liquid fuels and Nitro- 
gen Tetroxide indicated up to 
99% theoretical combustion effi- 
ciency. Additional advantages 
offered by NoOgare listed at right. 

We'll gladly supply technical 
literature which includes a 59- 
page Product Bulletin, a bro- 
chure entitled “Large Scale Han- 
dling of Nitrogen Tetroxide”, 
and “Nitrogen Tetroxide as an 
Oxidizer in Rocket Propulsion”. 
And, of course, Allied Chemical 
technical service is always avail- 
able to users of NoOs. 
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fuels—ignition is simple and 
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carbon. 

@ Eliminates rough starts—which 
are due to accumulation of un- 
reacted propellants in thrust 
chamber. 

@ Motors are throttleable—when 
N20s4 is used as the oxidizer. 


e Immediate availability— Allied 
Chemical Nitrogen Tetroxide is 
shipped in 125- and 150-Ib. steel 
cylinders, 1-ton containers, 
and 50-ton tank cars. 
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flea GREAT advances which have been made in the de- 


sign of space missiles in the past few years have empha- 
sized the feasibility of establishing a space observatory for 
tlhe measurement and study of radiations from the celestial 
bodies. Moreover, the problem of the design and construc- 
tion of instruments that are functional and efficient in a space 
eivironment has already been largely solved. The develop- 
nent of dependable and accurate pointing controls for small 
aid large instruments has proceeded to the point where 
minor, though difficult, details remain to be worked out. 
I, addition, one cannot overlook the importance of obtaining 
scientific data from measurements at vantage points outside 
Farth’s atmosphere. We have only to follow the accelerated 
progress of upper air physics and astrophysics since the time 
o: rocket and satellite measurements, crude though some of 
them are, to realize the impact of the new dimension that has 
been added to man’s means of studying the physical universe. 

The motivations behind this movement into space are 
many and varied. Certainly one of them is based on the 
natural curiosity of man about the world around him. In 
this age of practicality and purposeful thinking it is probably 
more necessary than ever that we make concessions to the 
space scientists who wish to further our basic understanding 
by pursuing the accepted methods of scientific research in a 
direction which might seem at the moment to be unncessary 
and even some what naive. 

The purpose of this paper is to present a realistic picture 
of the instrumentation of a possible space solar observatory 
and to sketch some of the engineering and scientific problems 
behind the functioning of such an observatory. Emphasis 
will be placed on optical devices for the study of solar electro- 
magnetic radiation from the x-ray region of the solar spectrum 
through the ultraviolet, visible and infrared. We are not 
concerned with the missile techniques needed for the success- 
ful launching of such an observatory or for its stabilization. 
Problems (1, 2)! of pointing the instruments and their tempera- 
ture control are not discussed. Some time will be spent 
indicating how the data will contribute to the advancement 
of astrophysics. Mention will be made of some current 
results about solar radiation as obtained from rocket and 
satellite flights. In conclusion specific instrumentation of an 
imaginary space solar observatory will be briefly considered. 

This paper is a modified version of ARS Preprint no. 795-59, 
presented at the ARS Semi-Annual Meeting, June 8-11, 1959, 


San Diego, Calif. 
' Numbers in parentheses indicate References at end of paper. 
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Photon Counters 


A photon counter is (3) a modified Geiger-Mueller device 
adapted to the detection of ultraviolet and x-ray radiation 
over either a broad or a narrow band. It has been success- 
fully employed as a rocket instrument for the measurement of 
solar radiation. The essential features are shown in Fig. 1. 
Its action is as follows: The window cuts off radiation below 
a certain wave length and permits the remaining photons to 
enter a chamber which is filled with an appropriate gas mix- 
ture at a suitable pressure. The photons may cause ioniza- 
tion in the gas, or they may strike the cathode material and 
produce photoelectrons. Either event may initiate a Town- 
send discharge and cause a count. The voltage across the 
anode and cathode, the gases used, their partial pressures, 
and the nature of the cathode material are important factors 
in determining the detailed response. Low energy photons 
which get through the window cannot ionize the gas, but can 
cause the emission of electrons from the surface of the cathode; 
for these longer wave lengths the quantum vield of the tube 
may be of the order of one count per 10° photons. As the 
energy of the photons increases, the internal photoelectric 
effect occurs at the cathode, and the yield rises rapidly unless 
the gas is a strong absorber in this wave length region. By 
the time photo-ionization of the gas sets in, the quantum 
yield may be as high as one count per 5 X 10° photons. A 
rise to about one count per 100 photons is achieved as soon 
as the energy of the photons is sufficient to ionize the gas. 
Experiments indicate that electrons emitted by the ionization 
process are much more effective in inducing discharges, and 
hence counts, than those released by photoelectric emission 
at the cathode. 

Other factors besides the value of the ionization potential 
are important in determining the kind of gas with which the 
tube is filled. For example, if a trace of a halogen gas is 
present along with the bulk gas (usually a rare gas such as 
neon or helium) the quantum yield is known to be greatly 
reduced in certain wave length regions, usually at the very 
short (gamma ray) and the relatively long (1250 A) wave 
length regions. The partial pressure of the halogen gas has a 
great deal to do with this quenching effect; in general this 
effect increases with increase in partial pressure. HCl, NO 
and other gases which can acquire electrons easily have 
effects similar to those of the halogens. It is believed that the 
molecules of gases which exhibit the forementioned type of 
quenching effect combine with some of the electrons created 


: William A. Rense studied engineering physics at Case Technology in Cleveland, where he received 
the B.S. degree in 1935. His M.S. and Ph.D. degrees in physics were received at Ohio State Uni- 
versity in 1937 and 1939, respectively. In the following years he taught at several universities, in- 
cluding Louisiana State University, Rutgers University and the University of Colorado, where he is, 
at present, Professor of Physics and Director of the Upper Air Laboratory of the Physics Department. 

; For the past 10 years he has been active in research in the fields of solar ultraviolet spectroscopy and 


APRIL 1960 


upper air physics. 


“6 


WINDOW 


ANODE SUPPORT 
INSULATOR 


by the incident photons and carry them to the ‘anil Ww ithout 
inducing a discharge. In the case of the halogen molecules 
the electron capture is accompanied by dissociation, and the 
atoms carry the electrons. 

A knowledge of these principles and the parameters of the 
available cathode materials, windows and gases enables one 
to design a photon counter which will detect a desired spectral 
band, narrow or broad. 


Normal and Grazing Incidence Grating Spectro- 
graphs 


For dispersion and image formation involving radiation of 
wave lengths less than 1100 A, prisms and lenses cannot be 
used because of the absence of known bulk materials that 
transmit in that spectral region. Recourse must be made to 
reflection optics; reflection gratings replace prisms, and 
mirrors replace lenses. For rocket and satellite measure- 
ments, grating spectrographs have been designed which are 
basically the same as those previously used in laboratory 
research in vacuum ultraviolet spectroscopy. 

There are two main obstacles which one faces when dealing 
with reflection optics: First, the reflectivity of most coatings 
is very low in the far ultraviolet except when angles of in- 
cidence are large, and second, the elimination of aberrations, 
especially when grazing incidence is involved, is a difficult 
task. Both of these problems have been tackled since the 
advent of space solar ultraviolet spectroscopy, and have been 
partially solved. Some of the recent interesting information 


concerning the solar spectrum has been obtained as a result of 
the new design of reflection-type spectrographs. 


LIGHT SOURCE 


Fig. 1 A diagram of a typical photon counter after Chubb and Friedman (3). 
through the gas in the tube, and finally strikes the cathode surface. 


RADIUS OF GRATING, G, =R 


Radiation enters the window, passes 
The action is similar to that of a Geiger counter 


had 


The concave reflection grating, because of its self-focusi: » 
powers (thus eliminating the need for additional mirrors r 
lenses), is especially suitable for ultraviolet spectrograp!| 
Jn addition to the concave grating, all that is required is as. t 
and a detector. Fig. 2 shows diagrammatically the Rowlan | 
circle-type mounting. The slit, the center of the gratin : 
and the spectrum receiving surface are all on a circle calle | 
the Rowland circle whose diameter is equal to the radius «i 
curvature of the grating. Light from the source pass: s 
through the slit, strikes the grating and is brought to a focus 
on the Rowland circle. Actually, because of astigmatisn , 
there are two focal points involved, one on the Rowland circ. 
and the other beyond the circle. The first corresponds t > 
horizontal focus for an image of the slit, and the second to 
vertical focus for the source. Referring to Fig. 2 for the 
meaning of the symbols, one can write the equations whic!) 
locate the two astigmatic images 


r= Reosa 


1 _ cosa, 1 
r 


One notes that for very small angles of incidence and diffrac- 
tion the two images are close to coincidence when the source 
is at the slit. Likewise for angle of incidence of about 45 deg, 
with a source at infinity, a stigmatic image can be found (4). 
Insofar as photographic detection is concerned, this fact helps 
decrease the length of the lines of a spectrum, and hence 
increases thesdensity of the image. A spectrograph utilizing 
this feature is the normal-incidence type (5); for such solar 


HORIZONTAL 
AND CIRCLE Focus 
VERTICAL FOCUS > 


Fig. 2 A diagram of the essential optical parts of a concave grating ultraviolet spectrograph. Some of the light from 
the source passes through the slit S onto the grating G. The spectrum (horizontal focus) is formed along the so- 
called Rowland circle whose diameter is equal to the radius of the grating. Vertical focus of a line would occur at 7” 
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spectrographs the sun’s image must be focused on the slit by 
an auxilliary mirror to insure stigmatic quality. On the 
other hand, it is a well-known fact that radiation in the 
extreme ultraviolet is not reflected with high efficiency at small 
angles of incidence. One would therefore like to set a grating 
for large angles of incidence, in particular at angles 80 deg 
or greater (grazing incidence) because under these conditions 
reflectivity is high in the far ultraviolet. For grazing in- 
cidence astigmatism is very great, and the spectral lines are so 
long that considerable loss in overall speed is experienced 
when photographic detection is required. Thus the dilemma 
is, for normal incidence, small aberration but low reflectivity 
in the far ultraviolet; for grazing incidence, high reflectivity 
in the far ultraviolet, but large aberration. 

lhe problem of low reflectivity in the far ultraviolet has 
ben tackled by several groups (6-8) with encouraging results. 
Fiims have been found which increase the reflectivity near 
normal incidence. Zinc sulfide films and magnesium fluoride 
co.tings on aluminum are examples. Unlike the case of 
aliiminum, whose reflectivity in the ultraviolet falls off rapidly 
with age, these films retain their high reflectivities for long 
periods. Thus at 1000 A reflectivity for aluminum at near 
normal incidence is 9 per cent, whereas for a zinc sulfide 
coiting the reflectivity is 18 per cent. The reflectivity for the 
magnesium fluoride coating is even higher, rising to about 
50 per cent. 

\ method has been worked out for reducing the astigmatism 
in a grazing-incidence spectrograph (9). Referring to Equa- 
tions [1], one finds that for r, equal to infinity (sun) the condi- 
tion that the horizontal and vertical focus both fall on the 
Rowland circle is impossible for grazing incidence angles when 
gratings with reasonable curvature are employed. In fact, 
light must be converging on the grating in order to eliminate 
astigmatism. Thus, substituting r.’ = r’ in Equations [1] 


one finds that ~~ 

[2] 
For a typical case where R = 50 em, a = 85 deg, 8 = 78 deg 
(corresponding to 300 A radiation when the grating constant 
is 6000 lines/em) one finds that the source distance for vertical 
focus on the Rowland circle is about —11 em, and therefore 
that a virtual source must be used. For obtaining a virtual 
image of the sun at the proper distance behind the grating 
from the slit, a toroidal mirror of suitable radii of curvature 
may be positioned in front of the slit. To insure high reflec- 
tivity in the far ultraviolet, the mirror, like the grating, is 
set for grazing incidence. From the equations for the dis- 
tances of the astigmatic foci of a spherical mirror it is possible 
to choose the radii of the toroidal mirror so that the horizon- 
tally focused image of the sun is at the slit and the verti- 
cally focused image is at the required point beyond the grat- 
ing. Thus, for the case referred to above at 300 A, the 
radii are 46 and 2.7 cm, respectively. 

A rocket grazing incidence spectrograph is shown in Fig. 3. 
In this instrument the film (SWR Eastman Kodak) is fed 
through the Rowland circle onto the roll-up spindle in the 
steel cassette. The exposure times for the successive spectro- 
grams during a rocket flight are determined by a timer 
which automatically opens a shutter after the film is 
moved between exposures. The grating is a replica on red 
glass, uncoated, made from a master of 6000 lines/em and 
radius 100 cm by Bausch and Lomb Optical Co. 


re = R/(cos a — sin B tan B) 


High Resolution Spectrographic Devices 


High resolution spectra of celestial sources yield very 
valuable data for the astrophysicist whose equations often 
call for detailed information on line profiles, distribution of 
intensity of light in the source, exact wave lengths and 
related parameters. Attempts have been made to develop 
high resolution spectrographs for rocket spectroscopy, 


APRIL 1960 


- alpha line was discovered on this flight (Fig. 4). 


Fig. 3 A rocket ultraviolet spectrograph of the grazing-incidence 

type. The slit, grating, Rowland circle and film winding mech- 

anism are visible. Film is automatically wound into a light-tight 

steel cassette. The photoeyes which were used in pointing the 

instrument at the sun are visible on the grating end of the spec- 
trograph 


though most of the research to date has centered around high 
speed, low resolution instruments. 

The resolving power for a given wave length ) is defined as 
the ratio of \ to AX, where Ad is the difference in wave 
length between two spectral lines which can just be resolved 


by the detector. Many factors help control the resolving 
power for a given spectrograph. The slit width is one; 
others are the grating itself and the kind of detector used. 
For rocket instruments employed recently, the slit width is — 
the real limiting element, along with the grating constants. 
When the grating determines the resolving power, one can _ 
valeulate the approximate value of \/AX by multiplying the | 
spectral order of the line by the number of rulings on the 
grating. A typical rocket spectrograph grating might have a 
constant of 12,000 lines/em and a ruled area 3 X 3em. For 
the first order the resolving power might therefore be as high — 
as 36,000. For solar hydrogen Lyman-alpha radiation | 
(1216 A) a spectral separation of 0.034 A could be achieved. | 
For higher orders a correspondingly smaller separation — 
could be effected. The possibility of the use of high orders | 
to increase the resolving power is limited by the ability of the | 
manufacturer of the grating to blaze the grating, that is, to — 
shape the reflection grooves so that they can concentrate the | 
light in a given direction. Other limitations arise from over-_ 
lapping orders (ultraviolet lines falling on the strong con-— 
tinuum in the visible) and scattered light. Some of these — 
difficulties were overcome in a rocket instrument recently — 
flown (10). The object was to obtain high resolution spectra — 
of the solar hydrogen Lyman-alpha line at 1216 A. Two 
gratings were used, both blazed, one for the first order and — 
the other for the thirteenth order. A magnesium fluoride — 
coating increased the reflectivity of the gratings which were | 
set for nearly normal incidence. The first grating focused — 
the line in the first order on the slit of the second grating 

which formed the final stigmatic image on specially sensitized 

film. In this way, overlapping orders and scattered light — 
problems were reduced. The absorption center in the Lyman- _ 


Perhaps the most interesting type of high resolution spectio-_ 
graph suitable for satellite and rocket observations is one | 
which employs a special kind of grating known as a reflection © 
echelle.2, The use of the latter enables the designer to con- — 


? Now made available by Bausch and Lomb Optical Co. as a = 
result of a technique of ruling developed by G. R. Harrison at — 
a. 
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Fig. 4 A high dispersion photograph of the bydrogen Lyman- 
alpha line at 1216 A taken by a Naval Research Laboratory group 
(10). The instrument was a spectrograph employing a grating 
set for the thirteenth order at Lyman-alpha. A photometric trac- 


ing of the line is shown below the photograph 
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ECHELLE 


Fig. 5 


Illustration of the principle of the echelle grating 


_ propriate angle of incidence. 


ceive of a small compact instrument that can attain resolving 
powers of the order of several hundred thousand or even a 
million in the visible and ultraviolet. Optically there are two 
ways of regarding the echelle; one may consider it to be « 
coarse ruled grating with a strong blaze in one direction 
corresponding to a high order, or one may regard it as : 
multiple reflection interference device. With the aid of 
Fig. 5 we will describe it from the latter point of view. In 
Fig. 5a the echelle is represented as a steplike structure having 

= 0.6105 X 10-3 cm, t = 2s fora typical case. There ma 
be as many as 10‘ reflection strips for the whole echell 
Assuming that a plane wave strikes the strips normal to th: 
s surface, the diffracted light from each strip is strongly 
concentrated along the normal to the latter. A mirror, ss 
in Fig. 5b, concentrates the light from all strips to a commo: 
focus. If A = 2t, the phase difference between adjacent 
rays is 360 deg and the order m is one. However, ¢ may b: 
about 1.221 X 10-* em, so that for the ultraviolet (1216 A 
the order is about 201. At normal incidence, then, th» 
resolving power at 1216 A for a typical echelle having 104 ste) 
would be (104) (201) or 2 X 108. Two lines having a separ: - 
tion of less than 0.001 A could be resolved in this ultraviolet 
region. For angles of diffraction not zero, as shown in Fig. 
5c, the same order calls for smaller wave lengths, so that th» 
spectrum of the source will be spread along a direction pe:- 
pendicular to the lines of the steps. Also, at normal 
cidence one can find a wave length for which reinforcement 
occurs by shifting to a different order: (m + n)\X,. 
where n is 1, 2, 3, ete. To prevent overlapping of orders :: 
regular grating is ‘‘crossed”’ with the echelle so as to separat:: 
the various orders at right angles to the echelle dispersion 
A working arrangement is shown in Fig. 6. Here the light 
from the sun is focused on the slit and strikes a concave 
grating set for the Wadsworth-type mounting.* The grating 
disperses the radiation and collimates it. The echelle then 
introduces dispersion at right angles to that produced by 
the concave grating and redirects the light to the latter. A 
focus is effected near the entrance slit. To reduce scattered 
light the collecting mirror can be another concave grating 
which throws its spectrum in a selective manner on the slit. 
With a resolving power of 2 X 108 it would be possible to 
detect a Doppler shift corresponding to a relative motion of 
the source of about 0.15 km/sec. 

Echelle spectrographs suitable for solar ultraviolet de- 


- tection from rockets or satellites are now under construction. 


One employing a calcium fluoride prism for the near solar 


ultraviolet has already been successfully flown (11). 


3’ For such a mounting a distant monochromatic source is 
stigmatically imaged on the norm: il to the grating for the ap- 


GRATING 


ECHELLE 


3 Fig. 6 The optical arrangement of the essential parts of a space 
observatory high resolution echelle spectrograph 
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When the echelle is to be regarded as a grating (12), the 
regular grating equation applies 


md = d(sini + sin j) [3] 


where d = (s? + ¢*)'/*. When the incident ray comes in at 
right angles to the s surface, then ¢ = x where x = tan~'t/s; 
then for a diffracted ray also at right angles to the latter 
surface the equation above becomes md = 2d sin x, which re- 
duces to 2¢ in accordance with the previous results. It 
might also be mentioned that, when the incident rays strike 
the s surface at a grazing angle @ in the horizontal plane, 
the usual formulas hold, except that s must be replaced by an 
effective value s’ equal to s cos 6, and that the slit must be 
rotated by an angle 7 given by tan 7’ = 2¢/s sin 6. The 
fact that @ need not be zero deg without deleterious optical 
effects means that the echelle may be utilized in a variety of 
mountings. For example, by setting @ large (grazing in- 
cilence) good reflectivity in the ultraviolet is possible and 
high resolution of the solar 303.8 A Hell line may be achieved. 


Detection of Infrared 


Infrared radiation is emitted by all bodies whose tempera- 
tures are not absolute zero. For this reason its measurement 
poses special problems because a detector receives the radia- 
tion from many sources besides the one being studied. Con- 
sequently, except for such cases as strong, narrow band, 
infrared sources which filters can help isolate, one goes 
about detecting infrared radiation from a given source by 
first measuring the total radiation from the source and its 
surroundings and then the radiation with the source ‘cut 
off” as effectively as possible. It is just as though one were 
to attempt to measure a visible light source with an optical 
instrument whose inner walls and baffles were coated with 
luminous paint. 

By infrared is usually meant that part of the electro- 
magnetic spectrum which extends from the visible at 7000 A 
to the microwave region where wave lengths are of the order of 
1 mm or greater. The customary unit for measurement of 
wave length in the infrared is the micron (1 wis 104A). Black 
bodies emit infrared radiation in increasing amounts as the 
temperature is raised; although only for relatively low 
temperatures does the wave length corresponding to the 
peak of the radiation continuum fall in the infrared. 

An infrared detector is characterized by two quantities, 
the responsivity and the detectivity. The former is the 
voltage per watt per square centimeter on the detector; 
the latter is basically the signal to noise ratio of the device. 
Most infrared detectors (13-15) are solids of the photo- 
conducting type. The electric conductivity changes when 
infrared radiation strikes the sensitive surface. In many cases 
the radiation is chopped before striking the receiver so that 
a-c amplification techniques can be utilized for detection. 
When cooled, these detectors have increased sensitivity. 

A very sensitive detector with unusually rapid response is 
the germanium “gold-doped”’ photoconductor developed at 
Westinghouse. One part of gold is added to about one 
hundred million parts of germanium. The detectivity is 
high between 1 and 9 p. 


Between 1.8 and 4 yu the class of infrared photoconductors ° 


which includes lead sulfide, lead telluride and lead arsenide is 
effective, especially when cooled. In fact, miniature cooling 
systems have been designed for them. 

Only bolometer-type infrared detectors are effective in the 
fur infared, and these have been employed in rocket instru- 
ments (14) to measure radiation of wave lengths between 9 
and 30. The response time for bolometers is not fast unless 
the elements are very small; on the other hand bolometers 
are not good black bodies in the infrared unless they are 
relatively thick. In most instances the desired response of 
the bolometer to infrared radiation is that of a change of 
resistance with heating, or the generation of an emf (16). 
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A typical bolometer (14) is the flake thermistor bolometer 
with a 1-mm? receiving area and a 1.5-millisec time constant. 
The unit consists of two flakes connected in series. Bias is 
supplied by a d-c potential, and the signal is drawn from the 
junction point. One flake acts as the compensating flake. 
This detector has a high output voltage for a given intensity 
of incident radiation and gives nearly constant spectral 
sensitivity over the infrared region. No cooling is required. 
However, compared to other infrared detectors, the time 
constant and detectivity are not always satisfactory. 

Dispersion in the infrared region may be effected by re- 
flection gratings or prisms made of suitable materials. Mirrors 
are satisfactory for collimating or focusing purposes, since 
mirror surface coatings may be made with high reflectivity 
in the infrared. Lenses are also available for infrared optics. 
For compact, fast instruments the spectral range to be studied 
is best isolated by means of filters. Both band pass and low 
pass infrared filters (17-20) have been designed and made. 
They consist of alternate layers of materials of high and low 
index of refraction vacuum-deposited on a supporting sub- 
strate. Ordinary electromagnetic interference theory such as 
that applicable to the Fabry-Perot interferometer explains 
the action of the filter. Analog computer methods (19) 
have been developed for carrying out the calculations in- 
volved in the design of a multilayer filter. 


X-Ray Spectrographs 


Solar x-rays and ultraviolet radiation have been detected 
in the range 100 to 1 A by photon counters (3, 21) and by 
grazing incidence grating spectrographs (22). The latter 
are not effective below about 10 A. Rocket and satellite 
x-ray spectrographs employing crystals have been designed by 
several groups including Ball Brothers Research Corp. 
(Rosendahl) and Air Force Cambridge Research Center 
(Miezaika). The general principle of such spectrographs is 
similar to that of the grating spectrograph with a bent crystal 
acting as the diffracting element. Bragg’s law gives the 
diffraction angle. For mica in the first order the crystal 
plane spacing is 9.845 A. Mica is suitable to about 18 A 
beyond which wave length its dispersion becomes unmanage- 
able, For the range extending from this point up to about 
100 A no known inorganic crystals have sufficiently high values 
of d to overcome the difficulty. However, organic crystals 
such as that of melissic acid used by Dauvillier are suitable if 
they do not have too low a melting point. 


Image Formation 


In order to transmit images of the sun or parts of the sun 
from the space observatory to Earth one must solve the 
problems of (a) image formation by optical devices, (b) 
amplification of the image if it is intrinsically faint, and (c) 
transmission of the image by electromagnetic waves to Earth. 
We will be concerned here only with (a) and (b). 

Of interest is the formation of an image of the sun’s disk in 
monochromatic ultraviolet radiation. Both lithium fluoride 
(23) and grating cameras (24) have been used for this purpose. 
In the former (25), a prism separates the radiation corre- 
sponding to the various solar ultraviolet emission lines by 
dispersion; a second prism then neutralizes the dispersion 
before a lens forms the final image. In a form of the grating 
camera (24) a stigmatic image of the sun is formed in radiation 
corresponding to a given solar ultraviolet line (such as the 
intense hydrogen Lyman-alpha line at 1216 A) by a concave 
grating mounted in the Wadsworth-type arrangement. A 
small aperture encompassing this image prevents scattered 
light from entering the final image region. The final image is 
formed by a second grating ruled on a suitable toroidal 
surface. (The toroidal surface may be produced from a 
concave grating by mechanical stress of the mount.) Aberra- 
tions are present which prevent the achievement of high — 


) 
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abilities power. Other forms of grating cameras have been 
designed (26). 

In the visible and infrared region where filters can be 
found, the problem of good image quality is much more easily 


solved. In the extreme ultraviolet a camera must employ 
gratings set at grazing incidence in order to get high reflect- 
ance. Here the reduction of aberrations is even more serious, 
but some success can be realized by proper choice of grating 
constants, spacing and radii. 

A faint image formed by the optical elements cannot be 
transmitted to Earth unless an image amplifier of some 
kind is used. Of course, if the faint image were photographed 
and automatically developed, the image could be transmitted 
by an electronics scanning device. A better procedure 
would be to utilize an image amplifier first, and then to 
transmit the picture by television techniques or some modi- 
fication thereof. 

Several manufacturers are experimenting with image 
amplifiers. One of the latter made by Radio Corporation 
of America and based on Lollemand’s method will be briefly 
described in order to display the principles and kinds of 
problems involved (27). A multi-alkali photocathode receives 
the faint image from the optical component part of the 
camera. The ejected electrons are focused by an electro- 
static lens on a phosphor screen separated from another 
photocathode by a thin film of glass. The electrons ejected 
from the second photocathode are then focused by another 
electrostatic lens on a second phospher screen to form the 
final image. The arrangement is called a cascaded image 
converter. The glass sandwich in the center is the intensify- 
ing membrane. At present, the main disadvantage of such 
tubes lies in the low resolving power introduced by graininess 
and nonuniformity of the phosphor screens. The gain is not 
high (practically, of the order of about five), and the useful 
field size is smali. These tubes have been used, however, to 


photograph the planet Mars at Lowell Observatory in times of 
the order of 5}, sec, in order to avoid atmospheric smearing. 

Recently progress has been made in the transmitting of in- 
frared images. One of the ideas that appears to be leading to 
successful results centers around the use of an infrared sensi- 
tive material for conversion of the infrared image into a 
temperature matrix which can be scanned like a television 


Fig. 7 Photos of the solar spectrum taken with rocket instruments designed at the Upper Air Laboratory, University 

a ei a of Colorado. The top two regions include the near ultraviolet. The emission cores of the MgII absorption lines 

near 2800 A are easily seen. The third section covers the range from about 2000 to 1000 A; its most prominent 
feature is the hydrogen Lyman-alpha line at 1216 A on the left and numerous other emission lines 


camera tube. Lateral heat conduction is reduced in such « 
converter by the use of extremely thin films such as aluminum 
oxide films about 40 molecules thick developed at Westing- 


house Research Laboratories. Another sensitive material 
that can be used as a target in an errr ision camers 
is gold-doped silicon (28). fe 

Solar Rocket Spectroscopy _ 


The astrophysicist is interested in devising a model of the 
solar atmosphere in order to account for the observation 
which astronomers have made. Another of his problems ha 
to do with the interior of the sun; theories of the latter are tie 
in with the atmosphere of the sun in the sense that tl 
atmosphere represents a kind of boundary condition for tl 
interior, and that the nature of the exchange of energy betwec 
the interior and the atmosphere is a critical factor in reviewin : 
both problems. 

Until the time of rocket observations the sun’s “ees \ 
was limited to about 2800 A on the short wave length side | 
oxygen and ozone absorption and to about 2.5 w on the lon : 
wave length side by water vapor absorption. The ultraviol: t 
region has been extended considerably since rocket instru - 
ments were first flown in 1946. 

The ultraviolet spectrum of the sun may be divided ini» 
three regions: The normal extension of the visible continuu 
with its characteristic Fraunhofer lines to about 1500 
the line emission region which begins around 1900 A an | 
extends well into the extreme ultraviolet region at 84 
and finally the soft x-ray region from 100 to 1 A or les:. 
Sections of the solar spectrum photographed from rockets aie 
shown in Fig. 7 

The use of the expression “x-ray emission” from the sun his 
often led to a misunderstanding; one cannot speak of radi:- 
tion in a given wave length range as being x-ray radiation 
unless the origin of this radiation is typical of that for x-ray 
namely inner shell transitions or electron decelerations. 
the radiation originates as a result of optical transitions in a 
highly ionized atom, that is, as a result of outer electron 
jumps, then the radiation might be correctly referred to as 
extreme ultraviolet radiation, even though it may have 
identical wave length values with x-rays, and even 
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———- SOLAR SPECTRUM RANGE 
Fig. 8 A diagram of the solar spectrum with emphasis on the ultraviolet region. The divisions indicated are ar- 


(He II 303.6A LINE) 


bitrarily made for convenience. The scale is not linear 


tiough no distinction exists between the two radiations, 
once created, insofar as their proton behavior is concerned. 
Tig. 8 is a diagrammatical representation of the solar spectrum 
vith some of the more prominent ultraviolet regions indicated. 

In the past 15 years, astrophysicists have changed some of 
their views regarding the origin of solar radiation in the 
various spectral regions. One may still speak of the photo- 
sphere, chromosphere and corona, but the detailed concepts 
have changed. The photosphere is the visible surface bound- 
ary of the sun and is actually several scores of miles deep. 
Its kinetic temperature is about 6000 K, and the radiation 
emitted corresponds roughly to black body radiation for this 
temperature, being of the nature of a continuum from the 
near infrared to the near ultraviolet. It has been fairly well 
established that most of the solar continuum is caused by 
transitions involving the ion H~. However, the part below 
1800 A may originate (29) as a consequence of the presence of 
appreciable amounts of molecular hydrogen just above the 
photosphere where temperatures drop to 4500 K. Collisions 
of hydrogen molecules with electrons can excite the former to 
levels where transitions to a virtual level involving dissocia- 
tion is possible. The energies of these transitions are con- 
tinuous, and a continuum results. 

The photosphere is also the region of sunspots and the small 
granules, or poxlike structures. The latter have only re- 
cently been photographed (30), under high resolution, from a 
balloon. The granules are as small in some cases as 180 
miles across. They might represent the appearance of 
vertical columns of gas flowing upward from beneath the 
photosphere at speeds of a mile or two per second, and persist- 
ing for only a few minutes. 

The temperature of a shallow strip of the chromosphere 
immediately above the photosphere is less than the tempera- 
ture of the latter by about 1500 K. Some of the lines in the 
Fraunhofer absorption spectrum are formed when light from 
the photosphere passes through this region (sometimes called 
the reversing layer). Higher in the chromosphere the 
temperature rises irregularly again to the order of 30,000 
K before it blends into the corona or outer atmosphere of the 
sun. Many of the emission lines of the solar spectrum in the 
ultraviolet (Fig. 7) originate in the chromosphere. 

Conditions are such that thermodynamic equilibrium 
cannot necessarily be assumed. The process of line emission 
is very complex; a knowledge of many parameters is necessary 
for the formulation of the problem. Ionized atoms are 
common, especially in the upper regions where temperatures 
are high. The problem is further complicated by the fact 
that in some areas, often of large extent, called active 
regions, the chromosphere is considerably disturbed. These 
areas (plages or flocculi) are frequently, though not always, 
associated with the sunspot groups, and they follow the 
sunspot cycle. Only when the sun is photographed in light 
of the wave length corresponding to a transition in a given 
atom are the plage areas distinctly visible. Photographs of 
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the sun’s disk in the radiation of hydrogen Lyman-alpha 
(1216 A) have been successfully obtained from rockets by the 
use of both lithium fluoride and grating optics (25, 26). 
In these photographs, the strong preference for 1216-A 
emission in plage areas is revealed. Sometimes a very bright 
but relatively small chromospheric area called a flare develops 
near a sunspot region. The flare is especially bright in the 
light of the hydrogen Balmer-alpha line. It is invisible 
unless viewed in nearly monochromatic light. Flares may 
last from a few minutes to several hours. They are often 
accompanied by disturbances of various kinds in Earth’s 
upper atmosphere. Photon counters have proved that there 
is considerable enhancement of solar x-ray radiation in a 
flare. Probably enhanced ultraviolet radiation occurs as 
well. Flares occur in the chromosphere of the sun; their 
nature is not yet well understood. 

The solar continuum extends from the visible into the 
infrared. At sea level no observations have been made 
beyond about 2.5 wu because of strong absorption caused by 
water vapor in Earth’s atmosphere beginning at this region 
(31). However, solar radiation is again found in the micro- 
wave region at wave lengths around 1 mm. It is interesting 
to note that the effective black body temperature of the 
infrared continuum of the sun reaches a minimum at 
about 2.5 u; beyond this point the effective temperature 
rises. 

Efforts to measure infrared radiation at higher altitudes 
from airplanes and balloons have been successful (32-34). 
In most cases the resolving power has been quite low, but 
results clearly show the diminishing effects of absorption by 
water vapor and carbon dioxide with increase in altitude. 
There is very little trace of water vapor above about 46,000 
ft, but bands of carbon dioxide persist to much greater 
heights. In these instruments a grating is used to effect 
dispersion. Lead sulfide and lead telluride cells have been 
employed. No measurements have been made beyond 
about 8 uw to date. A large aperture high resolution infrared 
grating spectrometer is scheduled (by NRL and Nat. Bur. 
Standards) for a balloon flight in the near future for study of 
the solar infrared radiation between 18 and 38 p. 

In Fig. 9 is shown the infrared spectrum of the sun from the 
visible to beyond 13 yu (the deflection is not directly propor- 
tional to the energy). The large absorption due to water 
vapor, carbon dioxide, ozone and other constituents of 
Earth’s atmosphere is clearly shown. 


Origin of Coronal Lines 


In the corona of the sun the atoms are highly ionized and 
are greatly outnumbered by the free electrons with which 
they make frequent collisions. The electron temperature is 
between 1 and 2 million deg. The radiation field in the 
corona is mainly determined by the relatively low energy 
photons from the cool photosphere and chromosphere below. 
Thermodynamic equilibrium does not prevail; i.e., collisional 
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processes and radiation processes do not balance within 
themselves. The equations which describe a given phenom- 
enon in the corona are based on equilibrium of a statistical 
nature that is concerned with the balance of the various kinds 
of processes. 

At any one time there is a certain number of atoms of a 
given species in a particular ionized state. Besides abun- 
dance, temperature is the main factor in determining this 
number. Computing 1/2kT for coronal temperatures, one 
concludes that there are few atomic ions with ionization 
potentials less than 200 ev. The coronal emission lines there- 
fore probably are caused by ions that have ionization poten- 
tials greater than this value. 

The mechanism for the emission of a given line in most cases 
is as follows (29). An atomic ion collides with an electron and 
is excited to a higher state. Since for ions the ground state is 
considerably below the first excited state, the chances are that 
the outer electron of the ion will be moved only to the first 
states; higher states are much less apt to be reached by the 
outer electron. The astrophysicist must compute a cross 
section for such a mechanism which he calls a collisional 
excitation cross section. This cross section will of course 
depend upon the energy of the colliding electron so that at 
this point a knowledge of the distribution of the electron 
velocities in the corona is important. When the electron 


reaches the first upper state it almost immediately returns 
to the ground state with the emission of a photon whose 


frequency is given by the usual relationship E = hv, where F 
is the excitation energy for the state. If this transition of 
the electron to the lower state does not occur fast enough, 
the ion may suffer another collision before the electron has 
chance to return, in which case the electron may be raised to 
still a higher state, or may return to the ground state without 
the emission of a photon, if the excess energy is carried off by 
the ion and the colliding particle (radiationless collision). 
However, when the probabilities of the several events ar 
computed, the chances are overwhelmingly in favor of th« 
electron’s returning to the ground state before anything 
happens to the ion, so that the photon is emitted. Such 
photons produce the resonance emission line of the ion. 

It may sometimes happen that the collision of an ion witl 
another particle (electron) excites the outer electron of th 
ion toa so-called metastable level. Transitions of the electro: 
to the ground level are said to be forbidden in such cases 
Actually this means that the electron takes much more tim: 
than usual for returning to the ground level, so that an inte: 
vening collision similar to those described previously is mor: 
likely to occur. Even so, for the conditions which prevail ir 
the corona, the transition from a metastable level is in most 
instances effected, giving rise to a photon corresponding to 
so-called forbidden line in the spectrum. The coronal emis 
sion lines which have for some time been observed in thi 
visible region of the spectrum are all forbidden lines in this 
sense, except that here the transitions are between uppe! 
levels lying close to each other, rather than between a highe1 
state and the ground state. 

This mechanism accounts nicely for most of the permitted 
and forbidden coronal emission lines. However, there is 
another way by which emission may be achieved that is very 
important in relation to certain ions (29). A case in point 
is the oxygen atom with five electrons removed, namely 
OVI. This ion has an ionization potential of 138 ev, hence 
should be almost completely ionized in the corona, so that 
very few of them would be available to cause emission. 
Nevertheless, when one examines the solar ultraviolet spec- 
trum one finds that several lines due to O VI are present, 
and that these represent many more than just the resonance 
line (22). The explanation of these unexpected results 
centers around the high ionization potential of O VII (over 
700 ev). Apparently most of the oxygen atoms in the corona 
should theoretically be in the ionized state O VII, because 


there is not enough energy in the free electrons to ionize them 
further; whereas, the ionization of O VI to O VII is so easily 
~ effected because of the relatively small energy required, that 
- little oxygen should remain in ionized states lower than O 
VII. The O VI lines are emitted after an O VII ion captures a 
_ free electron (recombination) and when the electron cascades 
to the ground state. The O VI ion then makes a collision 
and loses an electron to become O VII again. Such a process 
can be theoretically analyzed only after recombination coeffici- 
ents for the ion have been computed. 


me 
a It can be shown that the coronal gases also emit a con- 
sia tinuum in the soft x-ray region, although most of the emission 
occurs in the coronal lines. The continuum is expected to 
a : be that of a gray body at a temperature of between 1 and 2 
million deg. 
_ An Imaginary Space Solar Observatory 
° 
8 We will here describe the instrumentation which might be 
s| present in a stabilized space solar observatory equipped with 


suitable power and pointing devices, and with telemetering 
and television transmission equipment. The basic principles 
involved for each instrument have been discussed elsewhere in 
this paper. 

Instrument No. 1 is a photon counter with windows so 
chosen that x-rays of wave lengths less than 20 A are meas- 
ured. The gas is neon with some ethyl formate, and the 

- window and aluminum are of such thickness as to give good 
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Fig. 9 The infrared spectrum of the sun from the visible to be- 
yond 13 uw. The diagram was constructed by Gates of the Na- 
tional Bureau of Standards. The deflection is not proportional 
to the intensity, and the scale changes at various places. The 
strong telluric absorption due to bands of HO, CO, and O; is 
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transmission in the 10-A region (21). The upper wave length 
cutoff is at 25 A. In this x-ray region solar radiation is 
usually small, variable or absent. However, at times of solar 
flares the intensity increases remarkably. This latter effect 
is what is expected to be observed by this space observatory 
instrument. 

Instrument No. 2 is a grazing incidence spectrograph de- 
signed to cover the region of the solar spectrum between 
20 and 900 A. It is similar to the one of Fig. 3 except that in 
front of the slit a prefocusing grating is used at grazing in- 
cidence for selecting a wave length range to fall on the slit in 
order to reduce scattered visible light along the Rowland 
cirele and to separate orders (35). The detector is not a film 
but a photomultiplier tube with its glass window coated with 
the phosphor ZnO: Zn which has a relatively high quantum 
efficiency for exciting the radiation in the far ultraviolet (36). 
The tube automatically sweeps the Rowland circle. Trans- 
mitted signals give the position of the tube (wave length) and 
the response when solar radiation of sufficient intensity is pres- 
ent. The prefocusing grating must be rotated and translated 
4s the detector tube moves in order to focus the appropriate 
wave length region on the slit of the spectrograph. In this way 
the spectrum may be observed between 20 and 900 A. Of 
interest would be changes in intensity of solar ultraviolet 
lines as a function of time and the correlation of this fact with 
solar activity. The measurement of the absolute intensity 
of solar ultraviolet radiation as a function of wave length is no 
less an important problem. Were the photomultiplier in the 
spectrograph calibrated, the latter measurement would also be 
possible (37). 

Instrument No. 3 would be a spectrograph of the normal- 
incidence type with optics coated with magnesium fluoride 
on aluminum to insure high reflectivity in the ultraviolet. 
The spectral region of interest would be the hydrogen Lyman 
series and the continuum beyond the series limit. One can 
use a photoelectric-type detector, set behind a movable slit 
which sweeps along the Rowland circle. The device makes 
use of the internal photoelectric effect (38, 39) of tungsten, 
which serves as the cathode. There is negligible response to 
light of wave lengths longer than about 1300 A so that 
scattered visible light in the instrument is no problem. The 
photoelectrons are directed by a magnetic field into a Bendix 
manufactured secondary emission photomultiplier, whose 
output is ultimately telemetered to Earth. Because hydrogen 
is by far the most abundant element in the solar atmosphere, 
and because its behavior there has been the subject of many 
an astrophysical theory, the data obtained by the instrument 
would be of great significance. For the lines and the con- 
tinuum absolute intensities, relative intensities and variation 
of intensity with time are examples of important data that 
may readily be obtained. 

Instrument No. 4 is a high resolution echelle spectrograph 
similar to the one shown in Fig. 6. Because angles of in- 
cidence are small, special optical coatings such as zine sulfide 
or magnesium fluoride to increase reflectivity are needed. 
The observations to be made are the measurement of the 
intensity of hydrogen Lyman-alpha line at 1216 A and of its 
profile under high resolution. The first grating forms a 
stigmatic image of the sun in 1216-A radiation on the slit. 
The final high resolution line image, much dispersed, is formed 
near the entrance slit where it is scanned by a narrow slit 
located in front of a photomultiplier coated with a suitable 
phosphor such as sodium salicylate. The output of the tube 
is amplified and telemetered back to Earth along with the 
signals locating the sweep-slit so as to fix the wave length of 
the part of the line pro“le being examined. In this way the 
line profile, wave length vs. relative intensity, may be re- 
constructed. With care in calibration even absolute intensity 
can be measured. The results are significant in astrophysical 
problems. Continuous observations of the absorption center 
of the hydrogen Lyman-alpha line, Fig. 4, ean yield, provided 
resolution is high enough, critical data on the physical condi- 
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tions which prevail in the chromosphere, or on the nature of 
hydrogen gas clouds in space between Earth and the sun. 
If the focusing grating were repeatedly made to move slightly 
so as to scan the sun by relocating the image on the entrance 
slit, the variation of the profile and intensity of the line over 
the disk of the sun may be studied. 

Instrument No. 5 consists of two parts, one a camera for 
imaging the sun in the visible, the other a camera for imaging 
the sun in the radiation of the hydrogen Lyman-alpha line 
(1216 A). In the first case, direct television techniques may 
be employed for sending the image to Earth; in the second, 
an image amplifier similar to one described previously must be 
used before transmission to Earth. The optical design of the 
1216-A monochromatic camera has already been discussed. 
It is one utilizing but two gratings. These cameras will 
enable an Earth observer to follow important changes of solar 
activity (such as development of flares or large spot regions) 
and, provided electronic communication between Earth and 
observatory is possible, to monitor instruments Nos. 2, 3 and 4 
so that radiation from pertinent areas of the sun’s disk enters 
the instruments. In any event these pictures may be useful 
in identifying the parts of the sun’s disk being scanned by 
each of the instruments and in indicating which regions are 
active. 

Instrument No. 6 is an infrared camera. A high speed, 
grating-type infrared monochromatic camera, similar to the 
ultraviolet one of instrument No. 5, forms an image of the 
sun in infrared radiation from the center of an absorption 
region of one of the bands of a solar compound. The image 
falls on a gold-doped silicon target and is televised to Earth. 
The picture will give information about the distribution of 
the molecular compound over the solar disk, and about the 
temperature of various regions of the solar atmosphere 
where the molecule can exist. 
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Stationary supersonic nozzle flow of a reacting gas mixture, in thermodynamic equilibrium or in 
states between equilibrium and frozen flow, was investigated experimentally and analytically. 
With this fully determined flow system, it was possible to check various criteria to predict departures 
from equilibrium. A method proposed by Bray for Lighthill’s ideal dissociating gas was extended 
to our real gas system, and the results showed agreement with the experiments. | 
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IGH speed flows with large changes of temperature pro- 

vide interesting problems in thermodynamics, since the 
usual notion of treating the medium as a thermally and 
salorically perfect gas becomes invalid. However, so long as 
it is permissible to assume that thermodynamic equilibrium 
prevails everywhere, certain flow fields are amenable to cal- 
culation. The problem becomes more difficult if it is found 
that departures from thermodynamic equilibrium occur in 
the flow, and rate processes must also be taken into account 
to compute parameters of practical interest, such as the heat 
transfer to a blunt body in hypersonic flow, the temperature 
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of a flame front, the thrust of a rocket nozzle or the electro- 
magnetic properties of a shock wave region. Such calcula- 
tions depend on the knowledge of the relaxation times of all 
transformation processes involved, and they are usually diffi- 
cult or impossible to perform. It is at this point that one 
wishes to resort to a “criterion’”’ with which to predict depar- 
tures from equilibrium in a given flow situation. Such a 
criterion may require an estimate of the dynamics of the 
process and some understanding of its mechanism, in addition 
to an estimate of the flow parameters of the limiting condi- 
tions for frozen or equilibrium flows. Also, numerical answers 
should be rapidly available. The latter requirement is im- 
portant, since, with the information outlined above, the 
equations of motion would, in principle, be approximately 
solved, and we expect the application of equilibrium criteria 
to effect a simpler computation. Finally, the simplification 
of involved rate processes inherent in a workable criterion may 
iced to a deeper physical understanding of flow situations. 
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We wish to discuss the rapid expansion of a gas or gas mix- 
ture from a reservoir with the medium undergoing simul- 
taneous dissociation and recombination reactions. It will be 
assumed that the flow is in thermodynamic equilibrium in the 
reservoir, and that, furthermore, all external and internal 
degrees of freedom of the participating species remain in equi- 
librium during the expansion. Practical applications of this 
process are the flow in a rocket nozzle, the flow around a 
blunt body moving at hypersonic speeds, and the expansion 
of a gas mixture in a hypersonic wind tunnel from a reservoir 
of high stagnation enthalpy such as that found in plasma jets 
or arc heated hypersonic wind tunnels. 

Many of the chemically important features of such flows 
may be modeled by studying experimentally and analytically 
the stationary flow of a reacting gas mixture of nitrogen 
dioxide NOs and tetroxide NO, carried at low concentration 
in inert nitrogen Ns. This reaction may be written 


7 kp 
N2 + N2O, Ne + [1] 
kr 


where kp and ke are the rate constants for dissociation and 
recombination, respectively. With a suitable experimental 
environment it is possible to produce expansions of this mix- 
ture through a supersonic nozzle, anywhere from near frozen 
flows to flows in chemical equilibrium. Such a system, in- 
cluding the dynamics of the reactions, is fully determined, 
permitting the use of this flow as a testing ground for the 
validity of some of the criteria for thermodynamic equilibrium 
found in the literature. The procedure will be to compute 
the terms in the criteria without resorting to approximations 
and to compare the predictions with the experimental obser-— 
vations. 

These are the criteria for near equilibrium (or near frozen) — 
flow initially proposed by Penner and co-workers (1,2,3)* in 
conjunction with performance calculations for rocket nozzles.4 
The near equilibrium criterion embodies the calculation of a 
small temperature difference representative of a small differ- 
ence between the actual concentration of the reactants and 
that concentration predicted for chemical equilibrium. The 
method was used by Logan (4) who applied it to hypersonic 
flow around blunt bodies and in a nozzle. A different ap- 
proach to the prediction of departures from chemical equilib- 
rium is based on a comparison of the numerical value of the 
rate of change of degree of dissociation with distance in a 
nozzle with the magnitude of the dissociation rate term in the 
equation for the reaction mechanism, with both terms evalu- 
ated for chemical equilibrium. This method was proposed 
by Bray (5,6) who applied Lighthill’s theory of the ‘ideal 
dissociating”’ gas (7), with Freeman’s extension to nonequi- 
librium conditions (8), to the flow in a hypersonic nozzle. 
Bray’s suggestion can easily be extended to real gases, and it 
will be shown that this method provides remarkable insight 
into the flow in expansions.° 


Analytical Methods and Experiments 


Equations of Motion for Frozen and Equilibrium Flow 


The equations of motion for the stationary, one-dimen- 
sional flow of the mixture of Equation [1], neglecting friction, 
heat conduction and diffusion, have been derived elsewhere 
We define 


(10). The equations will be merely listed here. 
the degree of dissociation of the reactants by 
Q = WNo,/(WNO, + WN20,) [2] 


where the mass flow fraction of the ith species is given by w; = 


3 Numbers in parentheses indicate References at end of paper. 

‘ The references given are those most readily available rather 
than those that trace the historical development. 

5 Another scheme proposed, but not tested with our system, 
was recently given by Rudin (9). 
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m;:/m, and symbols without subscripts refer to the gas mix- 
ture. Throughout this work we will discuss only flows where 
dissociation in the nozzle supply (subscript 0) is complete and 
(oxz0,)o = 0. From Equation [2] 


wno, (NO2) 


(WNOs)o N¢ 
where 
) = concentration (mole/em*) 
F = frozen flow without composition change 


We now write the continuity and momentum equations in 
the customary forms 


Apu =m 


—pudu = dp 


where 
u = velocity 
p = density 
p = pressure 
A = cross-sectional area of the nozzle 


With ap = 1, the integrated energy equation is given by 


Wn An, + (WNoO,)o hNo. — Wn,01 Ahn,o, + 


(Wn )o + [6] 


where 


me” = enthalpy per unit mass of the 7th species 
= the heat of dissociation of per unit mass of 


Generally, hyo, is greater than hy,o,; and we deal with an 
exothermic reaction if NO. recombines to N2Ox during the 
nozzle expansion. Finally, if we assume all participating 
species to be thermally (not calorically) perfect gases, the 
equation of state assumes the form 


1 — (wno,) (WNO.)o — WN,O. 
MN» MNO: HN20, 
where 
‘a 
7 = temperature 
R = the universal gas constant 


Mi = molecular weight of the 7th species 


For given A(z) we may now compute from Equations 
[4-7] all flow parameters for the frozen flow for which 
(wn.0,)0 = wnx.o, = 0 in Equations [6 and 7]. The needed 
thermodynamic properties of the reactants are well known 
from the sources critically reviewed by Gray and Yoffe (11). 
The opposite limiting expansion process of flow in chemical 
equilibrium exhibits the greatest possible change of composi- 
tion during the expansion. Its properties for known A(x) may 
be computed by adding to Equations [4-7] the law of mass 


action 
(WNO-)o” 
Kp = 2—— p| —— — 2wnoy)o + wry, [8] 
MNO, e 
where 
wal 
yp = equilibrium constant in units of pressure ; 
K; juilil tant ts of 
Bw = the molecular weight of the mixture vs 
e = chemical equilibrium 


In all succeeding examples, values for the two limiting flows 
will be shown with concentrations computed from (7) = 


Wip/ Mi. 
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FLOW DIRECTION 


-10 > 10 15 
x, cm 
Fig. 1 Nozzle contour for exit Mach number of 1.8 in nitrogen 


Calculation of Relaxation Time 


Frozen and equilibrium flow may also be viewed kinetically 
as exhibiting zero or infinite reaction rates, respectively. In 
addition to the equilibrium thermodynamic properties of the 
species of the gas mixture of Equation [1], the rate constants 
are known (10,12,13); the reaction mechanism of Equation 
[1] may be expressed by 


= 2kp(N204)(N2) — [9] 
T 


Equation [9] describes a second-order dissociation reaction 
opposed by a third-order recombination process. In the 
small temperature interval of interest, the recombination rate 
constant is independent of temperature whereas the dissocia- 
tion rate constant follows an Arrhenius temperature depend- 
ence given by kp ~ exp (—B/T). It ean be shown (10) by 
integration of Equation [9] that the reaction time 7 elapsing 
for a change of reactant concentration from (NO.); to (NO2) 
at constant temperature is given by 


in NO») = C] [=(Ke/2) + Cl 
[—(K./2) — 2(NOz) + C] [=(K-/2) — 2(NO2) C] 


where 


= [(K.2/4) + 4KA(N.Oy), + {11] 


and A, is the equilibrium constant in units of concentration 
given by K. = kp/kr = K,/RT. As long as the NO, con- 
centration at the end of a chosen time interval is not too close 
to that prevailing in chemical equilibrium, either of the terms 
on the right-hand side of Equation [9] may be neglected, 
depending on whether dissociation or recombination pre- 
dominates. In the latter case Equation [9] becomes 


(ONO:/dt)p = —2ke( NO2)2(N2) = k’(NO2)? [12] 


List of experimental conditions 


Table 1 
Departure from 
chemical 
To; equilibrium Number of A 
K xr, em experiments (see Eq. [27]) Remarks 
400 4.5-5 5 1 least reliable owing to low 
4,54 2.21 408 >15 chemical 
F 5.00 2.16 402 >15 1 equilibrium throughout 


with &’ constant at constant temperature. The integration 


of Equation [12] leads to 


2(N2)(NOs)kr L(NOs) 


_ which is a simpler formula for computing the reaction time. 


Apparatus 


The experimental equipment described elsewhere (10) con- 
sists primarily of a continuously operating supersonic wind 
tunnel of 2 X 3.2-em exit cross section. The reactants are 
fed at a known rate into a dry nitrogen flow with supply 

_ pressure pp of about 2 atm and supply temperature 7 of about 
400K. The reactant concentration @ in mole per cent in the 
supply is defined by 


= [14] 


with values up to 5 per cent employed practically. For these 
typical supply conditions, a» is equal to 1 as postulated in the 
previous analysis. The expansion takes place through a 
converging-diverging, water-cooled nozzle designed to give 
continuously varying flow parameters. The nozzle outline is 
shown to scale in Fig. 1 with the sonic throat located at x = 0 
in reference to the succeeding figures. The flow in this nozzle 
is one-dimensional (14), and the function A(z) entering the 
equations of motion is empirically determined from a nitrogen 
‘alibration giving the effective cross-sectional area. The 
Mach number at the nozzle exit is about 1.8 with the NO. 
largely recombined to N2O, for chemical equilibrium. The 
local NO» concentration can be determined directly from an 
absorption measurement with light of about 4300 A (10,14), 
since NsO; and N» are transparent. The N.O; concentra- 
tion finally may be found from 


(N204) = (1/2)[(NO2)r — (NO2)] [15] 


At the higher initial reactant concentrations, static pressure 
measurements in the flow in conjunction with the known area 
and Equations [4-7] permit the determination of all other 
flow parameters (10). The conditions of selected experi- 
ments, identified alphabetically, are given in Table 1. Two 
contrasting examples are shown in Figs. 2 and 3, where the 
experimentally found dimensionless NO» concentration in the 
nozzle is plotted with the curves predicted from the previous 
analysis for frozen and equilibrium flow. In Fig. 2 it is seen 
that the flow of experiment C follows the curve for equilibrium 
up to x ~ 3 em, beyond which the observed flow departs 
markedly from chemical equilibrium and exhibits reactant 
concentrations at the nozzle exit about halfway between the 
two limits. (These stations of departure from equilibrium 
are also noted in Table 1.) In contrast to this behavior, 
experiment F on Fig. 3 is « flow which is in chemical equi- 
librium throughout. Departure from chemical equilibrium is 
more clearly seen if degree of dissociation is plotted rather 
than concentration, as shown again for experiment C in Fig. 4 
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Fig. 2. Dimensionless NO, concentration in the nozzle with the 
flow observed to depart from chemical equilibrium 


‘Temperature Lag or Characteristic Time Criterion _ 


Prediction of Departure 
From Chemical Equilibrium 


Penner and co-workers (1,2) considered a residence time of 
the chemical species in a given temperature interval in the 
expanding flow of a rocket nozzle. Comparing composition 
computed for the end state of this residence time with that 
predicted for chemical equilibrium, it is then possible to deter- 
mine if equilibrium was maintained or not. In principle this 
method may be applied to the entire expansion in one step, 
or one may subdivide the flow at will (1,4). The treatment 
may be applied to a reaction of the type of Equation [1], which 
is formally similar to the reassociation of atoms in the presence 
of an inert third molecule, and Altman and Carter’s dis- 
cussions of the method (3) as extended to our own reaction 
will be followed. We obtain for chemical equilibrium, from 
Equation [9], by setting the left-hand side equal to zero 

kp ( Js 


K. = [16] 
kr (N Ya de 


Equilibrium in a flow will be maintained as long as the reac- 


tant concentration follows . 
> 

din Ke d\n(NOs)? — d [17] 

dt dt dt : 


where the quantities d In (_)/dt refer to the total change of 
concentration of the reacting species with time. In a nozzle 
expansion the concentration change results from two separate 
processes: The flow itself and the chemical reaction. Each 
term on the right-hand side of Equation [17] may then be 
split according to these two processes, and the total rate of 
change of concentrations may be expressed directly. In this 
expression one can replace the concentration ratio at equi- 
librium of Equation [16], corresponding to K, (7), by an as- 
sumed equilibrium constant K,'(T’). In chemical equi- 
librium T’ = 7, whereas for flows whose composition is near 
equilibrium ir — T\, called the temperature lag, is a small 
positive number for expansions. Near equilibrium K,’ 

about equal to K,; and expanding the equilibrium constant in 
a Taylor series about T, retaining the first two terms and 
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Fig. 3 Dimensionless NO, concentration in the nozzle with the 
flow in chemical equilibrium 


| 


making a number of substitutions and approximations from 
thermodynamic relations, the temperature lag for our reaction 
may finally be estimated (3) by 


(1 — CrT/Ahy,o,) (—dT/dt) _ ( dT 

kp(No) [20NO.) + ) 
If the temperature lag computed for a given expansion is 
small, equilibrium is expected. If (7" T) is large, the 
flow has departed from equilibrium; however, since Equation 
[18] is valid only near equilibrium, no quantitative statement 
on the extent of the departure may be made. The expression 
given by 7 has the dimension of time, and (—dT’ dé), the in- 
verse residence time per unit temperature interval, may also 
be viewed as the cooling rate in the nozzle. This cooling rate 
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Fig. 4 Degree of dissociation in the nozzle with the flow ob- 
served to depart from chemical equilibrium 
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in our experiments varies from 10° to 10° K per sec; in smaller 
rocket nozzles it may go above 10’ K per sec. It is even 
higher in hypersonic nozzles, and in each instance the maxi- 
mum value occurs near the sonic throat. 

The criterion of Equation [18] may be viewed as the product 
of a reciprocal characteristic flow time and a maximum reac- 
tion time Tmax. Logan (4), following a further suggestion by 
Penner (2), assumes the flow to be in equilibrium as long as 
the temperature lag is much smaller than the local tempera- 
ture. Equation [18] may, then, be written as a condition for 
equilibrium 


where the choice of the value 10~* is admittedly arbitrary. 


dT 


x < 10-3 19 
dt | Tmax S [ ] 
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Fig. 5 Temperature lag criterion computed for several 
experiments 
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Fig. 6 Characteristic flow time and relaxation time in the range 
of all experiments ecw 
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For 7max Logan computes a relaxation time; the reaction time 
implied in Equation [18] or that given by Equations [10] or 
[13] may be used. The initial reactant concentration in all 
instances ought to be chosen consistently to maximize the 
reaction time. 

In our examples, we find at each nozzle station (NO»)r > 
(NO.) >(NO,),. Ata given position, therefore, we decide 
to calculate the maximum possible reaction time from Equa- 
tion [10], choosing (NO2); = (NO). Using the integral o! 
the complete equation, we must determine a final reactant 
concentration (NO2) = (NO2), > (NO:)., to be reached at 
constant temperature, by 


(NOz), = (NOs). + (1/e)[(NOz)r — (NO2)e] [20 


in keeping with the usual definition of relaxation times. I 
other words, we define the maximum reaction time in th 
nozzle as the time it would take the reaction at T = constan' 
to proceed from (NO:)r, the maximum possible reactan 
concentration, to (NO»), defined by Equation [20]. Th 
actual choice of Tmax does not appear to affect the criterion o 
Equation [19] significantly; the choice may merely shift th: 
arbitrary value of 10~* as a limit for equilibrium.® 
The temperature lag computed in this manner from Equa 
tions [19, 10, 11 and 20] is shown in Fig. 5. The departur: 
from equilibrium (Table 1) is indicated by the change fron 
solid to broken curves. This break is not seen at a consistent 
(T’ — T) value. Fig. 6 shows the ranges of the character- 
istic flow times and the maximum reaction times for al 
experiments with the minimum values occurring near thi 
throat. Finally, in Fig. 7 we see values for the criterio: 
Equation [19] plotted directly. The departures from equi 
librium indicated are all appreciably above the value 107°; 
however they are, in addition, not consistent with any par- 
ticular value. It is believed that another interpretation o/ 
initial reactant concentration in Equation [18], or any othe: 
consistent choice of Tmax, would only lead to a shift of th« 
6 This point is borne out by a numerical example for experiment 
Batz =5cm. The following are the values of the concentra- 
tions in units of 10-7 mole/em*: (NOz)r = 1.78, (NO2) experi- 
ment = 1.38, (NOs) = 1.32, and (NO2). = 1.05. The time of 
reaction, in microsec, from (NOz2)r to (NOz)r is 150 (Eq. [101), 
120 (Eq. [13]), 120 (Eq. [18]). The half life (Eq. [13]) is 340 
starting at (NO.)r. To go from (NO) experiment to (NO2), 
(Eq. [13]) gives 140, and all these numbers are of about the same 
magnitude. 
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Fig. 7 Characteristic time criterion for several experiments 
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curves in Figs. 5-7 without appreciably affecting their basic 

trends. The same finding is given by Bray (6), who _— 
that the approximate solution of Equation [18] leads to — 
appreciable errors if compared with this exact calculations on : 
hypersonic nozzle flow. It is in this context that Logan’s — 
predictions (4) are to be viewed with care. Penner heel 
implies (1) that such detailed results as hoped for by Logan, 7 
and in our calculations ought not to be expected. He sug- 
gests’ that higher order terms in the Taylor series must be | 
retained to obtain more realistic results. However, some 
important correct qualitative facts emerge from Figs. 5-7. 
Ali experiments show minimum temperature lag values ini- 
tially in the expansion. This finding agrees with the trend 
o! the curves for N recombination, O recombination, etc., in 
hypersonic nozzles given by Logan (4). Therefore, qualita- 
tively, we find in the regions of the largest temperature gra- 
dients (or shortest characteristic flow times) near the nozzle 
throat, that the flow is closer to equilibrium than at the exit. 
This behavior can also be seen in the complete calculations on 
hydrogen recombination in a rocket nozzle by Krieger (15). 
Another complete calculation by Heims (16), who determined 


the recombination of atomic oxygen in hypersonic stagnation 
: ig. issociati recombination terms for experimen 

point flow (or nozzle flow), shows the same result. It is ane’ recombination terms | 

therefore clear that, in recombination reactions in expansions, yy 

the initial regions of near sonic flows are those in which t 

thermodynamic equilibrium is likely to prevail, as shown 

experimentally or theoretically in many instances (5,6,10,14, 

15,16). Unfortunately, the opposite statement is occasion- Equation [22] is positive in expansions, since a decreases with 


ally found in the literature. 


flow becomes frozen. 


gen and oxygen through a hypersonic nozzle using the model 
of Lighthill’s ideal dissociating gas (7) with Freeman’s rate 
parameter (8) chosen to adapt numerical results to the cur- 
rent range of uncertainty of the actual rate constants. 
Bray finds that the hypersonic nozzle flow may be qualita- 
tively described by observing three flow regions. Initially 
the flow is in equilibrium and follows the prediction of the 
law of mass action. Next, there is a short transitional regime 
in which the flow departs markedly from equilibrium, and 
finally the flow is essentially frozen. The choice of a range of 
rate constants affects only the axial location of the departure 
from equilibrium, not the described three stages of the ex- 
pansion. Based on the finding of the apparent ‘sudden 
freezing”’ of the flow, computed with the full equations, Bray 
proposes a simpler approximate solution in which the flow is 
initially in equilibrium; the point of ‘sudden freezing’ in 
the nozzle is determined by a criterion based on the rate 
equation, and from then on we have frozen flow. This be- 
havior is remarkably similar to that observed in our system. 
as shown on Fig. 4. Bray (6) suggests that it may be possible 
to apply his criterion to a real gas, and we will do this with 
our system. The rate Equation [9] with Equation [3] 
gives 


inequality 


ment B. 


dax (Nz) 
( dt ), (NOs) [kr(NOz) p( N20 [21] 


gible values. 


where the subscript A signifies that at a given nozzle station of N.O, is practi rally equal to R, which, however, is also a 
for which the equation is evaluated numerically, the area is small value, a result which leads to practically frozen flow. 
to be held constant so (NO.) » = constant and d(NO.)»/dt = 0. With the inequalities of Equations [23 and 24] borne out by 
In this manner conditions at various stations may be compared the experiment, it is logical to ask for the approximate condi- 
directly. With this restriction in mind, Equation [21] may tion of departure from equilibrium or sudden freezing. Since 


be shortened by 


3 


proaching zero. 


as a condition for near frozen flow. 
the expansion, we may compute the terms D and R from the 
complete equations, and they are plotted in Fig. 8 for experi- 
Here we see the number of NO. molecules pro- 
duced by the recombination term or destroyed by the dissoci- 
ation term per unit volume and unit time. 
chemical equilibrium on the scale shown, R is about equal to D, 
and both are large, and the net production R — D is small. 
After departure from equilibrium, D drops rapidly to negli- 
In fact, at the nozzle exit the net production 


> D 


x, com 


These statements applied 


—da/di< D 


as a requirement for the flow to be near equilibrium. 
down the nozzle, with concentrations and temperature de- 
creasing, both R and D fall rapidly. 
tion region or region of departure from equilibrium, D ap- 
proaches zero faster than R, owing to the strong temperature 
dependence of the dissociation rate constant. 
yond the departure from equilibrium, the flow may be viewed 
as nearly frozen with the left-hand side of Equation [22] ap- 
In almost frozen flow, D will be negligibly 
small; the net rate of the reaction is about equal to the re- 
combination term R, and we obtain from Equation [22] the 


15 


time (or distance), and the rate approaches zero if the nozzle 
Initially, in the nozzle we have equi- 
Criterion Based on Terms in the Rate Equation halt dar librium; the net recombination rate is small; R and D are large 
he fl terms, and they are about equal. 
3ray (5,6) calculates the flow of initially dissociated nitro- to Equation [22] lead to the inequality 
23] 
Farther 


Beginning in the transi- 


Finally, be- 


[24] 


o verify this view of 


Initially in 


the net production term is much smaller than D in near 
equilibrum, and much larger than D for near frozen flow, 


99 
sg atl at Bray assumes that somewhere in the region of departure the 
where R and D are the recombination and dissociation rate two terms may be of the same order of magnitude, setting 
functions implied in Equation [21]. The left-hand side of bill il [25] 
where it is hoped that A is a constant of order unity. The 


7 Private communication from Dr. 8. Penner. 
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great practical advantage of this criterion for predicting de- 
parture from equilibrium is that both sides of Equation [25] 
may be evaluated at equilibrium since the flow is in this 
condition up to the departure point. From Equation [22] we 
may then write 


— (da/dt), = AD. ~ AR, [26] 


since D ~ R near equilibrium; and R, is simpler to compute 
since kg is constant for our reaction. With da/dt = uda/dz, 
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EXPERIMENT B 
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EXPERIMENT F 
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ox 
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Fig. 9 Criterion based on comparison of dissociation rate term 
with rate of change of degree of dissociation for three 
experiments 
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Fig. 10 Test section length needed in terms of nozzle length L to 
return a esate equilibrium from nozzle exit conditions for 


and with Equations [21 and 26], the criterion for sudden freez- 
ing for our reaction finally becomes 


da, 2(N2)(NO2),? 
= Aker 27 
dx u(NOs)r (37) 


The results of three calculations with this criterion are shown 
in Fig. 9. They were obtained by inserting in Equation [27 
values from a complete calculation for chemical equilibrium 
with Equations [4-8], which include the variation of (NO.), 
with z. Also shown in Fig. 9 are the experimentally observed 
departure points from chemical equilibrium (Table 1). The 
agreement of prediction and observation for experiment B i: 
remarkable (i.e., A ~ 1), and for experiment F departure is 
neither observed nor predicted. From Equation [27], 
may be expressed by choosing the value of the other terms at 
the observed departure point. The results are 0.42 < A < 
1.25; i.e., the departure ordinarily occurs upstream of the 
point of prediction with A < 1. 

Of particular interest to hypersonic tunnel design is the fact 
that once the flow departs from equilibrium, it will not readily, 
return to this state in the region of the test section. In par- 
ticular, it is fallacious to assume that mere tunnel size wil! 
materially affect this behavior. The additional length of test 
section, in units of nozzle length L, needed at constant flow 
conditions to return the flow to near equilibrium is shown it 
Fig. 10. Calculations with Equations [10 and 13] are given 
and the limits of the simplified calculation are clearly seen 
More than 10 nozzle lengths must be added to the nozzle to 
return to near equilibrium ; beeause of the very low rate 
of recombination, this i ‘is not a ie result. 


The fully determined nozzle flow system with two reacting 
species lends itself, indeed, to a check of approximate methods 
to determine departures from thermodynamic equilibrium 
The important point is that these checks can be made without 
resorting to approximations. It was found that methods 
based on characteristic flow and relaxation times, as developed 
for rocket nozzles and also applied to hypersonic flow prob- 
lems, do not furnish sufficient detail for predictions of points 
of departure from equilibrium. They are valuable, at most, 
to generally indicate in which regions of the flow such de- 
partures are likely. A method recently proposed by Bray 
using Lighthill’s ideal dissociating gas is based on an estimate 
of the magnitude of the dissociation rate term computed for 
chemical equilibrium with the change of degree of dissoci- 
ation, and it was shown to give satisfactory results when ex- 
tended to a real gas flow. Aside from a reasonable prediction 
of points of departure from thermodynamic equilibrium, the 
method also sheds light on the physical understanding of 
flows in expansions with reactions that are primarily those 
of three-body recombinations. Contrary to statements 
occasionally found in the literature, it can be shown that 
chemical equilibrium in expansions is most likely to prevail 
in the regions of high temperature gradients initially in the 
expansion. It appears that this fact is independent of the 
area ratio of the nozzle or the substance used, within wide 
ranges covering rocket nozzles, as well as hypersonic wind 
tunnel nozzles. It is further found that once departure 
from equilibrium takes place, no practical test section length 
of a hypersonic nozzle will bring about the return to equilib- 
rium, and that the mere invocation of large-scale equipment 
occasionally proposed is not a guarantee of sustained chemi- 
cal equilibrium. 


Conclusions 
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Nomenclature 


A = cross-sectional area of nozzle 


A,B = constants 
Penner, 8. S., ‘‘Chemistry Problems in Jet Propulsion,’’ Pergamon 


Cc = defined by Equation [11] Press; London; 1957. 
h = enthalpy per unit mass 3 Altman, D. and Carter, J. M., ‘“Nonequilibrium Effects,” in ‘‘High 
Ah = heat of dissociation per unit mass ie ei yo ype re Propulsion,” Princeton University Press, 
N.J., 1956, vol. II, part 1-B-4. 
ky dissoc tation rate constant 4 Logan, J. G., ‘‘Relaxation Phenomena in Hypersonic Aerodynamics,” 
kp ee recombination rate constant TAS Preprint no. 728, Jan. 1957. 
ie K» = equilibrium constants 5 Bray, K. N. C., “Departure from Dissociation Equilibrium in a 
L = nozzle length Nozzle,” Aeron. Research Council, London, rep. no. 19,983, 
March 1958. 
me = mass flow 6 Bray, K. N. C., “Atomic Recombination in a Hypersonic Wind- 
p = pressure Tunnel Nozzle,” J. Fluid Mech., July 1959, vol. 6, pp. 1-32. 
? = universal gas constant r 7 Lighthill, M. J., ““Dynamies of a Dissociating Gas, Part 1: Equi- 
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x = distance on nozzle center line 9 Rudin, M., “Criteria for Thermodynamic Equilibrium in Gas Flow,” 
= degree of dissociation Phys. Fluids, Sept.—Oct. 1958, vol. 1, pp. 384-392. 
lar 10 Wegener, P. P., ‘Supersonic Nozzle Flow with a Reacting Gas Mix- 
MOIE war weigh ture,” Phys. Fluids, May-June 1959, vol. 2, pp. 264-275. 
p = density 11 Gray, P. and Yoffe, A. S., ‘‘The Reactivity and Structure of Nitrogen 
= reacts . ati j gzle « 12 Carrington, T. anc avidson, N., ‘Shoe Yaves in Chemical Ki- 
the nozzle supply netics: The Rate ot Dissociation of N2Os,"’ J. Phys. Chem., April 1953, vol. 
= mass fraction, 1.€., = 57, pp. 418-427. 
) = concentration 13. Wegener, P. P., ‘Measurement of Rate Constants of Fast Reactions 
. in a Supersonic Nozzle,”’ J. Chem. Phys., April 1958, vol. 28, pp. 724-725. 
A M - 14. Wegener, P. P., Marte, J. E. and Thiele, C., ‘‘Experiments on Chemi- 
Subscripts ; cal Kinetics in a Supersonic Nozzle,’ J. Aeron. Sci., March 1958, vol. 25, pp. 
205-206. 
) = supply conditions ’ 15 Krieger, F. J., “Chemical Kinetics and Rocket Nozzle Design,” Jour- 
= chemical equilibrium NAL OF THE AMERICAN Rocket Socrery, Nov. 1951, vol. 21, pp. 179-185. 
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Transport and Thermodynamic 
Properties in a Hypersonic Laminar 
Boundary Layer 
Part 2 Applications 


General Electric Co. 
Philadelphia, Pa. 


® 


An analysis is presented in which the transport and thermodynamic properties of dissociated air, 


which were presented in Part 1 of this paper, and the Curtiss, Hirschfelder and Bird method of pre- 
dicting the thermal conductivity of a chemically reacting gas. are utilized in the determination of 
stagnation point heat transfer to a hypersonic vehicle. New theoretical results are obtained for 
stagnation point heat transfer, and these are compared with earlier theoretical results, and with 
available experimental data. It is found that at both high and low altitudes the authors’ results 
are larger than the theoretical values predicted by Sibulkin and Lees, but smaller than those of 
Romig, and Fay and Riddell. It is also found that the authors’ results are consistently lower than 
the mean value of the shock tube data of Rose and Stark, and Vitale et al. However, it is noted 
that the scatter in the experimental data on heat transfer is of the same order as the variation in the 
different theoretical predictions, 


HE PROBLEM of accurately predicting the stagnation air. [See Part 1 (1).4]> The second part consists of the 
point heat transfer rate for hypersonic flow is treated in solution of the boundary layer equations utilizing these prop- 
two parts. The first is the determination of the transport and erties. Considered in this manner, a given solution of the 
thermodynamic properties of high temperature dissociated boundary layer equations is as “exact’’ as the particular set 
_——__—_—_—___— ‘ Numbers in parentheses indicate References at end of paper. 
Received Feb. 18, 1959, and revised Feb. 10, 1960. 5 Regarding Part 1 of this paper, note the following corrections: 
1This paper is a revised version of General Electric Co., Equation [3] 


MSVD, TIS no. R59SD306, Feb. 1959. m + ft? 
? Manager, High Altitude Aerodynamics, Space Sciences Labora- [Diezh = 1.171 x 10-679" [ 2mym2 | a2? sec q 
tory. Member ARS. 
Member ARS. Table 1: Species No, «/k = 79.8 K. 


’ Physical Chemist, Space Sciences Laboratory. 
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of gas properties utilized in the analysis is ‘‘accurate.” ai 

Earlier studies of stagnation point heat transfer by Sibulkin 
(2), Lees (3), Romig (4) and Fay and Riddell (5) have in- 
corporated assumptions such as those given in Equations 
[12]. In this paper all of the transport and thermodynamic 
properties of dissociated air were evaluated during the solu- 
tion of the differential equations, as a function of the local 
gas composition and temperature within the boundary layer. 

Since there is no simple way of distinguishing between the 
separate effects of variable Prandtl number, Lewis number, 
viscosity, thermal conductivity, etc., our results are compared 
with earlier work on the basis of the gross stagnation point 
heat transfer to the surface, for given environmental condi- 


tions. 
le 


Boundary Layer Calculations 


The gas in a hypersonic boundary layer is sae a multi- 
component mixture of the products of dissociated air; it is 
convenient in first approximation to treat the gas as a binary 
mixture of “air atoms” and “air molecules” (3), whose ther- 
modynamic and transport properties are then averaged in 
some appropriate manner. 

Consider the laminar flow around a blunt axially sym- 
metric body; the boundary layer equations which represent 
the conservation laws for a binary gaseous mixture may be 
written (6) as follows: 


Global conservation of mass 
(pure) > (per) =0 


Conservation of species 7 


( Oc; re) (oo OC; 4 + 
u— = — ; [2 
ox oy oy ‘> T oy wi 


Conservation of momentum neglecting centrifugal effects — 


Or oy Or oy 


Conservation of energy neglecting the Dufour effect 


(. or dp (2) (x or 
Ox oy Ox oy oy oy 


> | (oo Dit | or 


4) 
i 
It is noted that the diffusion coefficients have been defined 
so that the binary diffusion coefficient 9;; and the thermal 
diffusion coefficient D;” satisfy the following equation for 
the diffusion flux of the ith species 


—pD; Ve; — D;TV [5] 


piV; = 


Although the effects of thermal diffusion have been largely 
ignored, a hypersonic boundary layer exhibits steep thermal 
gradients and large differences in the molecular masses of the 
gaseous species present, which tend to promote the occurrence 
of this phenomenon. Consequently, one of the authors has 
usually included thermal diffusion in his analyses (7). It is 
also noted that a special case exists when the Lewis number is 
unity and thermal diffusion is neglected. 

Upon multiplying Equation [2] by h; and taking the sum- 
mation over 7, the heat release term due to chemical reac- 
tions can be eliminated from Equation [4] to yield 


oh oh Op ou \? 


Oc; 
abe hi — + Dit 


sits 


This form does not contain the chemical source term. 
However, nothing has yet been gained since the com- 
position gradient Oc;/Oy still appears in the last term on the 
right-hand side. Thus, the diffusion equation is still required. 
However, noting that h = 2 cih:, the (0/Oy)(KOT/dy) term 


1 
may be eliminated from Equation [6] 


oh oh re) ou \? oh 


” oy oy Oy 


Pa) K\ oc; D;? oT 
2 - T dy 7] 


— the Prandtl and Lewis numbers, the latter be- 
comes 


ox oy oy oy \Pr ody 
ro) oc; Di? oT 
oy Ls Pr oy 


If the Lewis number is unity and thermal diffusion can be 
neglected, the energy equation becomes 


oh ied oh OP 19] 
or oy ox oy oy \Pr oy 
This form constitutes an enormous simplification over 
Equations [4 or 6], since now the diffusion equation is com- 
pletely uncoupled, and the effects of the chemical heat release 
— term Baa ; are included implicitly without recourse to a 


[ chemical rate equation. 
Note that a very similar type of equation can also be written 
in terms of the stagnation enthalpy H = h + u?/2 


P oy oy Pr oy 
vir 


dy Le dy \Pr 


so that, upon assuming the Lewis number to be unity and 
neglecting thermal diffusion, the latter becomes 


p Oe oy Oy Pr oy 


oy LPr oy (11) 


Therefore, the utilization of Equations [9 or 11] uncouples 

the diffusion equation even for the case of a nonequilibrium 

flow without resorting to any specific constraint on the chemi- 

cal source term. That is, the ¥ wih; term is accounted for 
t 

implicitly; i.e., it adjusts itself. 

The difficulty is that now the fluid properties cannot be 
properly evaluated, since there is no direct step from the 
enthalpy to the properties which depend on both the tem- 
perature and the composition. Of course, in first approxima- 
tion, one could use the properties of a gas which is in thermo- 
dynamic equilibrium between the same limits of free stream 
and wall enthalpy, but this approach may not prove to be 
satisfactory. It appears that for a precise treatment of 
boundary layers containing chemically reacting components, 
even if the )> wih; term is eliminated from the energy equa- 


(Le — y= 


110] 


t 
tion, the diffusion equation which contains the w; term is still 
required. The diffusion equation must, in fact, be utilized 
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simultaneously with the other equations in order to evaluate 
the composition and temperature, so that the transport and 
thermodynamic properties can be properly calculated. It is 
noted that in an approximate solution, the use of assumptions 
such as 


pu = constant 
Pr = constant ae 

Le = unity te 

€ = constant 


w 
‘liminates the need for evaluation of the transport and ther- 
modynamic properties, but conclusions must be drawn with 
‘onsideration of these assumptions. 

In the analysis, in which dissociated air was treated as a 
binary mixture of air atoms and air molecules, the properties 
of the gaseous mixture were determined during the course of 
the solution of the boundary layer equations in terms of the 
local gas composition and temperature. 

For the binary mixture approximation, the following arbi- 
trary assumptions were utilized 


= 15 my = 30 aan. = 
=> Don, 
= (1/2(mo + 
= 
= Cpair atom 
= Cpair molecule 

[13] 


and a mass weighted heat of formation was taken as 


LA IN IN 


In order to include the effects of thermal diffusion, which 
may be important in the nonequilibrium boundary layer, the 
thermal diffusion ratio defined by (8) 


kr 


was introduced, where, if ky > 0, component 1 diffuses from 
the hot to the cold region and vice versa. It is also noted that 
D," = —D,.", where the first approximation to the thermal 
diffusion ratio is given by (8) 
XxX X, — 
(6C\2* — 5) [16] 


kr= 


Xx + Yr 


The quantities ¥,x, ¥x and [Kj2]; have been previously in- 
troduced in Part 1, and S® and S® are unique functions of 
m, m2, Ky, Ky, Ke and Ay*, where and are func- 
tions of the reduced temperature 7'* (8). 

Upon applying the Mangler and Dorodnitsyn transforma- 
tions (3) 


Pelle i [17] 
7=- rody 7 
7 
x 
f, [18] 


to Equations [1-4], one obtains the following set of ordinary © 
differential equations: 
Conservation of atoms 


Wa 
2pldu,/dz), 
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Momentum 


Energy 


z l 


Le? 
4 (16 CAL An 


These have been utilized recently in a study of catalytic ef- 
fects in the nonequilibrium hypersonic laminar boundary 


layer (10). 
These equations must satisfy the following boundary con- 


ditions. At the surface 


3 
= fry = = 0 4 [22] 
f(0) = fo = —vy |2— ( - =0 [23] 
De 
00) = 0 = T./T. [24] 
ZA(O) = Zaw = CAw/CAc [25] 


2s 


(assuming no slip or mass transfer at the surface). At the 
outer edge of the boundary layer 


lim fy = lim @ = lim Z4 = 1.0 [26] 


Once these equations are solved, the resulting heat transfer 
to the surface which is given by 
of 
i w 


becomes 


w 
= = [220+ 
) 


V 2Pubu(due/dx 8 
Le 


in terms of the transformed variables. The minus sign indi- 
cates that the heat is transferred to the surface, and the right- 
hand side of Equation [28] includes the contributions of con- 
duction, ordinary diffusion and thermal diffusion. In these 
calculations, the authors have assumed that the stagnation 
point velocity gradient is Newtonian, i. ¢ : 


dx Rep Pe 


Hence the effect of body size appears as the square root of 


the nose radius V/ Rp. 

It is also noted that the results which are obtained in the 
transformed plane may be transformed back into the physical 
plane by means of the inverse Mangler-Dorodnitsyn trans- 


formation 
w 2 w 1 € alti € 
y= Pu ( ] Pe dy [30] 
Pe \dx /s 0 p 


Discussion of Results 


Before discussing the general results which were obtained 
on an IBM 704 digital computer, attention will be given to a 
typical boundary layer solution obtained during hypersonic 
flight. Ata Mach number of 20 and an altitude of 120,000 ft, 
the stagnation point conditions of pressure, temperature and 
atom concentration at the outer edge of the boundary layer 
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(Ufan)n + ffin + > = [20] 
p 
b 
| 
ma 
DAM 
MA 
Mu 
Cpa 
D Le CA,Z An + - On + fcaZan + a 


are estimated as p, = 2.4 atm, 7, = 12,400 R, ca, = 0.482, 
based on free stream conditions predicted by the 1956 ARDC 
atmospheric model (11). 

A surface temperature of 800 R was selected, and, since the 
atom recombination rates are still not known within several 
hundred per cent, calculations were performed for a wide range 
of atom recombination processes. Thus, in Fig. 1 are shown 
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Fig. 1 Boundary layer profiles at Mach 20, altitude 100,000 ft, 
and a surface temperature of 800 R 
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the nondimensional boundary layer profiles obtained at the 
stagnation point when it is assumed that the boundary layer 
is either in “local thermodynamic equilibrium,” or when atom 
recombination is “frozen” in the gas phase and allowed to oc- 
cur on four surfaces having differing finite catalytic efficiencies 
(7). 

Other results of the computations are shown in Figs. 2-5 
where the gas properties are plotted as a function of the local 
temperature in the boundary layer. It is seen that when it is 
assumed that the boundary layer is in thermodynamic equi- 
librium, the contribution of the atomic species to the behavior 
of the properties of the gaseous mixture is a minimum. Since 
the wall has a finite catalytic efficiency, when the reactions are 
frozen the atom concentration builds up on the surface, and 
the atomic species contribute in larger measure to the prop- 
erties of the mixture within the boundary layer. 

Some of the more interesting results indicate that for low 
wall temperatures, ¢, is not constant (Fig. 2), and, as Fay and 
Riddell have shown (5), the product of density and viscosity 
is not constant (Fig. 3). The Prandtl and Lewis numbers 
were not assumed constant as in previous analyses, but were 
found to be variable and, furthermore, exhibited gradients 
which affect the temperature and concentration gradients at 
the wall. (See Fig. 4 in which the ratio of Prandtl to Lewis 
number, i.e., the Schmidt number, is shown.) Further, the 
effects of thermal diffusion were included; the thermal dif- 
fusion coefficient is shown in Fig. 5. Note that since the ther- 
mal diffusion coefficient is negative, the atoms, which are the 
lighter particles, diffuse thermally from a low temperature 
region to a high temperature region. Thus, thermal diffusion 
tends to oppose concentration diffusion and reduces the dif- 
fusion flux of atoms toward the wall. 

If one defines the nominal boundary layer thickness in 
terms of the point at which the velocity in the boundary 
layer achieves 99.9 per cent of its free stream value, then one 
obtains the results given in Fig. 6. Since the boundary layer 
thickness is directly related to the density of the gas in the 
boundary layer, it is seen that the nominal thickness of the 
boundary layer increases with decreasing free stream density. 
Thus, the prime effects are due to changing the altitude. 
There are slight Mach number and surface temperature effects 
as well. Higher flight speeds V. tend to compress the gas 
and decrease the boundary layer thickness. On the other 
hand, higher surface temperatures 7, act to decrease the 
mean density of the gas in the boundary layer and hence 
increase the thickness of the boundary layer. 
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Fig. 3 Normalized product of density and viscosity of binary 
mixture air in boundary layer 
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Fig. 4 Variation of Schmidt number of binary mixture air in 
boundary layer 
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Fig. 5 Thermal diffusion coefficient of binary mixture air in 
boundary layer 


Of greater interest are the results given in Figs. 7 and 8. 4 
These curves show the dependence of hypersonic stagnation 
point heat transfer Q,, (in the absence of mass transfer) upon 
environmental conditions. The two surface temperatures 
selected are typical of heat sink vehicles and ablating vehicles. 
As is readily seen, the heat transfer increases with flight speed 
and free stream density and decreases with surface tempera- 
ture. The calculations for heat transfer were performed as- 
suming that the gas was in dissociation equilibrium at the 
wall. It is noted that when one makes this assumption, the 
results are virtually independent of the rates of homogeneous 
atom recombination (3,5,10). 

As stated earlier, these results were obtained for variable 
Prandtl, Lewis and thermal Lewis numbers. Thus, one can- 
not readily compare our results with analyses in which arbi- 
trary and constant values were selected, except on an absolute 
basis. This is done in the next section. 

At this point it must be noted that the heat transfer results 
presented in Figs. 7 and 8 are based on the use of the frozen 
thermal conductivity recently derived by Curtiss, Hirsch- 
felder and Bird (9) and mentioned by the authors in Part 1 
(1). For purposes of comparison, the simpler relationship 
between the thermal conductivity of a binary mixture and the 
thermal conductivity of the pure species given in Part 1 was 
also employed. It was found that use of the latter resulted 
in heat transfer rates which were of the order of 10 per cent 
lower than those based on the more rigorous Curtiss, Hirsch- 
felder and Bird expression. This may be attributed to the 
fact that the relation, as derived in (8) for the thermal con- 


Aprit 1960 


250 

200 

> 12* Vay FT /SEC. 

2 Tw*800 °R 
—— 3500 °R 

100 

/ 


Oo 005 O15 020 .025 030 035 


Wh, FT? 


Fig. 6 Nomina! thickness of the boundary layer at the stagnation 


point; 
ALTITUDE FT. 
L @ 20,000 0) 
9+ @® 30,000 
40,000 
« 8- Tw = B00°R 
50,000 +. — 
O aL © 60,000 4 
w 7 
SF 80,000 Z 
or @ 90,000 | 
° @ 100,000 
x 
3 
a 
L a® 


14 16 18 20 22 24 


FT./SEC. 


Fig. 7 Heat transfer vs. flight speed—low altitude 
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ductivity of “binary mixture’’ dissociated air, yields values 
which are of the order of 10 to 20 per cent larger than those 
obtained using the approximate expression presented in Part 
1. Hence, although the simpler method may be considered 
adequate, the more rigorous equation should be used wherever 
possible. 


Comparison With Earlier Work 


In some of the earlier theoretical studies (2-5) approxi- 
mate solutions, exact numerical solutions and empirical cor- 
relation formulas have been offered which represent relation- 
ships between the stagnation point heat transfer rate, the en- 
vironmental conditions and the gas properties. Thus, the 
heat transfer rate can be calculated at a given flight speed and 
altitude if one has available properties of the standard atmos- 
phere (11), normal shock tables (12,13,14) and transport and 
thermodynamic property data (1,15). 

Sibulkin (2) suggested a heat transfer rate equation based 
on an extrapolation of a low speed boundary layer solution in 
the form 


Q = 0.763(Pr)*(h. — hw) 


Lees (3) published a formula which is similar to Sibulkin’s 
result, in which the Lewis number was taken equal to unity 
and hy << h, 


Q = 0.71(Pr)~?*h. [32] 


Romig (4) derived a “reference enthalpy’ method for find- 
ing the heat transfer rate. Beginning with a form of Sibul- 
kin’s equation, she obtained an expression for the heat trans- 
fer rate in terms of flight Mach number, nose radius and am- 
bient pressure, but wall temperature effects were not included 

Q = [33] 

Fay and Riddell (5) have published a complex correlation 

formula based on extensive numerical solutions for which con- 


stant values were assumed for the Prandtl and Lewis num- 
bers 


whw 1 1/2 
Q = 0.763(Pr)*(he — hw) ) 
Pefe dx}. 


+ — 1) [34] 


The magnitude of the heat transfer rate given by each of 
these four equations clearly depends on the gas properties 
which are utilized in the correlation formulas. For a fixed 
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value of the Prandtl number, the form of these four equations 
is such that the viscosity of the gas at the outer edge of the 
boundary layer enters the final result explicitly to the 0.5 
power in Equations [31 and 32], implicitly to the 0.5 power in 
Equation [33] [see (4)], and explicitly to the 0.4 power in 
Equation [34]. Thus, the uncertainty in the heat transfer 

rate is 50 per cent of the uncertainty in the value of y, ap- 

pearing in the correlation equations given by Sibulkin, Lees 

Bar Romig, and only 40 per cent of the uncertainty in the 
value of u, appearing in Fay and Riddell’s equation. 

It is therefore interesting to compare the values of the vis- 
cosity of dissociated air predicted by a number of different 
theories. In Fig. 9, we have shown the value of viscosity pre- 
dicted by the Sutherland law (curve 4), the value of pix 
(curves 7 and 8) predicted by the Buddenberg and Wilke 
relation (Part 1, Eq. [5]), and the value of the viscosity 
(curves 2,3 and 5) predicted by the recent work of Hansen (15). 
For reference purposes, the viscosity of pure atomic oxygen 
(curve 6) and pure molecular oxygen (curve 1) given in Part | 
are also shown. It is noted that curves 2 and 7 were calcu- 
lated for the environmental conditions simulated by Rose and 
Stark (16), who kept the shock tube initial pressure p, con- 
stant. Curves 3 and 8 were calculated for the environmental 
conditions simulated by Vitale et al. (17) who kept the stagna- 
tion pressure p, constant. 

It is seen that for the particular choice of the rigid sphere 
model for the atomic species and the Lennard-Jones model for 
the molecular species, the predictions of the Buddenberg and 
Wilke expression utilized by the authors diverge significantly 
from Hansen’s predictions with increasing temperature. 
(Compare curve 8 with curve 3, and 7 with 2.) Also note 
that the values of pmix fall below the Sutherland law values, 
whereas Hansen’s results lie above. 

The effects of dissociation upon the viscosity coefficient are 
particularly pronounced where the stagnation pressure is low 
and the concentration of atomic species is high, for a given 
range of hypersonic flight speeds. Therefore, the deviation 
of the viscosity coefficient from the Sutherland law (which 
neglects dissociation) appears more strongly in high altitude 
shock tube simulation, e.g., p, = 0.034 atm, than in low altitude 
simulation, e. g., p: = 0.01316 atm (see Fig. 9). Thus, the 
use of different viscosity laws should lead to larger differences 
in the theoretical heat transfer rates at high altitude than at 
low altitude. 

In the new theoretical calculations presented here, the 
transport data given in Part 1 were employed. However, for 
purposes of comparison, the authors’ transport data were 
not used in the empirical correlation formulas [31-34]. 
Instead, the Prandtl and Lewis numbers were taken as con- 
stants, and the Sutherland viscosity law was utilized. Thus, 
although the Prandtl number is actually a variable, it was 
taken as a constant (Pr = 0.71) in the correlation equations. 
Although Buddenberg and Wilke’s equation is probably more 
representative of the behavior of a gaseous mixture at high 
temperatures than is Sutherland’s law, the latter was em- 
ployed by Fay and Riddell. Hence, for consistency with 
their numerical calculations, it was also employed here ex- 
plicitly in Equations [31, 32 and 34], and implicitly in 
Equation [33], since Romig utilized Sutherland’s law in the 
derivation of Equation [33]. 

The Lewis number, which is actually a variable, was taken 
as unity by Lees in the derivation of Equation [32] and as a 
constant by Fay and Riddell. Therefore, for comparison, the 
Lewis number was taken as Le = 1.0 and Le = 1.4 in Equa- 
tion [34]. 

Shock tube data on stagnation point heat transfer are avail- 
able only for a limited number of altitude conditions, e.g., 
Rose and Stark, and Vitale et al. Consequently, representa- 
tive high and low altitude conditions were investigated, and 
the theoretical predictions were evaluated for the environ- 
mental conditions which were simulated in the shock tube ex- 
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Fig. 10 Comparison of theoretical and experimental heat 
transfer—low altitude 


Fig. 10 shows a comparison between theory and experiment 
it simulated altitudes in the range 65,000 to 75,000 ft. It is 
seen that the predictions of Fay and Riddell, Romig, and Scala 
and Baulknight fall within the scattered experimental data 
points of Rose and Stark, whereas the values predicted by 
Sibulkin and Lees fall below. 

Fig. 11 shows a comparison between theoretical and experi- 
mental results at simulated altitudes in the range 200,000 to 
240,000 ft. Here again it is seen that the differences in the 
theoretical predictions are of the same order as the scatter in 
the experimental data. It is noted, however, that at the higher 
altitudes, as anticipated, there is a somewhat more pro- 
nounced spread in the different theoretical predictions than at 
lower altitudes. 

It is to be noted that one could, of course, by judicious 
selection of the transport properties to be inserted in the 
various correlation equations, appreciably change the magni- 
tude of the individual theoretical predictions of heat transfer, 
a posteriori, until “agreement”? with available experimental 
data is obtained. By way of illustration, Hansen’s viscosity 
data could be used to increase the magnitude of the theoretical 
heat transfer rates given by Lees’ and by Sibulkin’s formulas, 
which would then bring their predictions closer to the avail- 
able experimental data. Since the same technique can be 
employed with any correlation formula to obtain a better 
empirical fit, it is doubtful whether this type of approach can 
be used to rationalize the superiority of a particular theoretical 
treatment upon which the correlation formula is based. One 
still requires independent knowledge of the transport prop- 
erties of the gaseous system. miei 


Conclusions ga 

The use of the approximate expression for the thermal con- 
ductivity of a chemically reacting binary mixture gas discussed 
in Part 1 [i.e., Eq. [6] of (1)] may be considered adequate 
for heat transfer calculations, since comparison with the more 
rigorous expression derived in (8) shows that the error in heat 
transfer introduced by this relation is only of the order of 10 
per cent. However, the Curtiss, Hirschfelder and Bird expres- 
sion should be used, wherever possible. 

When our estimates of the viscosity of dissociated air are 
compared with those predicted by Sutherland’s law and by 
Hansen (15), it is found that the results diverge significantly 
with incre¢ sing temperature. At low temperatures where 
dissociation is negligible, one finds, as anticipated, that the 
three methods yield practically identical results. However, 
as the level of dissociation rises, Hansen’s predictions yield 
values which are of the order of 25 per cent larger than 
Sutherland’s law and 50 per cent larger than our results for 
the cases treated here (see Fig. 9). 
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Fig. 11 Comparison of theoretical and experimental 
transfer—high altitude 


Since the uncertainty in the theoretical prediction of heat 
transfer is of the order of 50 per cent of the uncertainty in 
the coefficient of viscosity, the marked disagreement between 
the different theoretical predictions of viscosity at high tem- 
perature is in itself sufficient justification for an experimental 
program leading to the determination of high temperature 
transport properties. 

If one accepts the premise that the binary mixture approxi- 
mation adequately represents the behavior of dissociated air, 
then the accuracy of the new theoretical heat transfer calcula- 
tions is not limited by the mathematical model utilized here. 

A comparison between the various theoretical predictions of 
stagnation point heat transfer indicates a maximum deviation 
of approximately +20 per cent from a mean value. The ex- 
perimental heat transfer data indicates that the scatter ob- 
tained at high gas enthalpy levels (V. >20,000 fps) is a 
maximum of approximately +25 per cent from a mean value. 
On the basis of the relationship between heat transfer and 
transport properties discussed earlier, it may therefore be con- 
cluded that the uncertainty in the experimental data is 
actually at least equal to, if not greater than, the uncertainty 
in the high temperature transport properties. This situation 
at present precludes the possibility of inferring the inter- 
molecular potential function from gross heat transfer data 
taken at high gas enthalpy levels 

Thus, although the rigid sphere and the Lennard-Jones 6:12 
models are undoubtedly only zeroth order approximations for 
the high temperature behavior of the atomic and molecular 
species, little evidence exists that any other arbitrary choice of 
high temperature intermolecular potential function can be ex- 
clusively recommended at present. 

Finally, a comment should be made with regard to the in- 
ternal consistency of the calculations of the transport proper- 
ties of the gaseous mixture. Of the four transport properties 
utilized in boundary layer calculations, three appear exclu- 
sively in the Prandtl, Lewis and thermal Lewis numbers. 
These are the thermal conductivity coefficient A, the ordinary 
diffusion coefficient 2,2. and the thermal diffusion coefficient 
D,’. Thus, the assumption of specific values for the Prandtl 
and Lewis numbers appears to preclude the necessity for in- 
dependently computing these transport coefficients. How- 
ever, since the coefficient of viscosity appears not only in the 
Prandtl number but also in the shear stress term in the 
momentum Equation [3], and in the viscous dissipation term 
in the energy Equation [4], one must decide on a specific 
model for the viscosity law. 

Since any viscosity law is in fact related to a potential en- 
ergy function for the intermolecular interactions, the assump- 
tion of a particular model implies that all other transport 
coefficients have also been implicitly specified. Thus, cal- 
culations which are based on constant values of the Prandtl 
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and Lewis numbers may be subject to internal inconsistency, 
unless it can be shown that the assumption of constant dimen- 
sionless groups very nearly represents the actual behavior of 
the gas. Since the assumption of constant Prandtl and Lewis 
numbers considerably simplifies the numerical calculations, it 
is recommended that these dimensionless groups be evaluated 
at the wall conditions and then held constant across the 
boundary layer. This would constitute a better approxima- 
tion than the assumption of arbitrary constants. 
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Nomenclature 
Cy* = functions of 7'/(«2/k) 


pi/p, mass fraction of species 7 
Cp; = (Oh;/OT)p, specific heat at constant pressure of 


s 


species 
Cp > “frozen’”’ specific heat of mixture 


1 
Di = binary diffusion coefficient 
D;" = thermal diffusion coefficient 
enthalpy of mixture 


h; —~—s = ‘static enthalpy of species 7, including heat of for- 
mation 

Ahy? heat of formation of species 7 
mechanical equivalent of heat 
K; thermal conductivity 
ky = thermal diffusion ratio 
Le = ptpD»/K, ‘frozen’? Lewis number 
l = pu/pwltw, dimensionless 
Le™ = ¢,D,;"/K, frozen” thermal Lewis number 
mi = molecular weight of species 7 
M Mach number 

= molecular weight of mixture 


moles of species 7 per unit volume 
static pressure 

Cou /K, * ‘frozen’”’ Prandtl number 
Schmidt number Pe 
temperature 
T/(«/k), reduced temperature 
x-component of velocity, parallel to surface 
diffusion velocity of species 7 

y-component of velocity, normal to surface 

net mass rate of production of species 7 per unit 
volume per unit time by chemical reaction 


= coordinate system 
Xi = mole fraction of species 7 
Z; = ¢i/C,, referred mass fraction of species 7 
p = density 
T = p Ou/dy, shear stress 
n & = similarity variables 
0 = 17'/T,, temperature ratio 


atomic species 
interaction function for an atom-molec 


ture 

D = dissociation 

e = outer edge of the boundary layer 

i = ith chemical spec ies 

M = molecular species 

t = total 

w = wall 

12 = interaction function for a binary mixture of 
species 1 and species 2 — 

[ h = first approximation 


References 


1 Seala, 8S. M. and Baulknight, C. W., “ 


Transport and Thermodynamic 
Properties in a Hypersonic Laminar Boundary Layer, Part 1 Properties of 


the Pure Species,”” ARS Journat, vol. 29, no. 1, Jan. 1959, pp. 39-45. 

2 Sibulkin, M., ‘‘Heat Transfer Near the Forward Stagnation Point of 
a Body of Revolution,” J. Aeron. Sci., vol. 19, no. 8, Aug. 1952, pp. 570-571. 

3 Lees, L., “Laminar Heat Transfer Over Blunt-Nosed Bodies of 
Revolution at Hypersonic Flight Speeds,’ Jer Propuxsron, vol. 26, no. 4, 
April 1956, pp. 259-269. 

4 Romig, M. F., “Stagnation Point Heat Transfer for Hypersonic 
Flow,” Jer Proputsion, vol. 26, no. 12, Dec. 1956, pp. 1098-1101; Adden- 
dum, vol. cig ~~ 12, Dec. 1957, p. 1255. 

5 Fay, . and Riddell, F. R., “Theory of Stagnation Point Heat 
— in Met iated Air,” J. Aeron. Sci., vol. 25, no. 2, Feb. 1958, pp. 

3-85. 

6 Seala, S. M., ‘‘The Equations of Motion in a Multicomponent Chemi- 
cally Reacting Gas,’ ’ General Electric Co., M.S.V.D., R.M. no. 5, T.I.S 
document no. R588D205, Dec. 1957. 

7 Seala, 8. M., ‘‘Hypersonic Heat Transfer to Catalytic Surfaces,” 
J. Aeron. Sci., vol. 25, no. 4, April 1958, pp. 273-274. 

8 Hirschfelder, J. O., Curtiss, C. F. and Bird, R. B., “ 
of Gases and Liquids,” John Wiley and Sons, N. Y., 1954. 

9 Curtiss, C. F., Hirschfelder, J. O. and Bird. R. B.. * 
Transport Properties,”’ in ‘Transport Properties in Gases,” 
Univ. Press, Jan. 1958, pp. 3-11. 

10 Seala, S. M., “Hypersonic Stagnation Point Heat Transfer to Sur- 
faces Having Finite Catalytic Efficiency,”’ Proc. Third U. 8S. Nat. Congress 
Applied Mechanics, June 1958, pp. 799-806. 

11 Air Foree Cambridge Research Center, ‘ARDC Model Atmosphere 
1956,”’ General Electric Co., M.S.V.D. document no. 568D233. 

12 Logan, J. G., Jr., ‘Normal Shock Calculations in Air,’’ Cornell Aero- 
nautical Lab., Jan. 1957, General Electric Co., Subcontract SDP-1001. 

13. Feldman, ‘Hypersonic Gas Dynamic Charts for Equilibrium 
Air,”’ Aveo Research Lab., Jan. 1957. 

14. Huber, P. W., “Tables and Graphs of Normal Shock Parameters at 
Hypersonic Mach Numbers,’’ NACA, T.N. 4352, Sept. 1958. 

15 Hansen, C. F., ‘Approximations for the Thermodynamic and Trans- 
port Properties of High-Temperature Air,” NACA, T.N. 4150, March 1958. 

16 Rose, P. H. and Stark, W. I., “Stagnation Point Heat Transfer 
Measurements in Dissociated Air,”’ J. Aeron. Sci., vol. 25, no 2, Feb. 1958, 


Molecular Theory 


Theories of Gas 
Northwestern 


17 Vitale, A. E., Kaegi, E. M., Diaconis, N. S. and Warren, W. R., 
“Results from Aerodynamic Studies of Blunt Bodies in Hypersonix Flows of 
Partially Dissociated Air,’’ Heat Transfer and Fluid Mechanics Institute, 
June 1958, pp. 204 215. 


i 
"7 
~ 5 
Subscripts 
* 
A 
: 
i 
n 
rr 
Sc 
T 
<A 
Ve | 
%.. ~ 
i 


Combustion of Highly 


tunnel. Aluminum borohydride, 


silane were tested. 


the tunnel walls. 
studies in larger wind tunnels. 


 igtenlr enka of heat to a supersonic airstream produces 
£&% aerodynamic effects which may have application in 
aireraft propulsion or control. Theoretical studies of hea 
addition to supersonic flow are reported by a number 3 
authors (1-4). Practical applications are considered in 
papers exemplified by (2). A study (2,5,6) was begun 
several years ago at the NASA Lewis Research Center, of 
the feasibility of adding heat to a supersonic airstream in 
in the vicinity of airfoil by combustion, in the absence of any 
disturbances except those caused by the airfoil and by the 
combustion itself. More recently, a number of other in- 
teresting closely related studies have been reported (7,8). 
It therefore seems appropriate to summarize some of the 
NASA work at this time. 

Aluminum borohydride was originally chosen as the fuel 
because it ignites readily in air even under extremely severe 
conditions. However, it presents many problems which de- 
tract from its practicality as an aircraft fuel. Jt decomposes 
slowly during storage, liberating large volumes of hydrogen. 
Its extreme reactivity requires that all parts of the fuel 
system be scrupulously cleaned and purged of air to prevent 
formation of materials which will plug orifices. It therefore 
appeared desirable to search for other materials which could 
be burned in available test facilities, so that further aerody- 
namic and combustion studies might be carried out with 
greater ease of fuel handling. 

The work was extended to include a study of the pos- 
sibility of burning other fuels in a small Mach 2 wind tunnel: 
Pentaborane, trimethylaluminum, diethylaluminum hy- 
dride, trimethylborane, triethylborane, alkylboron hydrides, 
vinylsilane and mixtures of aluminum borohydride with 
JP-4 (9). 

The fuels were injected from the top wall of the wind tunnel. 
Combustion phenomena were observed by means of high 
speed Schlieren and direct motion pictures. In addition, 
the gross nature of the flow field associated with the combus- 
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in Supersonic Airstreams 


The feasibility of adding heat to supersonic airstreams by combustion was studied i ina small wind 
pentaborane, hydrocarbon-aluminum borohydride mixtures, 


ignite or burned only in the diffuser. 


Reactive Fuels} 


EDWARD A. FLETCHER? 
ROBERT G. DORSCH® 
HARRISON ALLEN Jr.‘ 


Lewis Research Center, NASA 
Cleveland, Ohio 


trimethylaluminum, diethylaluminum hydride, alkylboranes, alkylboron hydrides and _ vinyl- 
The first three ignited easily and burned well. 
Trimethylaluminum and diethylaluminum hydride produced 
light and other evidence of heat evolution when water was simultaneously injected into the tunnel. 
Gross effects on flow were studied by observation of shock patterns and water sprays injected from 
Examples are given whic h illustrate the use of these techniques in aerodynamic 


The others either failed to 


tion of aluminum borohydride was studied by observation of 
the behavior of water streams inected into the stream. 

The techniques developed in the combustion work have 
been applied in a series of aerodynamic studies. They in- 
clude studies of the effect of the combustion of aluminum 
borohydrice on the pressure distributions around flat plates 
(10), a body of revolution (11), a supersonic wing in a small 
(1 X1-ft) wind tunnel (12), and various flat plates in a very 
large (10 10-ft) wind tunnel (12). 

This paper describes the techniques used to achieve com- 
bustion in supersonic airstreams with several fuels and some 
of the combtstion and flow phenomena which were ob- 
served in a small wind tunnel with aluminum borohydride. 
It also describes the behavior of some of the fuels which could 
not be ignited under the conditions of the experiment. 
Finally, examples are given of the aerodynamic studies which 
have been made in larger wind tunnels; and some of the 
more important details of the experimental technique used 
to achieve combustion are presented. 


Combustion and Gross Airflow Studies in a Small 
Wind Tunnel 


Apparatus and Procedure 7 


The experimental setup consisted of the wind tunnel, 
photographie equipment, fuel injection system and ignition 
system arranged as shown schematically in Fig. 1a. 


Combustion was studied in a 3.84 X 10-in. supersonic wind 
tunnel having plate glass sides. The windows ran the entire 
length of the tunnel down to the subsonic diffuser (Fig. 1b), 
permitting convenient visual observation and high speed 
(1000-4000 frames/sec) motion picture photography. The 
tunnel was operated over a range of Mach numbers with 
aluminum borohydride. However, most of the work was 
done at Mach 2. Tunnel stagnation pressures ranged from 
37 to 50 in. of mercury. The tunnel air had a nominal dew- 
point of —20F and was preheated to 80-L10F. —— 
test section conditions are listed in Table 1 eed 


Wind tunnel 
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MIGICTURE CAMERA 


(a) Viewed from top of tunnel, showing approximate camera 


positions 


Pressure measurements 


Strain gage differential pressure transducers and NASA 
standard base six capsule differential pressure manometers 
were used to measure static pressure ¢ — along the center- 
line of the top wall of the test section. A detailed descrip- 
tion of this pressure instrumentation is presented in (10). 

A boundary layer study made in this tunnel at Mach 2 
(5) indicates that the boundary layer flow at the point of 
fuel injection is as shown in Fig. 2a. Several stagnation 
pressure surveys made during the present investigation agreed 
with these data. Fig. 2b shows the thickness of the bound- 
ary layer as a function of the distance downstream of the 
throat. At the point of fuel injection, the boundary layer 
is about 0.6-in. thick; however, the subsonic — the 


boundary layer is only about 0.01-in. thick. if 


Fuel injection 


Aluminum borohydride was injected into the airstream by 
two different techniques. Flow rates of the order of 100 
ec per sec (55 gm/sec) were injected for short intervals (3 to 
8 millisec) from capsules which were mounted in the wall of 
the tunnel. The end was merely broken off a pressurized cap- 
sule to permit a 3 to 1 ce slug of material to be injected into 
the stream. 

This mode of fuel injection was unsuitable for aerodynamic 
study and was eventually abandoned in favor of the fuel 
injector illustrated in Fig. 3. This injector was charged with 
5 to 15 ec of fuel through the upper valve either by distilla- 
tion to the filling arm from a conventional high vacuum sys- 
tem or by direct fuel transfer to the injector in the inert 
atmosphere of a dry box. The method used depended on the 
properties of the fuel. The injector was mounted near the 
upstream end of the test section on the centerline, flush with 
the top wall of the tunnel. After steady supersonic flow had 
been established in the tunnel the injector was pressurized with 
helium through the top valve, the lower valve was opened 
remotely, and the liquid fuel was sprayed into the tunnel. 
A considerable variation in the flow rate could be achieved 
by varying the helium pressure and the injector orifice size. 
The intere hangeable orifice caps are shown in Fig. 3. The 
duration of the fuel injections varied from 1 to 3 sec. 

Water or a hydrocarbon fuel (JP-4) was injected during 
some aluminum borohydride runs through 0.028-in. diameter 
holes in the top wall of the tunnel. JP-4 was used in order 
to see if a hydrocarbon flame could be piloted with aluminum 
borohydride. Water streams were used to study the gross 
nature of the flow field in the heated region. The difference 


between the tunnel static pressure and atmospheric pressure 


forced the liquid into the tunnel. 
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Ignition system” 


Aluminum borohydride ejected from capsules usuall: 
ignited spontaneously. It was found desirable to introduc: 
a spark plug to insure prompt ignition in the runs where th 
injector was used because of the lower fuel flow rates (5 
The spark plug was mounted 25.25 in. downstream of th 
fuel injector, as shown in Fig. 1, and produced a 1-joule 
five-spark-per-sec repeating capacitance spark. After ignition, 
the flame, accompanied by a shock wave, quickly traveled 
upstream to the point of injection and remained there until 
the fuel was expe nded. 


. 
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Table 1 Typical test section conditions associated with 
various Mach numbers in the 3.84 X 10-in. supersonic wind 
tunnel 

Linear 
Static Static flow 
Mach pressure, tempera- velocity, Quality 
number n. Hg ture, F fps of flow 
1.5 12 —74 1443s‘ Fair; static 
pressure 


profile ir- 


regular, 
large 
= number of 
shock 
waves 
present 
Good; static 
pressure 


profile less 
= irregular, 
fewer 


shocks 
present 
oxcellent; 
fairly flat 
static 
"pressure 
profile,few 
disturb- 
ances 
2270 ~=Poor; flow 
slightly 
unstable, 
presstire 
profile ir- 
regular 


— 148 1727 


— 260 2082 


—327 
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DISTANCE FROM TOP WALL OF TEST SECTION 


(a) Velocity profile through boundary layer at point of fuel 
injection 


AIR VELOCITY, FT/SEC 
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INJECTOR POSITION 
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- 


BOUNDARY-LAYER THICKNESS, IN. 
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(b) Boundary layer thickness as function of downstream distance 
from throat of tunnel 


Fig. 2. Boundary layer flow in 3.84 x 10-in. supersonic wind tunnel at Mach 2 


tluminum borohydride 


Results 


A summary of observations made at Mach 1.5, 2, 3 and 4 
with slugs of aluminum borohydride from glass capsules 
has been presented in (5). Aluminum borohydride almost 
always ignited and burned vigorously. Flame usually ap- 
peared as a momentary, very intense, bright green flash, and 
was frequently accompanied by a loud bang that could be 
clearly heard above the very high noise level of the tunnel. 
The greater part of the combustion took place within a period 
of about 3 to 8 millisee. Occasionally there were less vigorous 
runs, where it appeared that the main body of the fuel had not 


ignited. In these runs the flame was weak and pale blue. 
Many of the burnings, however, took place with explosive 
violence, frequently completely choking the tunnel and driv- 
ing the flame at very high speed a considerable distance 
upstream of the point of injection, sometimes almost back to 
the tunnel throat. An example of such a run is shown in 
Fig. 4. 

No clearly defined structural details, such as a distinctive 
flame front could be discerned either in the direct or Schlieren 
photographs. In such runs, burning usually started about 1 
millisec after the start of the fuel injection and continued 
quite vigorously for about 8 millisee. By this time the main 
body of the fuel had been used up, but driblets from the 
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Fig. 4 Open shutter photograph of combustion of }-cc slug of aluminum borohydride of Mach 4. Note that flame extends almost 
back to throat of tunnel 


Fig.5 Open shutter photograph of the combustion of aluminum borohydride in Mach 2 wind tunnel 


capsule continued to burn for about 30 millisec more. The 
shock wave initially off the end of the capsule started moving 
upstream when combustion started, moving rapidly at first 
and then more slowly. Choking of the tunnel was virtually 
complete in about 8 millisec. Restoration of normal flow in 
the tunnel was quite slow. The normal shock did not return 
to the injector position until about 20 millisec after comple- 
tion of the vigorous burning, and normal flow was not estab- 
lished in this section of the tunnel for an additional 20 
millisec, at about the same time that the last driblets of fuel 
from the capsule had stopped burning. The violence and 
highly transient nature of these runs was due to the high 
momentary fuel concentration and precluded any sort of 
precise observation or measurement. 
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With the metal injector shown in Fig. 3, lower (of the order 
of 3 gm per sec) fuel flow rates were used than for the 
glass capsules. Successful steady-state burning of alu- 
minum borohydride was achieved in the small tunnel at 
Mach numbers of 1.5, 2 and 3. Combustion at Mach 
4 was not studied with this injector. A typical time 
exposure of the combustion is shown in Fig. 5. The fuel ig- 
nited at the spark plug, and the flame flashed back to the 
point of injection. Occasionally, the fuel ignited sponta- 
neously at the point of injection. Fig. 6 shows single frames 
taken from high speed direct and Schlieren motion pictures of 
a typical aluminum borohydride flame. From the motion 
pictures, the flashback velocities with respect to the tunnel 
can be estimated. With aluminum borohydride, they were 
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usually in the range from several hundred to 3200 fps. 
If the stream through which the flame was propagating was 
flowing against the flame, the spatial flame velocities in the 
gas were, of course, much higher. Flame velocities of this 
magnitude are usually associated with detonation waves. 
The combustion region is, indeed, associated with a shock 
wave which can be seen in the Schlieren photograph of Fig. 6. 
The static pressure rises associated with combustion along 
the top wall of the test section are presented in Fig. 7. The 
experimental points are average values taken from 16 runs. 
The curve goes through a maximum just behind the point 
of injection, falls off and goes through a minimum, and then 
rises again to a steady value. The forward pressure rise 
(shaded in Fig. 7) is unambiguously the result of combustion. 
However, the downstream rise undoubtedly is due both to Mp 
some combustion in this region and to wind tunnel effects. 
FUEL INJECTION 


These effects arise from the interaction of a shock wave with 
the subsonic portion of the boundary layer. The high pres- 
sure behind the shock wave compresses the boundary layer 
and is fed upstream through the subsonic portion of the 
boundary layer. The result is that the pressures exerted on 
the wall are higher than they would be in the absence of the 
shock wave interaction. In addition, the pressure gradient, 
which is now adverse, could be expected to induce the forma- 
tion of a recirculation zone. 

The shock wave does, indeed, interact with the hot zone. 
This is shown in Fig. 8 which is a composite artist’s sketch 
and photograph of the test section. Because of these wind 
tunnel effects the pressure profile on the top wall (except for 
a small region near the fuel orifice) is more indicative of the 
internal pressures to te expected during combustion in a 
supersonic stream in a duct than it is of the pressure profile 


in Fig. 6 Direct and Schlieren photographs of aluminum boro- 
S — hydride combustion using single fuel injection 
Aluminum borohydride-JP-4 fuel mixtures 6- ee 
2 FORWARD 
Mixtures of aluminum borohydride containing 22, 41 and FUEL 
59 per cent JP-4 fuel by weight were also studied. The 22- a INJECTOR | SPARKPLUG 
and 41-per cent JP-4 fuel mixtures were easily ignited, and A | 
burned well. With 59 per cent JP-4 fuel, two attempts to a - 
achieve ignition were not successful. Both attempts re- z . | 
sulted in weak burning downstream of the spark plug. ¥ PE ae 
= + 
| 
Tandem injections of JP-4 and aluminum borohydride & ; 
Attempts at piloting a JP-4 flame with aluminum borohy- = IF 
dride injected either upstream or downstream of the JP-4 vn . ma 
were successful. The motion pictures and the behavior of —olp 
the pressure instrumentation as well as visual observations DISTANCE FROM FUEL ORIFICE, IN. 
showed that the JP-4 burned only within the aluminum Fig. 7 Static pressure increase on top tunnel wall due to 
borohydride flame zone and only as long as the borohydride combustion of aluminum borohydride. Fuel orifice diameter, — 
flame was present. The JP-4 fuel was unable to sustain 0.0156 in.; pressure, 50 psig 


combustion by itself. When the JP-4 was injected upstream 
of the borohydride fuel, flame would not propagate upstream 
to the JP-4 fuel injection point. 


Pentaborane 


After some difficulty, pentaborane was also successfully 
burned at Mach 2. The lowest volumetric fuel flow rate with 
which pentaborane could be made to burn was about 1.4 
times that for aluminum borohydride. The combustion was 
somewhat different (Fig. 9) from that for aluminum borohy- 
dride. The flame was much larger and penetrated more deeply 
down into the test section. The increased fuel flow coupled 
with the high heat of combustion of pentaborane gave a heat 
release rate for this fuel that was about twice that achieved 
with aluminum borohydride. © 


Fig. 8 Composite sketch-photograph of entire test section of 
The following fuels were not suitable for burning at these wind tunnel showing combustion of aluminum borohydride and 
Mach 2 tunnel conditions: Trimethylaluminum, diethyl- water streams injected in heated region of flow 
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Fig. 9 Direct and Schlieren photographs of pentaborane 
combustion 


aluminum hydride, trimethylborane, triethylborane, alkyl- 
boron hydrides and vinylsilane. Although the fuel flow rates 
were varied considerably, none of these fuels was capable of 
maintaining steady stable combustion. Occasionally unsteady 
combustion occurred behind the spark plug with a few of the 
alkylboron and alkylaluminum compounds. This consisted 
of intermittent flashes of flame associated with the spark- 
ing of the plug. These flashes were all in the vicinity of the 
spark plug or in the tunnel diffuser. It should be noted that 
in flowing from the test section to the diffuser, the fuel-air 
mixture would go through strong shock waves in the diffuser; 
the higher static temperature and pressure in this region is 
more conducive to ignition and combustion of the fuel. 


Alkylaluminum compounds with water injection 


Trimethylaluminum and diethylaluminum hydride ex- 
hibited some evidence of combustion at these Mach 2 tunnel 
conditions when water was injected upstream or downstream 
of the fuel injection point. Light was emitted, and a slight 
rise of static pressure was observed. Alkylaluminum com- 
pounds react explosively with water. Hence the observed 
chemical reaction was probably the hydrolysis of the fuel and 
possible oxidation of the by-products. The reaction flame 
occurred at the point of intermixing of the water and fuel 
sprays. The reaction was mild in comparison to that of 
aluminum borohydride burning in air. 


Flow in the heated region 


Water streams injected into the heated region permitted 
study of the gross nature of the flow field during combustion. 
Fig. 8, constructed from Schlieren pictures taken at various 
stations along the entire length of the wind tunnel shows these 
streams. The sketches emphasize the observations taken 
from the Schlieren motion pictures. The streams are more 
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or less vertical lines originating at the top wall of the tunnel. 
When water was injected through the small holes prior to 
combustion of the fuel, the streams of water adhered very 
closely to the tunnel wall and were atomized very near the 
orifices. When combustion began, the water streams flowe:| 
vertically down into the heated region; they extended to the 
boundary between this region and the main airstream before 
they were deflected into the airstream and atomized. This 
distance is approximately 4 in. from the top tunnel wall in 
the downstream nonluminous portion of the heated region. 
In the upstream luminous region the penetration of the water 
streams was not very deep. The streams near the point of 
fuel injection followed a parabolic path which suggested mere]. 
a thickening of the boundary layer. Examination of the 
motion pictures from which Fig. 8 was constructed revea!s 
that the downstream water streams moved back and fort! 
from the true vertical position, but their preferred orienti- 
tion was in an upstream direction. This indicates that tho 
flow in this region was continually reversing direction, an:| 
suggests the existence of large recirculation zones in the regio: 
of shock interaction with the heated flow. The absenc» 
of shock or Mach waves off the water streams, the dee) 
penetration of the water columns and their upstream orienta- 
tion also suggest that much of the flow in the heated zon« 
downstream of the flame front is subsonic. 


Application of Aluminum Borohydride 
Combustion to Aerodynamic Studies 


Of all the fuels tested, aluminum borohydride was the 
most dependable and easiest to use in existing low temper::- 
ture facilities. With it, the effect of combustion on the aero- 
dynamic properties of the airstreams around several models 
was studied. In al 1-ft supersonic wind tunnel, flat plates 
having chord (distance from leading to trailing edge) lengths 
of 13 and 25 in. were studied at a Mach number of 2.5. In 
the same tunnel, blunt based supersonic wings having 13- 
and 22-in. chords were studied at Mach numbers of 2.5 and 
3.0. In addition several flat plates were studied in a large 
10 10-ft supersonic wind tunnel at Mach 2.4. In all of these 
investigations, the stagnation temperatures were of the order 
of 100 F, and the stagnation pressures were realistic only for 
relatively low Mach numbers at moderately high altitudes. 
In high speed flight, temperatures and pressures would be 
higher, and other fuels could be easily burned. 


Typical Apparatus and Procedure 


An example of the kind of research installation used in the 
aerodynamic studies is given in Fig. 10. This figure shows the 
installation of the 25-in. chord flat plate in the 1 X 1-ft wind 
tunnel. This particular model was mounted at zero angle of 


attack, with the reference surface, adjacent to which combus- 


tion took place, near the centerline of the tunnel. 


Fuel injection system 


The fuel injector used in the small tunnel studies was 
coupled to a }-in. copper tube connected to a fuel manifold 
inside the model. From the fuel manifold, three ;¢-in. 
tubes of approximately equal flow resistance were connected 
to three fuel orifices drilled in the reference surface. The 
fuel orifices were 0.016 in. in diameter and were located as 


shown in Fig. 10. 


The injector was pressurized with dry helium through the 
upper valve to about 50 psig. Then the fuel lines below the 
lower valve of the injector were flushed by a continuous 
stream of low pressure helium (20 psig) which flowed out th« 
three fuel orifices. The helium was turned off by a solenoid 
operated valve just before the lower valve of the fuel injector 
was opened by the drive motor. The fuel was injected 
through the three orifices for a 1- to 5-sec period until it was al! 
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expelled from the injector. Following fuel injection, the 
injector and the lines were flushed with helium to remove 
the residual fuel. 


gnition system 


The ignitor was a }-in. rod which extended from the tunnel 
wall to within spark gap distance (of the order of 35 in.) 
o! the plate surface. The tip of the rod was placed so that it 
was in line with the fuel stream from the lower fuel orifice and 
was about 2 to 3 in. upstream of the blunt trailing edge of 
the model. The ignitor was connected to a 1-joule, repeat- 
ing, capacitance-type power supply which provided a spark 
from the rod to the model surface five times per sec. With 
some of the models where the ignitor could interfere with 
combustion or flow, it was retracted immediately after 
ignition. 


Results and Discussion 


The detailed results of such aerodynamic studies are re- 
ported in a series of NACA research memoranda (2,10,11,12). 
Some pertinent observations made during the course of these 
studies may be of interest. > ald 

Ignition and propagation 

Either flush mounted spark plugs, which caused no dis- 
turbance to the stream, or spark electrodes, which sparked 
through the model to ground, could be used with equal 
suecess. The mixture usually ignited at the spark, quickly 
flashed back to the point of fuel injection, and burned steadily. 
However, two interesting points should be made concerning 
the location of the spark. 

The flame would not propagate upstream unless the mix- 
ture through which it was propagating was part of a bound- 
ary layer in contact with a solid surface. For example, 
in the early flat plate studies, the spark plug was placed 16 
in. downstream of the trailing edge of the flat plate in the 
stream tube into which fuel was being added. Ignition was 
achieved, but the flame would not propagate back to the 
point of injection on the model. Instead, it burned down- 
stream, being held by the spark plug itself or by instrument 
rakes in the downstream portion of the tunnel. 

The flame sometimes would not propagate back to the point 
of fuel injection over curved surfaces. Thus, with the wing 
models, positioning of the spark electrode near the trailing 
edge of the model sometimes resulted in flames which stabi- 
lized themselves 2 or 3 in. upstream of the trailing edge or 
occasionally in the wake just downstream of the base. This 
problem was easily alleviated by simply placing a retractable 
spark plug at the midpoint of the chord. With the spark 
thus located, achieving flashback to the point of fuel injec- 


tion presented no problems. 1" 


Examples of the kinds of pressure profiles which are ob- 
tained in such studies are shown in Fig. 11. In this figure is 
plotted the change in pressure coefficient resulting from com- 
bustion, as a function of the distance downstream from the 
point of fuel injection, for the various models tested in the 
| X 1-fttunnel. These are compared with a similar plot of the 
data for a 25-in. chord flat plate model in a very large 10 X 
10-ft wind tunnel. The change in pressure coefficient is the 
local static pressure change due to combustion divided by the 
(lynamic pressure of the undisturbed stream. The ordinate is 
therefore a measure of the pressure rise on the airfoil caused 
yy combustion in the adjacent stream. 

The profiles are qualitatively quite similar in their upstream 
regions. There is a rather sharp pressure rise in the vicinity 
of the fuel injection region, which reaches a maximum a few 
nches downstream of the point of injection; the pressure 
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Examples of pressure profiles obtained 
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Fig. 10 Flat plate model installation in 1 x 1-ft test section 


MODEL CHORD AVERAGE ANGLE OF 
anal LENGTH, We, ATTACK, 
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Fig. 11 Comparison of centerline chordwise profiles of change in 


local pressure coefficient along lower surface due to combustion 
of aluminum borohydride in airstream below various models 
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This 
The 
13-in. wing, at Mach 2.96, clearly shows this characteristic 


rise gradually decreases in the downstream direction. 
is the main characteristic of the pressure rise profile. 


profile. The second pressure rise on the longer chord models 
is caused by reflection of the flame and leading edge shock 
waves from the tunnel walls back into the downstream por- 
tions of the flame zone where they interact with the heated 
flow, causing a subsonic recirculation region and further com- 
bustion of unburned fuel with air entering this region. Fig. 
11 shows that, with the 25-in. flat plate in the large 10 < 10-ft 
tunnel, as with the short models in the 1 X 1-ft tunnel, there is 
no large second pressure rise region over the downstream 
portion of the model. In the large tunnel the wall was over 
6 ft from the flat plate surface. Therefore, the point where 
the reflected oblique shock wave impinged on the heated wake 
was so far downstream of the base of the model that there was 
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Fig. 11 suggests that, because of the “fattening” effect of the 
heated flow about the model, the ratio of tunnel to model 
size should be much larger than is usually considered neces- 
sary to avoid tunnel effects. 


A study of the combustion of various highly reactive fuels 
injected through the top wall of a 3.84 X 10-in. wind tunnel 
into a supersonic airstream disclosed that: 

1 Aluminum borohydride, pentaborane and mixtures of 
up to 41 per cent JP-4 fuel blended with aluminum borohy- 
dride could be burned in a Mach 2 airstream under the condi- 
tions of this experiment without the use of a flameholder. 
Combustion of these fuels gave associated pressure rises in 
the reaction zone. 

2 Tandem aluminum borohydride and JP-4 injections 
gave high heat release rates, which tended to choke the super- 
sonic flow in tunnels of this size. 

3 JP-4 fuel could be burned in a Mach 2 airstream under 
the conditions of this experiment only as long as a piloting 
flame of aluminum borohydride was present. 

4 Trimethylaluminum and diethylaluminum hydride 
could not be ignited at Mach 2 tunnel conditions. How- 
ever, when water was simultaneously sprayed into the 
tunnel, there was a luminous reaction which produced heat. 

5 Trimethylborane, triethylborane, alkylboron hydrides 
and vinylsilane could not be ignited or burned in the tunnel 
test section. These fuels frequently did ignite in the tunnel 
diffuser; therefore, they might be combustible under condi- 
tions where the recovery temperatures would be higher than 
those of the present study. 

6 Studies in which the heated region was probed by water 
injections, indicated that the flow downstream of the flame 
front was subsonic and recirculating. 

The combustion techniques developed in the 3.84 X 10-in. 
tunnel were applied to aerodynamic studies of the effect of 
heat addition by combustion of aluminum borohydride on 
supersonic airstreams. The detailed aerodynamic results 
have been reported elsewhere in the literature. The follow- 


obtained data. 


T IS generally assumed that the problems of thermal con- 


trol of satellites are solved. Explorer I had inside a com- 
fortable temperature of 20 C, maintained over an operating 
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efficiency of the chemical batteries). 


ing observations should prove helpful to those wishing to do 
further aerodynamic studies of this nature: 

1 The presence of a boundary layer is extremely helpful, 
if not essential, in achieving flashback and stable combustion 
even with the most highly reactive fuel of this investigation, 
aluminum borohydride. 

2 The technique is well suited for aerodynamic studies 0: 
external combustion. However, the test facility must b 
large enough, or the model short enough so that disturbanc« 
which reflect from the tunnel walls intersect with the ho‘ 
gases well downstream of the region of aerodynamic inter 


est. 
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Thermal control of the Explorer satellites is discussed. The theoretical studies that were made 
prior to the launching of these satellites are described, and examples of relationships are presented 
in graphs. The lower limit of the instrument temperature was zero deg C (determined by the 
The upper limit was specified as 65 C (based on long-time 
temperature limits of transistors in the electronic package). This paper will relate some of the 
studies, describe the measuring results of telemetered temperatures and evaluate expected and 


period of 3 months. The temperature variations around this 
mean value were not more than +20 C; i.e., the climate on 
board the satellite had fewer variations than the tempera- 
ture, for instance, in Huntsville, Ala. However, the ther- 
mal problems of satellites require a number of studies and 
prototype testing prior to launching. 

Control of Explorer I was maintained between 0 and 40 C 
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during 90 days of operation; Explorer III was maintained be- 
tween 0 and 40 C up to 43 days after launching. Then the 
transmitter ceased to operate; this was attributed to freezing 
of the batteries. 

Explorer IV was controlled between 19 and 55 C during 
the operational time of 60 days. 


Parameter Study of Thermal Problems ze 


The temperature of a satellite is dependent on: 

1 Orbital characteristics. 

2 Surface properties, especially emissivities, of skin and of | 
internal surfaces. 

3 The attitude of the momentum vector of a spinning — 
stellite to the direction of the sun. 

4 The attitude of the spin axis to the direction of the sun. 

5 The attitude of the satellite and attitude of the momen- 
tum vector to the center of Earth. ‘ 

6 Heat capacity of skin and of instruments. 


and instruments, and between parts of the skin. Convective 


7 Conductive and radiative heat transfer between skin F 


leat transfer was excluded for the Explorers by letting the air — 


escape through gaps. 

8 Internal heat release of instruments. 

9 Aerodynamic heating for satellites with a perigee of 
200 miles or lower. 

10 Environmental conditions, such as solar constant, — 
albedo, Earth radiation. 


The Explorer satellites are spin stabilized; therefore, 


temperatures on individual points of the surface equalize in 
any plane perpendicular to the spin axis. 

Fig. 1 shows a cutaway of Explorer I. Explorer I data 
and instrumentation have been published elsewhere, e.g., in a 
previous paper by the author (1).2. Therefore, only points 
that affect temperatures are described here. The satellite is 
shaped like a pencil. The empty fourth stage motor stays 
attached, but it is thermally and electrically insulated from 
the cylinder. The wire antennas were used only in Explorer 

The parameter study was performed prior to Explorer I, 
and, to some extent, prior to each Explorer firing. Each of 
the parameters was treated as an independent variable and 
varied in a range wide enough to cover variations of the 
missile characteristics and their effect on the injection data, 
firing date and time, satellite properties and environmental 
conditions. The thermodynamic and celestial mechanics 
equations were programmed on an IBM 704 computer. Each 
parameter was varied independently by choosing three, four 
or five different values. 

An important factor that influences the satellite tempera- 
tures is the time in sunlight. For the Explorer-type satellite, 
this value may fluctuate between 63 and 100 per cent. Ac- 
cording to Krause (2), the shadow angle is given by 


sin X = sin B/cos 6 {1] 
where 
X = angle between 0600-1800 hours line and crossover from 
shadow to sunlight or from sunlight to shadow : 
8 = angle at center of Earth between satellite and tangent to : 
Earth 
6 = angular distance of sun from orbital plane 


The following approach was taken to solve the temperature _ 


problem with so many parameters, and to arrive at design 
criteria. 

An analysis was made of the crossover points of the satellite 
in and out of the sunlight for each revolution as a function 
of the orbital characteristics and the position of the sun. 
The angle 6 is a function of the orbital parameter and the 


2 Numbers in parentheses indicate References at end of paper. 
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sun’s position. Parameters used to describe the orbit were 


= inelination of orbit to Equator 

2 = right ascension of ascending node measured from spring 
* point 

= argument of perigee 

e = eccentricity of orbit ; 

R, = radius vector at perigee 

@ = time angle measured in orbital plane from east direction 

to perigee 


The position of the satellite within this orbit is given by the 
true anomaly v. 

The position of the sun enters into the computation by its 
longitude LZ) measured along the ecliptic. A mean motion 
of the sun in the ecliptic plane was assumed. A check was 
made for some cases, using the true position from the Ephem- 
eris and Nautical Almanac, and the difference was found 
to be negligible. 

The quantities entering into Equation [1] can be expressed 
by the parameters above by coordinate transformation be- 
tween orbital plane and ecliptic plane to the equatorial 
coordinate system that is commonly used for orbit param- 
eters. 

For orbiting during an extended period of time, it must 
be considered that the orbital characteristics change. First- 
order perturbation theory as developed by Krause (3) and 
by Cunningham (4) has been applied to compute these 
changes. Only secular changes have been considered, be- 
cause the small periodic changes have only negligible effect 
on the temperature problem. For this purpose, the time 
dependent angles 2, @ and a have been derived as functions 
of the hour of launching Hy and the day of launching Dp. 
The changes with respect to time are then expressed by the 
linearized equations 


Q = & + ki (D — Dp) [2] 
= go + kx (D — Do) [3] 
= a + k3(D — Dp) 
where 
ki, ke, kj = constants 


D — Dy = number of days after launching 


For the conditions of Explorer I, the coefficients are 


k, = —4.35 deg per day (rotation of orbital plane) 
ke = 6.49 deg per day (rotation of line of apsides) 
ks = 0.9856 deg per day (mean motion of sun) 
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Fig. 2. Time in sunlight for a satellite as a function of day and 
hour of firing. Maximum values during first 60 days of orbiting 


From the crossover points of the satellite in and out of the 
sunlight, the time can be computed by introduction of the 
eccentric anomaly F and the use of Kepler’s equation 


T, — 
i= #2 


ratio of time in sunlight to period 
period of satellite a) 


where 


time when satellite enters sunlight ==” 
time when satellite enters shadow —— 
eccentricity 


Attitude angle of satellite axis with respect to the sun is 
determined by the original coordinates with respect to the 
horizon coordinate system. In this system, the altitude 
angle is zero and the angular direction is given by the Earth- 
fixed azimuth. The initial attitude can be computed from 
these coordinates and the position angles of the sun as a func- 
tion of the day Do and hour Hp of launching. 

Changes of the attitude angle to the sun with time are 
dependent on the movements of the sun with respect to the 
original launching plane and not on other time variable 
changes of this orbit. It is assumed that the initial momen- 
tum vector of the spinning satellite is maintained in space. 
Exchange of angular momentum from the longitudinal axis 
to an axis of greatest moment of inertia has been considered. 
However, this initial concept has been refined on the basis of 
results obtained from the Explorer measurements. 

The surface temperatures of the satellite depend on the 
emissivities of the skin with respect to short wave solar 
radiation and infrared radiation from Earth and from satellite 
skin, and on the coefficients of solar radiation (solar con- 
stant), albedo and Earth radiation. 

The following two simultaneous differential equations have 
to be solved (1,5,6) 


A,aSD, A,aBSD, A;eES AyoeT,' 
—-CT.+Q=0 


GT, Cx T:) + C3a( =) [7] 
= projected area with respect to sun direction 
absorptivity for solar radiation 


solar constant 


step function (D, = 1 for sunlight, D,; = 0 for shadow) 
effective area of satellite with respect to albedo 
= step function (D, = 1 for orbit in hemisphere +90 from 


= — [(F, — esin E,) — (FE, — esin [5] 


Pe sun direction, D, = 0 for orbit in hemisphere opposite 
sun direction ) 


= albedo of Earth 

A; = effective area of satellite with respect to Earth radiation 

€ = emissivity with respect to infrared radiation 

E = ratio of infrared heat flux from Earth to solar constant 

A, = total surface area of satellite 
| 

= Boltzmann constant 

7, = skin temperature 

C; = heat capacity of skin : 


7. = rate of temperature change of skin | 
C, = heat capacity of instruments = 
T; = instrument temperature 

7; = rate of temperature change of instruments 
Q = heat release of instruments 

C; = heat transfer coefficient by conduction 

Cs = heat transfer coefficient by radiation 


The areas are obtained from 


+ H? ) 
A, = A;° 1— 20s 
QRH + H? ) 
io} 
where 
A,° = maximum area of cylinder as seen from side *.. 

Y = attitude angle to sun 

f(©) = function of attitude angle ¢ of satellite axis to radius 


vector to center of Earth. (It was approximated by 
a five-power polynomial of cos ¢) 


Ry = radius of Earth at Equator 
= altitude of satellite 
6 = angle of radius vector and direction to sun 


The thermal design is based on thermal insulation between 
skin and instrument package. The skin temperature can 
fluctuate in a wide range during each sunlight and shadow 
cycle. The instrument temperature is affected very little 
by this cyele. It is mainly dependent on long-time changes 
of time in sunlight, attitude to sun and Earth, long-time 
changes of surface emissivity and changes of the environ- 
mental conditions. 

In order to eliminate, as far as possible, changes of surface 
emissivities by impact of micrometeorites and other environ- 
mental conditions whose nature and influence is still to be 
determined, surfaces were selected that are insensitive to 
such effects. The selection for Explorers I and III was 
sandblasted stainless steel, and 25 per cent of the area was 
covered by stripes of aluminum oxide known under the trade 
name of Rokide A. Explorer IV was sandblasted stainless 
steel only. 

The discussion of the parameter study up to now has been 
made under the assumption that all parameters are variables. 
However, if the thermal design has been determined on the 
basis of the analysis, all on-board parameters become fixed. 
In further analysis, it was assumed that the thermal proper- 
ties of the satellite do not change during its useful operating 
time of 2 to 3 months. At this time, usually, an additional 
variable becomes fixed: The date of the firing. 

The following parameters remain variables: Environ- 
mental conditions, orbital characteristics at the injection 
point of the satellite into orbit, hour of satellite launching. 
The environmental conditions are subject to changes that 


E are not completely predictable and will be investigated in 


future satellites. For the solar “constant,” a variation of 


_ +3.5 per cent was assumed. Albedo was varied between 


24 and 54 per cent and the ratio of Earth radiation to solar 
constant between 15.9 and 19.3 per cent. Change of orbital 
characteristics at the injection point is due to dispersion of 
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Fig. 3 Per cent time in sunlight (7’,) of Explorer IV vs. hour of 
firing and days after firing 


booster and upper stages. A 2o variation of all these param- 
eters wasassumed. The main influence was due to dispersions 
of the directions and the magnitude of the velocity vector. 
Variations of the injection values and of the environmental 
value determine upper and lower limits of temperatures. 


The hour of the satellite launching can influence the initial © 
plane of the satellite by 360 deg and the attitude of the satellite — 


axis in space in a wide range. 

The selection of the firing time for all Explorers has been 
used as a means to maintain temperature control. Ex- 
plorers I and IV had a time in sunlight between 65 and 80 per 
‘ent during their operating time. 7 


Results of the Parameter Studies. 

I will not describe all details of the parameter studies, 
especially since some of the results have already been pub- 
lished (1,7). Therefore, it will suffice to show the relation- 
ships in a few graphs and a table. 

The time in sunlight was computed for orbital characteris- 
tics of Explorer IV as a function of day of firing Do and hour 
of firing Ho. The 704 computer program was run for a simu- 
lated orbiting time of 60 days for each combination of Do and 
H,. The maximum values obtained at any point during this 
time are plotted in Fig. 2 as a function of Do and Ho. The 
percentage time in sunlight is used as the parameter-like alti- 
tude lines on maps. For the selection of a firing time, vari- 
ations of the initial orbital conditions had to be analyzed. 
The following effects were found. The areas of high time in 
sunlight values increase with the increase of: The eccen- 
tricity of the orbit; the injection altitude; inclination angle; 
and error angle in pitch, both plus and minus, and in yaw to 
the left for firing from Patrick Air Force Base. 

A better understanding of the maxima of time in sunlight 
T, can be obtained from Fig. 3 which shows 7, for Explorer IV. 
The day of firing is no longer a variable but fixed on July 26, 
1958. The maxima of 7; are plotted as a function of the hour 
of firing and the days after firing. The maxima shift with 
the hour of firing, with a constant slope which is given by the 
daily change of the ascending node and the daily change due 
to the relative motion of the sun. The magnitude and ex- 
tension of the maxima are determined by a combination of 
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Fig. 4 Time in sunlight (7',) vs. hour of the day of launching 
for Explorer I. Maximum values during first 60 days of orbiting 
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these two factors and the progression of the line of apsides. 
This effect does not exist for circular orbits; for elliptical 
orbits, it increases strongly with the eccentricity. 

The effect of eccentricity is shown in Fig. 4, which is a plot 
of the maximum time in sunlight vs. the hour of firing for 
Explorer I. In this case, the program was also run for a 
simulated orbiting time of 60 days. Time in sunlight 7’, vs. 
the days D — Dy after launching is shown in Figs. 5 and 6 
for Explorer V. Fig. 5 shows the variation of 7, for the 
hours Hy 0-10 as parameters, and Fig. 6 shows the same for 
the hours 12 to 22. aes 
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Fig. 7 Attitude angle (+) of Explorer IV payload to sun 
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Fig. 9 Temperature of Explorer I. Nose cone instrument 
temperature 


Table 1 Linearized influence parameters for Explorers 
I and III 
Change of 
Nominal instrument 
Parameter value Variation temp., C 
albedo 0.34 +0.20 +5 
ratioofEarth 0.175 +0.0175 +8 
radiation to 
solar constant 
solar constant 2cal/min-em? +0.07 +23 
timein sunlight 0.73 +0.1 +10 
projected area 90 percent of +10 per +8 
to sun max. area cent 
altitude of 1300 km  +£100km 
satellite 
ratio of emis- 1.4 5 
sivities 


The effect of variations of several parameters applicable 
to Explorers I and III is shown in Table 1 in the form of 
linearized influence parameters. The expected temperature 
_ changes are quite large. A more detailed analysis and evalu- 
ation of Explorer results showed that variations of albedo, 
Earth radiation and satellite altitude affect only the skin 
temperature fluctuation, but had little effect on the mean 
temperature. A better analysis of this problem will be pos- 
sible after results of the radiation balance are available from 
future satellites. Variations of solar constant are mainly 
due to the variable distance of Earth from the sun. For an 
operating time of 60 days, the variations are negligible. 

This leaves as the main influencing parameters the time in 
sunlight and the projected area to the sun. Variations of 
the ratio of emissivities can have quite a decisive effect on the 
temperatures, and overshadow other effects. However, with 
the precautions taken to select surfaces insensitive to changes, 
it is reasonable to assume that the ratio of emissivities stayed 


constant. 
* > 


Attitude of Satellite 


The initial attitude of the satellite to the direction of the 
sun is easily determined from the direction of the velocity 
vector at injection into orbit and the known direction of the 
sun. This holds true for a spin stabilized satellite whose 
spin axis is aligned horizontally to Earth at the injection 
point and whose angular dispersion is negligible. These 
conditions were fairly well satisfied for Explorers I, III and 
IV at the time of injection. Fig. 7 shows the initial attitude 
angle to the sun for Explorer IV as a function of the day Do 
and the hour of firing Ho. This initial attitude angle will 
change due to the relative motion of the sun, as shown in 
Fig. 8, as a function of the number of days after the firing. 
For the pencil-shaped Explorer satellites, the spin axis is not 
stable, since it has the smallest moment of inertia. The 
initial momentum vector was considered unchanged in space 
(pointing toward the same fixed star). However, an ex- 
change of momentum from the minor axis of inertia to a 
major axis of inertia had to be considered in the analysis. It 
can be shown that such a momentum transfer is possible, if 
the energy equation can be satisfied. Since the spin around 
a major axis of inertia means a reduction of the spin rate, this 
mode of rotation corresponds to a lower energy level. There- 
fore, the transfer of momentum takes place if rotational 
energy is dissipated by internal friction. It was assumed 
prior to Explorer I that this transfer takes place during the 
active lifetime of the satellite instrumentation. Fig. 9 shows 
theoretical satellite temperatures of Explorer I for the two 
extreme cases: 1 Spin axis unchanged around longitudinal 
axis of the satellite; and 2 spin around the initial axis with 
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the satellite aligned perpendicular to this axis (flat spin). 
The graph shows that temperatures in case 1 start out with a 
comfortable 295 K and stay above 273 K for 30 days, but 
then drop rapidly below the lower operating limit of 273 K. 
As the temperature measurements obtained by telemetry 
show, this did not happen, but all measuring points shown 
in Fig. 9 follow the curve for case 2 after a few days of orbit- 
ing. From an analysis of the antenna patterns, it was deter- 
mined that the transfer in Explorer I had occurred during 
the first revolutions. This is the first case in which tem- 
perature measurements were used to determine the behavior 
of the satellite in space. Curve 2 of the flat spin is shown 
as an extreme case. The actual attitude and mode of spin 
is between case 1 and case 2. Also, a random tumbling 
causes the temperature to equalize at 295 K. To distinguish 
between case 2 and a random tumbling, a higher accuracy 
was needed. Further analysis showed that due to the 
nternal damping, the flat spin was reached in all cases. In 
‘xplorers III and IV, it took about 10 days. These satellites 
lid not have the flexible antenna wires. The dissipation of 
nergy took much longer. The cause for this dissipation is 
issumed to be due to inelastic flexing and damped vibrations 
nside the satellite. External forces on the satellite to be 
-onsidered are magnetic forces and forces due to interaction 
with the atmosphere. Magnetic damping was found to be 
of minor influence. However, interaction with the atmos- 
phere during each perigee pass means a transfer of momentum 
which causes a change of the momentum vector in space. 
The long-time vector sum of changes can have a decisive 
effect on the thermal design and was, therefore, considered 


for the basic design philosophy of Explorer IV. 


Thermal Coefficients and Thermal Tests 
The thermal design of the Explorers was based on the 

following concepts: The control is passive; i.e., no variable 

surfaces or heat switches are applied; and the instrument 
package is insulated from the skin of the satellite to minimize 
short-time fluctuation during one revolution. 

The coefficients in Equations [6 and 7] have to be known. 
The Explorers consisted of two parts, the cylinder and the 
cone, separated by insulation. The surface properties were 
adjusted so that the average temperatures in both parts 
were about the same. However, if the insulation is too good, 
the internal heat release will cause the package to overheat. 
Without insulation, the instruments would take part in the 
strong temperature variations of the skin. 

A prototype testing program was carried out to determine 
the soundness of the thermal concept and to measure the 
coefficients entering into Equations [6 and 7]. 

Tests were performed by applying as sudden a temperature 
rise as possible (step function), and to measure the rise of the 
inside temperatures at various places, such as the transmitter, 
batteries, structure. Fig. 10 shows a typical test run for 
Explorer IV. 

The temperature range a satellite of the Explorer type might 
assume without thermal design is between 200 and 380 K, 
as shown in a previous paper (1). This wide range applies 
for a constant value of the ratio of emissivities a/e of 1.5. 
For variable a/e the range would increase. 


Results of Temperature Measurements 


Fig. 9 shows average daily temperatures obtained by te- 
lemetry from Explorer I. The relation between the mode of 
the spin and the measured temperatures has already been 
discussed. Fig. 11 shows the temperature changes as ex- 
pected during one revolution of Explorer I inside the instru- 
ment package and on the satellite skin. The dotted line 
shows the mean of measured values during the first 10 days 
of orbiting. The dip in the skin temperature during the 
sunlit portion occurs when the satellite reaches its apogee. 
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Fig. 12 Explorer III, per cent time in sunlight vs. days after 
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Explorer III is an example of a satellite that was allowed 
have 100 per cent sunlight during 4 days of its operating 
ne. Fig. 12 shows the time in sunlight vs. days after firing. 
1e Maximum temperature measured was 40 C 36 days after 
Fig. 13 shows the theoretical temperature com- 


puted with the following assumptions for the satellite attitude: 


0-10 days, spin axis unchanged. 

After 10 days, flat spin. 

Upper curve, momentum vector at zero angle to sun vector. 
Lower curve, momentum vector at 90 deg to sun vector. 
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| ‘Fig. 14 shows the average of daily measurements and a 
third-order least square fit. Reception of data ceased after 
43 days of operation, although the batteries were set up for 
90 days. Signals were received twice at irregular intervals on 
the 62nd and 71st days. A number of theories have been 
offered to explain this. I would like to add one that seems 
plausible to me, although I cannot prove it by measurements. 
After 43 days, the temperatures dropped below zero deg C. 
This caused the batteries to freeze, and no more signals were 
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received. During short intervals later on, the temperature 
rose above freezing, reviving the batteries and causing the 
transmitter to operate again. These later signals were too 
weak to enable telemetry evaluation. The temperatures 
dropping below zero prevented the reception of telemetry 
data to prove it. 

Fig. 15 shows skin temperature and instrument tempera- 
ture for a typical case of one daylight and night cycle of Ex- 
plorer IV. Time in sunlight is 75 per cent. The apogee 


mined from variations of the subcarrier oscillator frequencies 
The temperature dependency of the oscillators was calibrate: 
by Ludwig of the State University of Iowa prior to the launch 
ing. Fig. 16 shows the results of this temperature determina 
tion during the operating lifetime of 55 days. Two of th 
_ switching frequencies of channel 3 of the Explorer IV payloa: 
_were used. Results of both are consistent, one frequency lev: 
always yielding higher values. Data were taken at th 
ABMA experimental tracking station at Huntsville, Ala 
_ Measurements from this station were very good, and ha 
almost no scatter. Measurements from other IGY station 
- confirmed the general trend of the measurements. As : 
most likely value, the arithmetic mean for each day has beer 
used for the evaluation. The maximum temperature meas- 
ured was 55 C, and the minimum temperature was 19 C 
The overall average was near the design value of 39 C. A: 
Fig. 17 shows, parameters of Explorer IV were chosen so 
that the time in sunlight was subject to small changes onl) 
during 60 days. The strong variations of temperatures 
were hard to reconcile with this smooth sunlight curve 
The projected area as exposed to the sun is shown in 
Fig. 18 for the two extreme cases: 1 Spin about the origina! 
longitudinal axis; and 2 flat spin maintaining the original 
momentum vector. A transition from case 1 to case 2 can 
be expected within the first 5 to 10 days. The drop of the 
measured temperature from an initial value of 40 to 30 C 
seems to be connected with this changeover, as was observed 
in Explorer I. Fig. 19 shows a superposition of measured 
temperatures and theoretical temperatures, computed for a 
series of exposed area ratios. The minimum area for a flat 
spin, with the sun in a plane perpendicular to the momentum 
vector, is an average of 2/7 or 63 per cent. Fig. 19 shows 
that the measured temperatures fit very well into this grid 
of theoretical temperatures. The point where an area of 


100 per cent is exceeded corresponds to a temperature differ- 
ence of 4 C. This is well within the accuracy, considering 
variations due to measuring error, variations of emissivity, 
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Fig. 16 Temperature of Explorer IV. Based on channel 3 
frequencies received by ABMA tracking station 
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occurs during the shadow period. The instrument tempera- 
: ture variations are +5 deg, whereas the skin temperature 
varies between 255 and 330 K. Temperatures were deter- 
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Fig. 17 Per cent time in sunlight (7'r) for Explorer IV vs. days 
after firing 


variation of environmental conditions, accuracy of orbital 
characteristics and accuracy in determining time in sunlight. 
It can be concluded that the constancy and accuracy of these 
factors has been fairly close for Explorer IV. It is not in- 
tended to generalize this for other satellites, because of the 
niany factors involved. 

The temperature measurements allow an analysis of the 
spinning behavior and the direction of the momentum vector. 
The following interpretation is offered. 

At the time of launching, the angle between the longitudinal 
axis of the satellite and the sun’s direction was 75 deg. This 
was simultaneously the direction of the momentum vector. 
The area exposed to the sun was 97 per cent. 

In the first 6 days of orbiting, the mode of spin changed to 
a rotation around a major axis of inertia. The momentum 
vector remained essentially the same. The angle between 
the momentum vector and the sun vector would decrease to 
7l deg. It is interesting to note that the area did not reduce 
to the 63 per cent minimum, but to about 68 per cent, which 
corresponds to an angle of 77 deg between momentum vector 
and sun vector. 

After 6 days, the flat spin had fully developed. Inter- 
action with the atmosphere became strong and caused a rota- 
tion of the momentum vector. The rate of change was up 
to 15 deg per day. This caused the exposed area to fluctuate 
between the extremes of 100 and 63 per cent. After 26 days, 
the rate of change seemed to decrease considerably. How- 
ever, since this is a three-dimensional problem, this break in 
the area curve may indicate merely a change in the direc- 
tion of the rotation of the momentum vector. 

It is possible to draw conclusions from the temperature 
measurements on the behavior of satellites in space. Since 
three angles are involved and only one measurement made, 
two types of ambiguities remain: Determination of the 
quadrant in which the momentum vector is located, and de- 
termination of whether the plane of momentum vector and 
sun vector has rotated from its original location. To resolve 
such ambiguities, additional measurements are needed. An 
entirely independent method of attitude determination has 
been used and published elsewhere (8). 

This method uses the signal strength of received radio 
signals. Fig. 20 shows the correlation of the attitude ob- 
tained by temperature evaluation and the evaluation of radio 
signals. 


Conclusions 


The temperature of satellites can be controlled by passive 
means. The on-board instrument temperatures can be kept 
between 0 and 65 C during an operating time of 60 to 90 
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Fig. 18 Theoretical projected area (Az) of Explorer IV vs. days 
after firing 
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Fig. 19 Measured temperature of Explorer IV superimposed on 
theoretical curves 
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days. Temperature variations of the skin can be kept be- 
tween —30 and 90 C. 

The thermal control of the Explorer satellites was based 
on an analysis of all influencing parameters. Relationships 
were developed by numerical evaluation on an IBM 704 com- 
puter, taking all variable factors into account. 

Time in sunlight and attitude angle to the sun can be con- 
trolled by the time of firing. The launching date cannot be 
freely chosen, but the hour of firing allows an angular vari- 
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ation of the orbital plane by 360 deg. The time of firing has 
been used with the Explorer satellites as an effective control 
parameter. 

The thermal design of the satellite was based on insulation 
between skin and instruments. Heat transfer coefficients by 
conduction and radiation were determined theoretically and 
experimentally by prototype testing. 

The interval of +32 deg to which a satellite’s instruments 
can be controlled is determined by the external parameters, 
such as orbital characteristics, firing time, solar constant, 
albedo, Earth radiation, attitude, mode of spin, and their 
variations with time. 

The nominal value of the temperature which is the mean 
of the interval is determined mainly by the on-board param- 
eters: Absorptivity to solar radiation, emissivity in infrared, 
heat capacities, heat transfer coefficients and internal heat 
release. For a body without heat release and without heat 
capacity, the ratio of emissivities has the decisive influence. 
See Table 2. 

Explorers I and IV had times in sunlight between 65 and 
80 per cent, during the operating time of 90 and 60 days, re- 
spectively. Explorer III went to 100 per cent sunlight 
during 4—5 days, 30 days after launching. Measured tem- 
perature during this time was 40 C. 

Due to the low design temperature, the upper temperature 
limit was not exceeded. It is very likely that the tempera- 
ture dropped below the lower limit after 43 days, causing the 
transmitter to drop out. 
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A general condition is derived which must be satisfied by an aircraft or space vehicle in order to 
reach a predetermined point in a minimum time. 


Table 2 Explorer temperature control 


Mean design Measured temp., C 


Explorer temp., C Max. Min. 
I 20 40 0 

III 20 41 0° 
IV 40 55 19 


* Actual temperature lower, but no measurements be- 
cause transmitter ceased operating. 
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The thrust, lift and drag can depend, in an arbi- 


trary (differentiable) manner, on the instantaneous position and velocity coordinates (or Mach 
number) as well as on the angle of attack between the directions of the thrust and the velocity. In 
order to make the solution applicable to space vehicles as well as to aircraft, the cases of a spherical 


HE PROBLEM of determining brachistochronic trajec- 
tories has been investigated by numerous authors from 
different points of view and under different assumptions (2-6).? 
The most general treatment probably originates from Miele 
(1), who considers the present problem from the standpoint 
of a variational problem of Mayer type. In that paper, 
however, the minimum conditions (Euler equations) still seem 
Received April 5, 1959. 
1 Chief Mathematician. 
* Numbers in parentheses indicate References at end of paper. 
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and a flat Earth are treated separately. 


to be quite inconvenient since, in addition to the trajectory 
parameters, a number of Lagrangian multipliers have to be 
determined. Only in very special cases can these equations 
be integrated in a closed form. 

It is possible, however, to arrive at a rather simple minimum 
condition if the same problem is treated as a classical vari- 
ational problem. This minimum condition and the equations 
of motion together represent the necessary number of equa- 
tions for the complete determination of the desired trajectory 
and the required flight parameter, © 
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The Case of a Spherical Earth 


In aspace fixed system of coordinates z, y, the instantaneous _ 
position of the space vehicle can be described by the radius — 


) 


a 
y 


r cos g 


Xx 
rsin 


Here r and ¢ are the circular polar coordinates. The in- 


stantaneous velocity is given by 


(3) ( 


where a dot stands for the derivative with respect to the time 
it. The following forces are acting on the center of gravity of 
the vehicle: Thrust 7, lift ZL, drag D and gravitational force 
G (Fig. 1). If mis the instantaneous mass of the vehicle, the 

«quation of motion is 


Vi 


‘) 


V2 


— rgsin g 
isin g + rgcos 


) 


y 


= 


(3) 


faking into account that the drag vector D is pointing in the ~ 
direction opposite to that of the velocity vector X, that L 
s perpendicular to X, that T forms the angle of attack a with 
¥, and that G is pointing toward the center of Earth O, the | 
equation of motion [3] can be written in components as_ 
follows 


d 
T+D+L+e 


+ Y cos ¢) = (r cosa — 6 mv) 
v 
(7 sin + A) [4a] 
m(b2 + sin = (Tsina@ + d) +: 
v 
v2 
(r cosa — 6 — mv) -— [4b] 
v 


Here 7, 6, \ and y stand for the absolute magnitudes of 7, D, 
L and G/m, respectively, and the instantaneous velocity v is 
given by 


v=lyl = = Vet re [5] 


The quantities 7, 6 and ) are, in general, empirical functions 
of “altitude’’ r, Mach number / and angle of attack a. If 
the physical properties (density, temperature, etc.) of Earth’s 
atmosphere are assumed to be spherically symmetrical and, | 
therefore, functions of r only, the Mach number can be ex-— 
pressed in terms of rand v. For the mathematical treatment 
of this problem it suffices to assume these quantities to be 
arbitrary (differentiable) functions of the local coordinates 
r, g, the velocity components v, and ve, and the angle of 
attack a. The fuel consumption m is supposed to be a 
function of the same variables. Thus we assume 


6 = Or, U1, v2, a) 
Ar, Y, Vi, V2, a) 
m = m(r, Y, M1, V2, a) 
Y = golt?/r? 


Here FR is the radius of Earth, and go is the gravitational 
acceleration at Earth’s surface. 

Now, taking into account that the velocity components 1 
and v2 are functions of r, 7, g¢ and ¢ according to Equation 
{2], one finds that the equations of motion [4a and 4b] repre- 
sent two differential equations for the three unknown func- 
tions r(t), g(t) and a(t). If, therefore, one of these three 
functions is chosen to be a known function of time ¢t, the re-— 
maining two functions are already determined by the equa- 
tions of motion. This means that there is only one degree 
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of freedom which can be used for the determination of the 
minimum time trajectory in question. We now choose the 
function a(t) so that the time to reach a final point (71, ¢g) start- 
ing from a point (79, go) is aminimum. For this purpose it is 
necessary to derive an explicit expression for the flight time 


itself. This can be done by first excluding the case where 
ti + y cos gy =O 7a] 
te +ysin g = 0 [7b] 


These equations represent exactly Kepler’s equation of motion 

of a mass under the influence of a central force. Dividing one 

of the equations of motion [4a and 4b] by the other, one ob- 


tains 
i +yeos Av, — Bry 
vo + sin Bu + Avy 
where 
A(r, 11, @) = T Cosa — 6 — mv [9a] 
B(r, 9, 1, 2, @) = Tsina +X [9b] 
Rearrangement of Equation [8] leads to 
(Bo, + At») = bo( Ary [10] 


— Boz) sin g — (Br; + cos ¢] 


Here we introduce ¢ to be the independent variable and de- 
note the derivative with respect to g by a prime. It follows 
that 
01 '(Bu; + Avz) — — Bre) 
— Boz) sin g — (Buy, + Ave) cos ¢] 


dg =dt [11] 


where the time ¢ is separated, if, in the expression on the left- 
hand side, the quantities r, v, v1, v2 and @ are considered to be 


functions of g. Integration with respect to g from ¢o to ¢1 
leads to the relation 


— Bre) 


= } 


sin g — (Bu, + Ar) cos ¢] 


Fide =, — to 
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L 
vector 7 
om, 
Fig. 1 Illustration of the coordinate system . 
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On the right-hand side, the flight time t, — t needed to reach 
a point of the set (r, gi) occurs explicitly. The integrand 
F in [12] depends on ¢, r(¢), 11(¢), 11’ (¢), v2(¢), v2"(¢), o(¢), 
and a(g). It is clear that the time also occurs implicitly in 
the velocity v(¢), since the functions r(g) and v(¢) together 
completely describe both the geometrical and chronological 
motion of the vehicle. In order to obtain from Equation [12] 
a necessary condition for the minimum flight time, however, 
one can reason in the following manner: For arbitrary given 
flight path characteristics, which are represented by ¢, r(¢) 
and v(¢) [the velocity components »;(¢) and v2( ¢) are uniquely 
determined by r(g) and v(¢)], the integral in Equation [12] 
is a minimum only if the Euler equation with respect to a is 
satisfied. Since this is true for arbitrary r(g) and v(¢), it is 
also true for the optimum flight path, if it exists. In the 
present case, the Euler equation with respect to a is 


oF 
= 


=0 


+ Avs) — — Bos) | [13] 


re) 
da {(Av, — Bre) sin — (Bo, + Ave) cos ¢} 


Taking into account that a@ only occurs explicitly in the 
quantities A and B, one obtains for [13] 
— Br) sin — (Bu, + Ave) cos g] X 


[1'(Bati + — v2'(Aari — Barz)] — 
[((Agu: — Bate) sin g — (Bari + cos X 


+ Ar) — %(An — Boat =©0 [14] 


where N is the expression in the denominator of F and the 
subscript a denotes the partial derivative with respect to a. 
A suitable collection of the particular terms within the braces 
in [14] leads to 


(v;’ sin — cos — BAg)v? = O 


[15] 


Multiplication by dg/dt # 0 leads to the minimum condition 


sin g — cos 
y[(Avu — Bre) sin g — (Bu; + Avz) cos 


(ABz — BAg) = © [16] 


This condition has to be satisfied along with the equations of 
motion [4a and 4b]. Multiplying [4a] by sin ¢g and [4b] by 
cos g and subtracting the second from the first, one obtains 


sin — cos Y = 


1 


[((Au, — Bre) sin g — (Bu, + Ave) cos ¢] 
mv 


Equation [16] then can be written as 


(ABa — BAa)v? 


=0 
— Bre) sin — (Bu + Ave) cos ¢] 


from which the minimum condition in question becomes 
— BAg = 0 [19] 


The quantities A and B are functions of r, ¢, v1, v2 and @ as 
stated in Equations [9a and 9b]. From [19] one obtains a 
as a function of r, ¢, v; and v2 

a = ar, ¢, 1, [20] 


Replacing a in the equations of motion [4a and 4b] by [20] 
and using for v; and v2 the kinematic relations [2], one obtains 


354 


attack a 


the equations of motion in the form Fey 
— cos — + rg) sin = flr, [21a] 
+ rg) cos + (r — rg?) sin = for, [21b] 


From these equations the functions r(t) and g(t) can be deter- 

mined. Under rather weak assumptions on f; and f2 (Lip- 

schitz condition), system [21a, 21b] leads to unique solutions 

r(t) and g(t) which satisfy the initial conditions at t = to 
% 


r(to) = To = to [22a] 
Y(t) = Go = [22b} 


The entirety of the solutions issuing from the initial point (ro, 
¢o) forms a two-parametric family of curves where 7p and ¢ 
are the parameters. For given fixed initial velocity 


Vint + Po? (23) 


only the initial direction of the flight path can still be chosen. 


Measuring the initial flight path angle Yo from the line 
tangent at the circle with radius ro, this angle is given by 
tan Wo > [24] 

By asuitable choice of Yo a trajectory can be performed which 
_ passes through a given final point (mi, ¢:). Every point lying 
on that particular curve is then reached in a minimum time. 
Finally it may be mentioned that one arrives at the same 
minimum condition, if in Equation [10] the variable r is 


chosen to be the independent variable instead of ¢. ‘ 
The Case of a Flat Earth Wei. ou Pa | 
This case differs from that of a spherical Earth insofar as 
the force of gravity G is pointing toward the negative direction 
of the y axis (Fig. 1). The equations of motion are in the 


Cartesian system of coordinates z, y 


mi; = (7 cosa — 6 — mv) —(rsina +A) [25a] 
mM. + y) = (rT sina + dA) + (7 cosa — 6 —mv) [25b] 


where 


= y 
Vette 


[26] 


v= 


The quantities 7, 6, \ and m are assumed to be arbitrary 
(differentiable) functions of the horizontal distance zx, the 
altitude y, the velocity components 1; and v, and the angle of 


6 
= T(x, Y, M1, U2, @) “pi 
= (2, Y, V1, Ve, a) 

= y, V1, V2, 
m(x, Y, V1, V2, a) 


Y = (yo*g0)/y? 


where go is the gravitational acceleration at altitude y = yo. 
Excluding here the ballistic trajectories which result from 


and dividing one of the equations of motion [25a, 25b] by the 
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other, one obtains : Wet For given fixed initial velocity 
4 + Bu) _, to = Wao? + go? [37] 
a. ee ' and a suitable choice of the initial flight path angle 
Here A and B have the same significance as in Equations [9a vin di ale (38) 


and 9b]. Introducing x to be the independent variable and 
denoting the derivative with respect to x by a prime, it 
follows from [30] that 


+ Ave) — v9'(Av: — Bre) 


Integration with respect to x from 2x9 to 2, leads to the flight 
time t, — to needed to reach a point of the set (2, y) 


dx = dt [31] 


+ At») v2'(Ar, Boz) dz = 
xo — Bre) 


x1 
f Fdz to 
xo 


Yo make the flight time a minimum we can reason in the same 
manner as we did in the previous case. The minimum condi- 
tion is then 


[32] 


oF /da = [33] 
which is equivalent to iA 


a trajectory can be performed which passes through the given 
Every point lying on this particular 


final point (x, 4). 
curve is reached in a minimum time. 
j 


Remarks 


The Euler equations [18 and 34] represent necessary condi- 
tions for a stationary behavior of the correspondirg integrals 
[12 and 32]. As it can be seen, these equations involve also 
the necessary condition for the maximum flight time, namely 
v = 0. If, therefore, the function F in each case is continu- 
ous and not identically constant, the stated conditions for the 
minimum flight time paths are even sufficient. 

It may happen that for different sets of initial values the 
corresponding brachistochronic trajectories will both pass 
through the two given points (2x, yo) and (1%, y:). This, how- 
ever, is not surprising, since the Euler equations in general do 
not furnish unique solutions for extremals “‘in the large.” 
The question for the absolute minimum time trajectory in 
such a case can be answered after determining the flight times 
explicitly according to Equations [12 and 32]. 


— BAg)v?/y[Anu — Br]? = [34] 
Multiplying this equation by dx/dt ¥ 0 and using the equa- So 
tions of motion, one obtains the minimum condition 
1 Miele, A., ‘General Variational Theory of Flight Paths of Rocket 


— BAag = 0 [35] 


which is the same as that in the case of a spherical Earth. 
The only difference is that the quantities 7, 6, A, m and vy 
which occur in A and B have to be replaced by those given 
in Equations [27 and 28]. 

The minimum time flight path in this case is again deter- 
mined by first solving Equation [35] for a and then integrat- 
ing the equations of motion for z(t) and y(t) with proper 
initial conditions at t = ty 
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reactor is presented. 


DETAILED theore sidten analysis of combustion in- 
tensity in an actual combustion chamber in terms of 
elementary transport and chemical reactions will be very 
difficult in view of the tortuous flow patterns in the chamber. 
However, it may be hoped that chamber performance can be 
predicted to a satisfactory degree in terms of simple idealized 
flow patterns or simple combinations of these patterns. 

Theories of chamber performance can, in general, be divided 
into the one-dimensional and well-stirred flow patterns. 
In the one-dimensional or plug flow approach it is assumed 
that a combustible material is passed down a tube and pro- 
gressively undergoes combustion. Combustion of a gaseous 
material for simple plug flow has been analyzed by numerous 
investigators (1)? and for plug flow with recycle of part of the 
burned products by Spalding (2). Combustion of a hetero- 
geneous dust or spray in simple plug flow has also been widely 
investigated (3,4,5). Penner (6) and Hartwig (7) analyzed 
evaporation of liquid droplets during simple plug flow in a 
rocket motor, assuming that a droplet was surrounded by the 
hot gas products during its residence in the motor. 

In the well-stirred reactor approach (8) it is assumed that 
fresh material is instantaneously and perfectly mixed with a 
volume of burning material, and that the composition every- 
where in this volume is the same as the material leaving the 
volume. Homogeneous combustion of a gas in a single 
stirred reactor was analyzed by Longwell et al. (9) and others 
(10,11), who related blowout conditions to combustion 
kineties.* 

Two stirred homogeneous reactors in series were discussed 
by DeZubay (12). Hottel et al. (13) analyzed two stirred 
homogeneous reactors in parallel with part of the products 
from the first reactor being fed into the second reactor, and 
the products from the second reactor being mixed with fresh 
combustible gas and fed into the first reactor. An infinite 
series of well-stirred homogeneous reactors with fresh com- 
bustible gas added to each reactor was analyzed by Avery and 
Hart (14) and Rosen and Hart (15). 

Although the stirred reactor scheme seems intuitively to 
resemble closest the chaotic conditions in an actual com- 
bustion chamber, there appears to have been no theoretical 
analysis of combustion in a heterogeneous stirred reactor, 
where combustion intensity depends upon the burning rate 


~ Presented at the ARS Semi-Annual Meeting, June 8-12, 1959, 
San Diego, Calif. 

1 Senior Scientist. 

* Numbers in parentheses indicate References at end of paper. 

’ Spalding (2) surmised that the actual flow patterns in these 
experimental studies involved turbulent plug flow surrounded b 
an envelope of recirculating burned gas rather than true et 
stirred flow. 
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geneous Stirred Reactor 
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< ‘e droplets in the feed spray. 
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WELBY G. COURTNEY! 
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Richmond, Va. 


A eo analysis si piianibbesitun intensity of liquid monopropellant in a heterogeneous stirre: 
The argument consists mainly of deducing the steady-state concentratio: 
of partly burned drops with radius r which remain in the reactor when a polydispersed spray is fe: 
into the reactor. The resulting equation is intractable, but an approximation to this equatio: 
Se et @ shows that combustion intensity increases with increasing density of burned gas in the reactor 
increasing burning rate of a single droplet, and particularly with increasing number of smaller siz 
Quantitative prediction of combustion intensity is limited by sever 
une ertainty in the aera rates of droplets under the > turbulent rea reactor conditions. 


of a liquid droplet or a solid particle. A theoretical analysi 
of the combustion of liquid monopropellant in a heterogeneou 
stirred reactor is presented herein. 7 


Theory 


Consider a liquid monopropellant combustion chambe 
which is operated as a well-stirred heterogeneous reacto1 
Liquid monopropellant is fed into the chamber as a spray, 
which is instantaneously mixed with the chamber contents. 
A flame envelope forms around a droplet, and the droplet 
progressively burns. A proportion of the chamber contents 
(gaseous products plus partially burned droplets) then 
passes through the exhaust nozzle. We will assume: 

1 The rate of combustion in the chamber is limited by the 
rate at which a droplet burns. 

2 Essentially all drops leave the chamber before being 
completely burned. 

3 There is no droplet agglomeration or fracture in the 
chamber. 

4 Steady state. 

The mathematical analysis given here often follows the 
discussion of nucleation and growth of solid particles in a 
stirred reactor given by Bransom et al. (16). 

Let a monopropellant liquid with mass flow rate Wo gm 
per sec per unit chamber volume be fed as a spray into a de- 
composition chamber of volume V liters. Let the concentra- 
tion distribution of the unburned spray be mo(7ro), where 
No(ro)dro is the number of unburned spray drops per liter of 
chamber volume whose radii are between 7» and 7 + dro. 
The rate of addition of these drops is thus %(7o)dro, but for 
present purposes the derivative with respect to time may be 
neglected since the steady-state drop population in the 
chamber does not change with time. At a later time some of 
the drops added earlier will have been lost through the 
nozzle, and those which have remained in the chamber, n(70), 
will have been partly burned. The rate of loss of these drops 
through the nozzle is 


fires 


6 = age of drop in the chamber, sec 
v = rate of flow of burned gas from the chamber, liters per sec 


Integration of Equation [1] gives 


—v0/V [2] 
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n(ro) = 


y 
<7 
> 
» 
< 
» = 
ay 


where n(ro)d@ is the steady-state concentration of burned 
drops whose initial radii were between ro and ro + dro and 
which have age between 6 and 6 plus dé. 

Investigations of the burning behavior of single drops indi- 
cate an induction period for the formation of a flame envelope 
around a drop and then steady burning of the drop. Nishi- 
waki (17) concluded that the induction period decreased with 


decreasing drop radius. The specific value of the burning 
constant k (and probably the induction period also) depends 
upon the temperature (18), pressure (18) and particularly 
the turbulence (18,19) of the gaseous environment. For the 
present illustrative purpose we will neglect the induction 
period and assume that k is a constant, or 


—d(r?)/d0 =k 


dr/d@ = —k/2r 

6 = (ro? — r*)/k [5] 
Converting n(7o)d6 to n(ro)dr by Equations [2, 4 and 5], and 
multiplying by —1 (to make dr go from 0 to 7 instead of ro to 
0) gives 


[3] 


where n(ro)dr is the steady-state concentration of burned 
drops with radii between r and r + dr, and whose initial 
radii were between 7» and rp + dr. The total concentration 
of burned drops N(r)dr with radii between r and r + dr 
includes, for example, those with small initial diameter and 
short residence time plus those with large initial diameter 
and long residence time, or 


2r-explor?/kV]- exp[—vro?/kV] -no(r0) 
min 
k 


(7] 


Nr) = 


where fomax 2Nd 7omin are the maximum and minimum drop 
radii in the initial spray. 
The steady-state weight of drops per liter of chamber vol- 
cham! 


ume is 
W = f -dr [8] 
Tmin 
where 
wr) = N(r)/3 
p: = density of the liquid drops which is assumed to be 


constant 


The lower integration limit is the radius of the smallest drop- 
let in the chamber and can probably be replaced by zero 
without introducing appreciable error. Therefore 


3k J0 
dry |-dr {9 
exp kV dro r [9] 


Combustion intensity € in grams of monopropellant burned 
per liter of chamber volume per sec is equal to the mass rate 
of flow of monopropellant into the chamber minus the mass 
rate of loss of droplets from the chamber, or — : 


e=Wo- Wn 


v = [11] 
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where 
W, = rate of flow of gas from the chamber, gm per sec 
Pg = density of the outflowing gas 
Material balance gives e=W, [12] 


Finally, combustion intensity is obtained by substituting 
Equations [9, 11 and 12] into Equation [10], giving 
Wo 


871, 10 €r °max €ry [13] 
1+ r*-exp exp (- | dr 


Equation [13] appears to be intractable for explicitly relating 
combustion intensity to the other parameters, but could be 
solved numerically if the values of the parameters were 
known. Unfortunately, k is unknown and there would seem 
to be little hope of guessing a reasonable value of k on theo- 
retical grounds in view of the uncertain but undoubtedly high 
degree of turbulence in the combustion chamber. An experi- 
mental investigation with a heterogeneous stirred reactor is 
required to investigate k, and even then the value of k will 
probably vary with combustion intensity. 

However, if we assume for convenience a monodispersed 
spray with initial drop radius ro, the integral in Equation [13] 
is dominated by r4, and Equation [13] may be approxi- 
mated by neglecting the exponential factors, giving 


W, 


1 + 2ro2Wo/5kpy [14] 


In Equation [14], combustion intensity increases as k or 
Po increases and particularly as 7> decreases. Thus, com- 


bustion intensity with a polydispersed spray would increase 
with increasing number of smaller size droplets. 
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Technical Notes 


Estimation of Nonequilibrium Reaction The characteristic chemical time consistent with Equation [3] 
Flight Regimes for Blunt Bodies at 

T. C. ADAMSON Jr.? — 

University of Michigan, Ann Arbor, Mich ky = 
N IMPORTANT question which arises in the considera- ad 

tion of aerodynamic heating problems is that of ascer- where X is any third body and the subscript 2 refers to oxygen 
taining the flight regimes in which nonequilibrium reaction atoms. (w) recombination is that part of the reaction rate 
occurs. A general idea of this regime may be gained by an which holds for the recombination of oxygen atoms, since 
order of magnitude analysis where the characteristic time of this is the essential reaction in the decreasing temperature 
reaction is compared to the characteristic residence time, with field. It is calculated using the concentration of O and O, 
limiting values of the ratios obtained, arbitrarily assigned. existing behind a normal shock, but at the wall temperature. 
Moreover, such an order of magnitude analysis is useful in This is a very conservative calculation, since an average tem- 
those cases where the flow is very far from equilibrium; and perature between the wave and the body would certainly be 
recombination may take place, finally, along the surface, since higher than the wall temperature, and the recombination rate 
an idea of the wetted area of reaction may be obtained. This term is sensitive to changes in temperature. Since no accurate 
type of study is of value in those cases where either catalytic value for the average temperature was known and since it is 
or reaction poisoning surfaces are being considered. evident that a linear variation is unrealistic, the conserva- 
In this analysis, the flow field is assumed to be that around tive value for the reaction rate constant was chosen. For a 


a blunt body at hypersonic speeds such that a detached shock 
exists. The field is divided into two separate regions, each 
with its own characteristic time. One is the stagnation point 
flow between the detached shock and the body. The other is 
the laminar boundary layer flow along the body surface. Each 
region has a characteristic residence time which is to be com- 
pared with the characteristic chemical reaction time. 
According to Li and Geiger (1) the detachment distance 
between a detached shock and the stagnation point is, as in 


Boundary 
Layer 


Fig. 1 
-14+ V1 (1 — k)? 
kRsy (1) Blunt Body 
aty (1 — k) 
s = conditions immediately behind the shock 7 re 
R,, = radius of curvature if the shock at z = 0 Pt eel 


If, as assumed by Li, the radius of curvature of the shock is 
the same as that of the body, at z = 0, then b may be cal- 
culated. Also, from his calculations, the velocity in the y direc- 
tion outside the boundary layer is essentially proportional to 
y; since the v velocity in the boundary layer is well approxi- 
mated by a linear function, the average velocity in the region 
between the detached shock and the body along the stagnation 
streamline is, to good approximation 


Vavg = v,/' 2 [2] 


The residence time in this region is thus eae 
= 2b/v, [3] 


Original received Nov. 18, 1958 and revised Feb. 1, 1960. 

1This work was supported by Air Force Contract no. AF 
18(600)-1607 under a subcontract with Bell Aircraft Corp. . 

? Associate Professor of Aeronautical and Astronautical Engi- Fig. 1 Stagnation point region between a detached shock and a 
neering. blunt body 
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Fig. 2 Variation of time ratio function with altitude with flight 
Mach number as parameter for blunt nose body with radius of 
curvature Fp, ft, = 1500 K 


more complete comparison, one could use the values of the 
temperature behind the normal shock for another limit, al- 
though this was not done for the present calculation. 
From Equations [1, 2 and 3], the ratio of characteristic 
times R, can be written as 
ty 
Us (Ke)s 


1 


where relative mass fractions have been introduced for the 
mole fractions in the rate term. According to Hirschfelder 


[6] 


(2) 
5/2 ems 
k, = 1. 16% | — ———— 7 
(=) (moles)? sec (7] 


and as noted, is calculated using T'y. 


Using values given in (3 and 4), R, Ry, was computed and 
plotted as a function of altitude with Mach number as a 
parameter (Fig. 2). In these calculations, 7, was taken to be 
1500 K. It is interesting to note that according to Fig. 2, 
as the Mach number increases at a constant altitude Rp,R, 
decreases. Thus the characteristic chemical time decreases 
faster than the characteristic residence time. Although the 


50 100 150 200 250 300 
ALTITUDE IN 10> FEET 


Fig. 3 Variation of time ratio function with altitude with flight 
Mach number as parameter for downstream distance along body 
surface L ft, = 1500 K 


creases as the square of the pressure, and this effect is stronger 
than the velocity increase. Since the temperature behind the 
wave increases along with the Mach number, it is clear that 
the use of some average temperature behind the wave, rather 
than the wall temperature, would cause R, to decrease even 
more as the Mach number increases, if the reaction rate is 
proportional to T°”. 

For any arbitrary limiting value of R, with a given nose 
radius Rp,, the range of flight Mach numbers and altitudes can 
be found, below which the stagnation flow can be considered 
to be in chemical equilibrium. For example, if R, < 10? is 
chosen as a limiting value, then for a blunt body of 1-ft nose 
radius at flight Mach number 11 the flow is considered to 
be in equilibrium for altitudes less than 135,000 ft. As the 
altitude increases the existence of chemical equilibrium be- 
comes less and less probable until nonequilibrium reaction re- 
sults. Since this was a conservative calculation, the altitude 
in the example noted must represent the lower limit for non- 
equilibrium reaction. 

R, may also be used to find the approximate conditions at 
which surface reaction can take place. Thus when R; > 1, 
the oxygen atoms could reach the body surface before react- 
ing. In order to calculate the order of magnitude of the pos- 
sible wetted area of reaction, another characteristic residence 
time, the boundary layer residence time, must be compared 
with the chemical reaction time. In this case, a convective 
term and the reaction term in the species continuity equation 
may be compared to obtain the desired time ratio. 
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ty (We2)recombination [8] 


This ratio is then the mass of oxygen atoms convected per 
unit volume per unit time to the mass of oxygen atoms dis- 
appearing by recombination per unit volume per unit time. 
It should be noted that the diffusion terms, in the boundary 
layer equations, are the same order as the convection terms, 
so that this ratio should hold for the time associated with 
the cross stream diffusion of mass as well as the time 
associated with the convection of mass, as long as the Schmidt 
number u/pDis of order 1. If ~ (kz),/L and u ~ »,, 
then substituting for w. as in Equation [6] one can calculate 
R.L in asimilar manner. Thus 


1 v,(Ke)s 
"\Ro) Tw? + 1/, 


Fig. 3 contains the plot of R,L vs. altitude with Mach num- 
ber asa parameter. As in the previous calculation, the values 
of k, were approximated by those existing behind the normal 
shock, whereas the same temperature was used to calculate 
k,. 


Fig. 3 may be used to obtain an order of magnitude of the 
possible atom wetted area for these cases where (according to 
Fig. 2)R, > 1, i.e., where the characteristic chemical time is 
greater than the residence time, and hence atomic species 
enter the boundary layer. For example, for a blunt nose body 
of 1-ft nose radius traveling at Mach number 11 at 225,000- 
ft altitude, R, is 10, according to Fig. 2, so that atomic species 
enter the boundary layer. To calculate the order of the dis- 
tance back along the body where one might find atomic 
species, one finds from Fig. 3 that at 225,000-ft altitude, when 
M = 11,R,Lis0.85. Therefore < 1 if L > 0.85 ft; 
the atom wetted area extends downstream to the order of 
0.85 ft, so that surface reaction could extend back to this 
distance. 

It is evident that the same type of calculation could be 
carried out to find the conservative downstream limit of the 
region of nonequilibrium flow for those cases where R;, is 
greater than the arbitrarily assigned value of 10~? in the stag- 
nation point calculation. That is, if the initial computation at 


the stagnation point indicated that R, > 10~* the downstream 
distance L to the point where R, = 10~* could be calculated 
using Fig. 3. 

It should be emphasized that the values obtained in Figs. 2 
and 3 depend entirely on the chemical parameters used in the 
reaction rate term, and since there is some argument as to 
what these values should be, the results of these analyses are 
inconclusive insofar as absolute values are concerned. How- 
ever, they do illustrate a method of calculation which should 
become of more quantitative value when more accurate 
chemical parameters are available. 


Nomenclature 

b = distance between body and detached shock along stagna- 
tion streamline 

ky = forward specific reaction rate constant 

k, = reverse specific reaction rate constant 

L_ = distance along body surface downstream from stagnation 
point 

m = molecular weight 

Nz = moles of atomic oxygen per unit volume eG 

N = moles of mixture per unit volume 


P = pressure 

Ry = radius of curvature of blunt nosed body 

Ry = universal gas constant 

R, = ratio of characteristic chemical time to characteristi: 
residence time 

T = temperature 

w:, = rate at which mass of atomic oxygen is formed by chemi 
cal reaction per unit volume per unit time = 


p» = density of atomic oxygen ~ 
ke = relative mass fraction of atomic oxygen - 
Subscripts 

0 = at stagnation point 

s = immediately downstream of detached shock 

w = at wall 

o = flight conditions 
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Sloshing of Liquids in Cylindrical Tanks 
of Elliptic Cross Section’ 


WEN-HWA CHU? 


Southwest Research Institute, San Antonio, Texas 


An analysis is made of the sloshing of liquids in partially 
filled rigid cylindrical tanks of elliptic cross section. By 
application of potential theory, the forces and moments 
acting on the tank are determined for the case of trans- 
lational and rotational excitation along or about the major 
and minor axes of the cross section, and expressions for the 
resonant frequencies are found in terms of the roots of the 
derivatives of modified Mathieu functions. The lowest few 
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resonant frequency parameters, as computed from those 


expressions, are compared with those of a circular cylinder 
having the same cross-sectional area. 


LOSHING motions of fuel in partially filled tanks of mis- 
siles may be established by essentially lateral accelera- 
ions of the vehicle. Such liquid motions are of importance, 
for there is the possibility of extreme oscillations if the fre- 
quency of excitation is in the neighborhood of one of the na- 
tural frequencies of the liquid fuel. The forces exerted by the 
fuel on the tank walls may then lead to perturbations in the 
flight trajectory, or may even impose severely high stresses 
on structural components. 

The force and moment resulting from free oscillations of 
the fuel in a rectangular tank have been presented in (1),° 
and those resulting from horizontal (translational), pitching 
(rocking) and yawing oscillations of the rectangular tank are 
given in (2). The free oscillations of fuel in tanks of other 
cross sectional shapes are also given in (1). The same prob- 
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lem for circular tanks undergoing horizontal and rotational (13) is presented in this paper. Although only a limited 
oscillations is treated in (3 through 5). Sloshing in flexible amount of computational effort has been carried out, it is 
tanks, considering such elastic vibrations as breathing, bend- possible to infer some of the dynamic effects resulting from 
ing and rotationally symmetric modes, are treated in (6 even a small ellipticity. The theory is not, however, re- 
through 8). stricted to tanks with small ellipticity. 


The tank of circular cross section is simple and is generally 


employed in missiles; however, some out-of-roundness may Force and Moment Acting on the Tank 


exist in practice or in laboratory models, and it is conceivable The tank is assumed to possess translational oscillations 
that tanks of elliptic cross section may even be desired for x», y» and rotational oscillation 6, 6,, all of small amplitude. 
certain applications. It has been noted (9) that the wave The incompressible potential problem (3, 4, 5, 13) is solved 
motion can become very complex when the variation in tank by separation of variables and expansions analogous to 
diameter is on the order of 0.4 per cent or greater. Conse- those in the problem of vibration of elliptic membrane (10). 
quently it is of some interest to investigate analytically the The forces and moments are given as follows 

sloshing of liquids in elliptic tanks. The result of the analysis de 
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and 
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2 sinh 


where 


Try 1?(&) = Sen(E, sinh (lE)dE [4b] 
In the limiting case of zero ellipticity (e > 0, 3ce* > 1), 
the previous results agree with those of a circular cylindrical 


tank (4,5). 
Numerical Examples 


In order to obtain some idea of the change in resonant fre- 
quencies, two elliptic cylinders of 1 and 5 per cent out-of- 
roundness, respectively, are considered, the out-of-roundness 
being defined by (a — b)/(a + b). 

The lowest six parametric zeros entering the general prob- 
lem of these two cases are calculated by making use of the 
available asymptotic formulas (10) and are given in Table 1 


= {1 2 cosh gMh9, E + 


\ tanh — 1 + sech 


Linn® sech (fimn®h) 


[1 — 2 cosh Mh sinh £5 — x 


—— 


= | [4a] 
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Pid j 


Define a frequency parameter N = w?R/az, where R is 


the radius of a circular cylinder of the same cross-sectional 


area. The resonant frequency parameters are then 
bh 
Ninn = Amn lo tanh | Amn 
\ a a 


(Amn = 2V ann cosh &) 
As h/R— o (or sufficiently large), the lowest six resonant 
frequencies of the given examples are compared with the 
lowest two of the equivalent circular cylinder in Table 2. 
Discussion 


It is noted that although the calculated values of Amn 
and Amn‘? are higher than the corresponding \,,,, the reso- 


Table 1 Parametric zeros 
% (1a) (a — b)a + b) = 0.01 (1b) (a — b)/(a +b) = 0.01 
qu® = 0.03334 qu® = 0.03460 qu = 0.15743 qu® = 0.18678 
qu = 0.28143 = 0.28740 = 1.3630 = 1.5144 
qa = 0.17648 3 = 0.17648 qu? = 0.87887 qu® = 0.88058 


Table 2 Comparison of resonant frequency parameters h/R — © 


(2a) (a — b)/(a + b) > 0.01 i (2b) . (a — b)/(a + b) = 0.05 (a — b)/(a+b) =0 
D E 
Ny™ = 1.8259 Nu = 1.8599 Ny™ = 1.7569 Nu = 1.9303 = 1.841 = 
Ni; = 4.2008 Nix = 4.20088 Nis = 4.1872 Nis = 4.1912 4.201 — 
N2© = 5.3047 Na? = 5.3607 Ny = §.2145 Nx® = 5.4964 N2 = §.332 = ry 


. = Horizontal motion along the major axis or rotation around the minor axis 
= Horizontal motion along the minor axis or rotation around the major axis 
C = Same as A 
D = Same as B 
EK = Horizontal and rocking oscillation of a circular cylinder of the same cross-sectional area 
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nant frequency parameters Ny << < Ny and No < 
NO < No. Furthermore, there are additional resonant fre- 
quencies introduced, due to the small ellipticity, correspond- 
ing to modified Mathieu functions of odd integer order higher 
than one. The numerical values of Ns: and N2; are both 
slightly less than 

With some positive empirical damping (12), the maximum 
magnitude of the additional terms as well as the dominant 
terms can be calculated [see (13)]. In practice, it is often 
sufficient to consider only the lowest few modes. 

To obtain some idea of the amplitudes, assume qmn small 
and h— 0; then, the coefficient of 1/[(wm,/w)? — 1] of the 
additional modes due to ellipticity is roughly of the order 


= =0 
A,O(qu™) 4 ) 4 
!(n — 1)! 


as compared with the coefficient of the fundamental mode, 
or of the order 


BY (mn) 0 
4 


In the case of 1 per cent out-of-roundness, the values of 
ns“? /4 and q13/4 are slightly more than 4 per cent; therefore 
the actual amplitude of the lowest additional modes may be 
not really serious near the corresponding resonant frequency 
parameter. The case of larger qmn cannot be easily discussed. 
The qualitative conclusion remains to be confirmed by ex- 
periments. 
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Nomenclature 

2a = major axis of the elliptic cylinder 4 

2b = minor axis of the elliptic cylinder ‘- - 

c = location of the foci of the ellipse 

q = magnitude of the gravitational acceleration 

h = height of the undisturbed liquid level 

= mth parametric zero of Cen’(, q) = 0 

= mth parametric zero of Se,’(%,q) = 0 

a ae = tank fixed Cartesian coordinates with origin at the 
center of the tank bottom (z upward along the 
tank axis) 

Lb, Yo = z- and y-components of the small displacement of 


the cylinder, respectively, with respect to Earth 
fixed coordinate system 


An™ = coefficient of cos (mn) of the expansion of cen(n, q) 
B,™ = coefficient of sin (mn) of the expansion of se,(n, q) 
Ce,(€, q) = even modified Mathieu function of order n 
Ce,’(¢,q) = first derivative of Ce,(é, q) with respect to & 

cen(n, gq) = even Mathieu function of order n 

F,,F, = 2x- and y-components of force acting on the tank 
Lana = coefficient of the expansion of cosh £ cos 7 in modi- 


fied Mathieu functions, defined by Equa- 
tions [1b, 1c] 
Linn = coefficient of the expansion of sinh — sin » in modi- 
fied Mathieu functions, defined by Equations 
[2b, 2c] 


M,,M, = moment about z and y axis, respectively, acting on 
the tank, due to liquid sloshing 

M = total mass of the fuel = rpabh 

Ninn = (Wmn)*R/er, resonant frequency parameter 

Np» = mth resonant frequency parameter of a circular 
cylinder of radius R 

R = V ab, the radius of circular cylinder of the same 
area 

Se,(é,q) = odd modified Mathieu function of order n 

Se,'(é, q) = first derivative of Se,(t, q) with respect to 

sén(n, q) = odd Mathieu function of order n 

&,n, z = elliptic cylinder coordinates 

Or, Oy = small, angular displacement of the tank about the 
z’ and y’ axis, respectively, measured from the 
vertical upward direction (z’, y’ are Cartesian 
coordinates translating but not rotating with 
tank) 

p = liquid density 

w = frequency of oscillation 

ssa = resonant frequency corresponding to gmn 

= ar tanh (umaPh); = 1,2) 

Aman = 2Vqmn cosh 

Amn = mth zero of the first derivative of the Bessel function 
of the first kind of order n, i.e., Jn’(Amn) = 0 

= separation constant corresponding to 

ar = axial acceleration of the tank in the absence of 

= gravity, i.e., axial acceleration minus axial com- 

; ponent of gravity 

Subscripts 

= first time derivative of ( ) 

= second time derivative of( ) 


( )}© = quantities related to the cylinder of radius R, unless 
defined otherwise 

( )® = quantities related to gm, 

( )®= quantities related to 

( )» = for quantities related to the boundary © 
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Simple Formula for Prediction and 
Automatic Scrutation 


R. J. DUFFIN' 
Duke University, Durham, N. C. 
T. W. SCHMIDT? 


U.S. Army Office of Ordnance Research, Durham, N. C. 


Prediction from data obtained at regular time intervals 
is discussed. Attention is restricted to the case in which 
the more recent data should be given more weight than 
past data. Thus a polynomial is fitted to the data by ex- 
ponentially weighted least squares. It is found that the 
formula for predicting the next datum is quite simple. 
This formula needs to “‘remember”’ only a limited amount 
of data. The prediction formula can be used to auto- 
matically scrutinize data for gross error. A gross error 
or a gap in the data may be replaced by a predicted value, 
and the prediction and scrutation process can be con- 
tinued. 


HIS paper deals with data obtained at regular time inter- 

vals. From a knowledge of a regular sequence of such 
observation, it is desired to predict the next observation. 
This prediction is carried out here by use of a least squares fit 
of the data with a polynomial. 

We are concerned with a situation in which observations in 
the distant past should be given less weight than those more 
recent. More precisely, it is assumed that the observations 
should be discounted by a certain percentage in every time 
interval. This is taken to mean that the least squares 
formula has an exponential weight factor. It is found that 
this method of discounted least squares lends itself nicely to 
algebraic analysis. In particular, the constants of the pre- 
diction formula can be explicitly evaluated. 

It is supposed that the sequence of observations is infinite, 
because this actually simplifies the formal considerations. On 
the other hand, the resulting prediction formula is a linear 
expression in the infinite sequence of observations. Such an 
infinite series is not too convenient for numerical evaluation. 
It proves possible, however, to obtain a finite formula which 
is essentially equivalent. Thus suppose the fit is to be given 
by a polynomial of degree less than m. Then the prediction 
is given by a simple formula involving only the last m ob- 
served values and the last m predicted values. 

By “scrutation of data’”’ we mean rejection of obvious errors 
or blunders. In hand calculation this is often accomplished 
by visual inspection of a graph of the data. The trend in data 
processing is toward elimination of human judgment. The 
prediction formula can be employed as an automatic scruta- 
tor, to this end. 

In automatic scrutation, predicted and observed values are 
compared. If the “error” exceeds a prescribed tolerance, the 
observed value is rejected and is replaced by the predicted 
value. This process can be used also to replace gaps in ob- 
servational data. The programming of the prediction and 
scrutation formula for a digital computer is straightforward. 
The method and principal formulas are given in this paper. 
A more complete mathematical analysis and proof will be 


Let y1, Y2, Ys, .. . be a sequence of real ‘It may 
supposed that these numbers are a time series of observations 
and that y, is the observation at time t = —n. The develop- 
ments considered in this paper stem from a certain approxima- 
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Table 1 Straight line predictors 


n Q(n) n Q(n) 

1 0.4000 0000 51 0.0000 2997 

2 0.2800 0000 52 2455 

3 0.1920 0000 53 2010 

4 0.1280 0000 54 1644 

5 0.0819 2000 55 1345 

6 0.0491 5200 56 0.0000 1099 

7 0.0262 1440 57 898 | 
8 0.0104 8576 733 
9  @.0000 0000 59 
10 0.0067 1089 60 489 

11 0.0107 3742 61 0.0000 0398 


| 
12 128 8490 62 | —— 825 
13 137 4390 63 265 
14 137 4390 64 216 
15 131 9414 6 — 176 
16 0.0123 1453 66 0.0000 0143 
17 112 5900 67 117 | 
18 101 3310 68 95 
19 90 0720 69 77 
20 79 2634 70 63 
21 0.0069 1753 71 0.0000 0051 | 
22 59 9519 72 
23 51 6509 73 ay — | 
24 44 2722 
25 37 7789 75 22 
26 «0.0032 1121 76 0.0000 0018 
27 27 2008 77 15 
28 22 9696 78 a 12 
29 19 3428 79 Me 10 
30 16 2480 80 8 
31 0.0013 6173 81 0.0000 0006 
32 11 3890 82 5 
33 9 5074 83 4 
34 7 9228 84 3 
35 6 5918 85 3 
36 0.0005 86 0.0000 0002 
37 4 5433 87 
38 3 7644 88 1 
39 3 1154 1 
40 2 5754 l 
41 0.0002 1268 0.0000 0001 
42 1 
43 1 1 
44 1 0.0000 0000 
45 
46 0.0000 
47 
48 
4 

Long formula: m = 2, 6 = 0.8, yo* = * Q(n)yn 
n=1 
Short formula: m = 2, 6 = 0.8, y* = 0.4m + 1.6y* — 
0.36y2 — 0.64y2* 


tion of the observations y, by a polynomial 


m—1 


The degree of plz) ts is taken to be less than a given integer m. 
As a matter of notation, it is convenient to write p, if n is an 
integer. 

A measure of the error of the approximation is defined to be 


E = — un)? [2] 


n=1 


_ Here the constant @ satisfies the inequality 0<@<1. The 
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least square approximant is taken to be the polynomial p, of 
degree less than m which minimizes E. Because the weight 
function @ is an exponential, an approximation of this type 
can be termed discounted least squares. 

[t is required to choose the coefficients a; of p(x) so that p(x) 
minimizes Thisimplies that 0F/0a, = 0 fork = 0,1,..., 
m — 1. Carrying out the differentiations results in a sys- 
tem of m equations for the m coefficients. Solving these nor- 
mal equations gives the least squares polynomial p(z). 

The polynomial p(x) is used to extrapolate the sequence 
y,. Of main concern is the “next value” yo. The predicted 
value of yo will be defined as p(0). Since p(0) = a, it is 
simply necessary to solve the normal equations for a. This 
yields a formula of the type 
yo* = [3] 


> n=1 


Here yo* denotes the predicted value yo, and Q(n) denotes co- 
e‘ficients which depend on m and @ but not on the sequence 
y.. Weterm Equation [3] “the long formula of prediction.” 

It proves possible to solve the normal equations explicitly. 


Ii we write Q(n) = q(n), then it is found that 


z- 


Here (7) denotes the binomial coefficient. It is seen that 


q(x) is a polynomial of degree m — 1 in x. The polynomial 
g(x) may be termed an associated Laguerre polynomial of 
discrete type because of a certain difference equation it 
satisfies. 


Vinite Prediction Formula 

A practical difficulty with the long formula [3] is that if 0 
is close to unity, a great number of terms in the series must 
be retained. This might well put a burden on the computer’s, 
memory. The procedure described here avoids this difficulty 
because it is found sufficient to store only 2m properties of 
the infinite time series. 

Let y,* be the predicted value of y, in terms of preceding 
values. Then 


= QnYnte [5] 
n=1 


It is found that the following identity holds between the ob- 
served values y, and the predicted values y,* 


y* = >, [yi + — y;)] [6] 


This identity may be employed as a prediction formula, and 
we term it “the short formula of prediction.”” Thus the pre- 
dicted value is given as a linear function of the last m ob- 
served values and the last m predicted values. Of course the 
short formula and the long formula are basically equivalent. 


Aprit 1960 


It isseen that the short formula fory;*is 


m 
= [yite + O(yj+e* — [7] 
fu 

For k = 0 this reduces, of course, to Equation [6]. For k = 
—1 it would appear that Equation [7] cannot be applied 
because yo is not known. However, if we assume yo = yo*, 
then Equation [7] does yield a value for y—,*. It can be 
shown, moreover, that y-1* = p(—1) where p(z) is the poly- 
nomial minimizing Equation [2]. This process can be con- 
tinued to obtain y—.* = p(—2), and so forth. It is thereby 
seen that the prediction formula [7] is not limited to one-step 
predictions. 


Automatic Scrutation 


The programming of the finite prediction formula for a 
digital computer presents no difficulties. It is first necessary 
to select a starting index k = K. Thus y,x* is to be the first 
prediction. The values of yx+i*,.. ., Yx+m* appearing in 
Equation [7] are not known. Let us assume them to be equal 
tO Yx+1,-- respectively. This is essentially the same 
as assuming that the values of y, for n> K are given by a 
certain polynomial. This method of starting has the desirable 
property of giving perfect prediction if y, is actually defined 
by a polynomial of degree less than m. 

In the programming of [7], provision can be made for 
scrutation of the data. Suppose by theoretical or empirical 
means an estimate can be made of the standard deviation of 
Yn, Say (Yn). Then if 


lyn* — yn| > 30(yn) [8] 


Yn is considered to be a blunder and is replaced by y,*. The 
factor 3 in relation [8] is merely a conventional value. 

A common type of excess error is a missing observation. 
For example, it may be that at some periods a missile whose 
path is being observed is obscured by clouds. Thus if y, is 
lost, it is replaced by y,*, and the program proceeds. Pre- 
sumably a fair number of gaps in the data can be bridged in 
this way. 

We have been speaking of prediction as if it were occurring 
in “real time.’”’ However, it may be used to speed certain 
types of data processing. It could be used to locate the 
probable position of a missile in a series of photographs. 

For other applications this automatic program could be 
employed solely asa scrutator. Then the short formula would 
serve to reject blunders and to fill in gaps in a data sequence. 


Specific Example 

To obtain a quantitative feeling for the prediction prob- 
lem, it is of value to consider a specific example. To this end 
we give the simple case of straight line prediction, m = 2. 
We somewhat arbitrarily select 0 = 0.8. Table 1 dis- 
plays Q(x) = (0.8)? (9 — x)/16 to eight decimal places. 
Since > [Q(n) }? = 221/729, the variance of the prediction 
would be about three tenths of the variance in the observa- 
tions. 
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_ Thermally Induced Bond Stresses in 
Case-Bonded Propellant Grains’ 


Pa ERIC E. UNGAR? and 
BERNARD W. SHAFFER® 


New York University, New York, N. Y. 


ey roar are presented for the stresses induced by 
radial temperature gradients in assemblies consisting of 
cylindrical propellant grains which are bonded to thin 
walled casings. Very long and very short assemblies are 
studied under the condition of sudden environmental 
change. It is shown that under otherwise identical condi- 
tions the bond stress in a very long assembly always exceeds 
that in a very short assembly. An optimization study 


indicates the effect of geometry and material propertios 


on the maximum stress in the bond. 


HE PROBLEM of thermal stresses in an infinite, hollow, 

case bonded cylinder was analyzed in detail in a recent, 
very elegant paper by Zwick (1).4 As he pointed out, the 
solution of this problem may be used to approximate the 
stresses in a long but finite propellant grain and case assembly. 
The solution is not directly applicable, however, to the many 
practical designs where the length of the grain is of the same 
order of magnitude as the diameter. 

A completely general analysis, dealing with the thermally 
induced stresses in cylindrical propellant and case assemblies 
of all possible length to diameter ratios, appears to be so com- 
plex as to be infeasible at present. However, the problem of 
very short assemblies may be solved with relative ease. Once 
the stress solution for both very long and very short assemblies 
are available, they may provide a guide for the design of prac- 
tical units. 

This note presents the solutions applicable to very long and 
to very short assemblies, with emphasis on those assemblies 
having casings that are elastically stiff but thermally thin. 
In addition, the stresses in the bonds joining the propellant 
grains to their casings are discussed in detail, with a view 
toward determining design criteria which can be used to 
minimize these stresses. Further, the relation between the 


bond stresses in very long and very short assemblies is es- 
tablished in order to demonstrate the effect of length. 


Thermal Stresses in Very Long and Very Short 
Assemblies 


Let us consider a cylindrical propellant grain of inside 
radius a; it is bonded at its outside radius 6 to a casing whose 
inside radius is also 6 and whose outside radius is c. It is 
assumed that the stresses and deformations, resulting when 
the assembly is subjected to temperature changes, are small. 
enough so that the Hooke’s law may be used in the analysis of - 
the entire assembly. We also assume that the temperature. 
rise everywhere in the assembly (as measured from the tem- 
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perature corresponding to the stress free condition of the 
assembly) is a function of only the radial coordinate r, that 


ae. is, T = T (r). 


Throughout the discussion the subscript P refers to the 


propellant grain, and the subscript C to the casing. The 
bya expressions for the radial stresses o, and the circumferential 


stresses o@ resulting from analysis of very long cylinders 
(generalized plane strain) or from very short cylinders (plane 
stress) may conveniently be written (1,2,3) as 


(o,)p = kp[(1 — + Wr) — Tr] 
= kp{(1 + a*/r?)(Tp + Wr) + Tp — 27) 

= — + We) — To) 

= + + We) + Te — 2T 


where r denotes the radial coordinate, 7 the local temperatur« 
rise, and where 


= 2 
rT or) = T(p)pd 


— = Tr(b) 


Tc(1 — c*/b?) = Te(b) 


so that 7'p and 7c denote the mean temperature of the pro- 


pellant grain and of the case, respectiv ely.6 Further, for 


short assemblies 


and for long assemblies 


= aE/2(1 — p) 


where 
a = linear thermal expansion coefficient } 
E = modulus of elasticity 
= Poisson’s ratio 


Of course, properties of the propellant must be used to 
evaluate kp, and those of the case to calculate kc. Also, 
for short assemblies 


(acT¢ — apT p)/ap 


for long assemblies 


W (acTc — apT p)/ap 7 
where 
1 a? b? — a? 
1 + bec 
[4b ] 


The latter applies for long and short assemblies, provided the 
proper expressions for & are chosen from Equations [3]. 


5 The notation defined in Equations [2], introduced by Zwick 


(1), serves to simplify the analysis. 
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The stress in the bond may be determined by setting r = 6 
in the radial stress expressions of Equations [1] and in the 
auxiliary Equations [2]. One then finds that 

= kp(1 — a?/b?)Wp = kc(1 — c?/b?)We [6] 
which of course agrees with Equation [5], and thus indicates 
satisfaction of radial stress continuity at the propellant case 
interface at r = b. 

Assemblies With Elastically Stiff, Thermally Thin 
Casings 


For the many practical instances where the casing is much 
stiffer than the propellant, as exhibited by the fact that 


Ne > Ep, one may note that 


Ec | c? — b? Ec [ — 6? 
|>1 7 
Ep E A Ep E 7] 


0 that the expressions of Equation [4a] may be replaced by 
the approximations 


Wp = ————_ [8a] 
1 — a2/b2 
where for short assemblies - 
= | + Mp 


for long assemblies 


= 1/(1 — up) 
and the ratio expressed in Equation [5] iy 
We/We< 1 (9) 


Extreme Stress Conditions 


Zwick’s proof (1) of the monotonic variation of temperature 
with time throughout the grain of a propellant and case as- 
sembly with a stiff but thermally thin case, applies equally 
well to assemblies of all lengths, as long as purely radial tem- 
perature gradients are assumed. Consequently the maximum 
stress conditions must occur either immediately upon im- 
mersion of the assembly in a new temperature environment, 
or after a long period of immersion. If the environmental 
temperature is an amount 7> above the temperature cor- 
responding to the stress free assembly, the initial stresses may 
be evaluated on the basis that the temperature throughout the 
propellant grain has not yet had time to change, whereas the 
thermally thin case assumes the temperature of the surround- 
ing environment instantaneously upon exposure. 

The initial stresses may thus be computed by setting 7 = 
Tp = Tp =O in the first two of Equations [1] and in Equation 
[8a], and by setting 7 = T; = Ty in the last two of Equations 
[1], in Equation [2] and in Equation [8a]. The results may, 
in view of Equation [5], conveniently be summarized as 


[(Or)P]emo (1 — a?/r?)N 
= (1 + a?/r?)N 


oa” = E 4 E 


2 


where 


kpT,ac/ap 11] 


1 — (1/2)(1 — a?/b?)n 


and 7 obeys Equations [Sb]. 

Similarly, the stresses occurring after a long period of im- 
mersion may be computed by setting 7’ = Tp = Tc = To in 
Equations [1,2 and 8]. One finds that the final stresses under 
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all circumstances are simply related to the initial stress in 
accordance with the expression 


= (1 — 


i=r,0 


The bond stresses corresponding to the foregoing conditions 
= b in the 


may, of course, be computed simply by ‘setting r 
appropriate radial stress expressions. p. 


In many practical instances the bond between the propel- 
Jant and the case is considerably weaker than either of these 
two components, so that evaluation of the maximum bond 
stress becomes particularly important. Compressive failure 
of the bond is extremely unlikely and need not be considered 
further, but tensile bond stresses may be significant and war- 
rant further analysis. If the environmental temperature is 
greater than that at which the assembly is stress free, then 
T; is positive and the initial bond stress is positive (i.e., 
tensile) in view of the first of Equations [10] and Equation 
[11]. Inspection of Equation [12] then reveals that the 
final bond stress is negative if ap is greater than ac, and that 
it is positive but less than the initial stress if a» is less than ac. 
On the other hand, for negative values of 7) (which correspond 
to conditions where the environment is at temperatures which 
are lower that that of the stress free state) the initial and the 
final bond stresses both are compressive, unless ap is greater 
than ac. In the latter case the final stress is tensile. In sum- 
mary then, the maximum tensile bond stress is 


Bond Stress Optimization 


for T) >0 


for ap > Qc 
7 0 
for ap < for 


((1 — a?/b?)N 


Fbmax tensile = 0 
(1 — ap/ac)(l — a?/b?)N 


where N is given in Equation [11]. 
The designer may naturally raise the question as to what 
may be done to decrease the maximum tensile bond stress 
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that occurs with given environmental conditions without al- 
tering the physical or geometric properties of the propellant 
grain. Inspection of Equations [11 and 13] indicates that 
for positive 7 the maximum tensile bond stress may be re- 
duced by selecting a casing materia] with a lower thermal ex- 
pansion coefficient ac, whereas the final tensile bond stress 
(if it is nonzero) may be reduced by selecting a casing material 
whose thermal expansion coefficient approaches that of the 
propellant. Further, from Equations [6] used in conjunction 
with Equations [4], one may determine that the bond stress 
under all conditions may be decreased by the use of a casing 
having a lower modulus of elasticity, as may be expected 
intuitively. 

The effect of the Poisson’s ratio of the casing material is 
not obvious from inspection of the previous equations, but 
may be determined by differentiation of the bond stress 
expressions of Equations [6] with respect to wc. One finds 
that an increase in yc leads to an increase in the bond 
stress of long assemblies and to a decrease in the bond stress 
of short assemblies. The forementioned derivatives, how- 
ever, are found to have very small absolute values for assem- 
blies with stiff casings, so that the effect of uc on the bond 
stress is virtually negligible. 

Of course, practical assemblies are neither infinitely long 
nor infinitesimally thin, so it is necessary that the previously 
listed results be applied judiciously. In order to compare the 
values one obtains for bond stresses in units which are iden- 
tical except for length and which are subject to the same tem- 
peratures, one may use Equations [6] in conjunction with 
Equations [3 and 8] to determine the ratio of the bond stress 
in an infinitely long, to that in an infinitesimally short, unit 

_ 2 — (1 — + pp) (14] 
7 1 + a?/b? — 2up 


long 


short 


This equation, which is plotted in Fig. 1 as a function of 
a/b for several values of wp, reveals that this ratio approaches 


unity as up approaches zero, but that the ratio increases with 
Mp and is considerably greater than unity for practical values 
of wp. Calculations based on the assumption of an infinitely 
long assembly thus always lead to bond stress values which 
are greater than those calculated for a very short assembly. 


Concluding Remarks 


It has been demonstrated that the bond stress induced in « 
very long assembly always exceeds that induced in a simila: 
but very short assembly under identical conditions; this fact 
may be a useful guide for design purposes. 

It is interesting to note also that a number of conclusions 
were found to be equally valid for very short and for very lon 
assemblies. These conclusions pertain to assemblies wit! 
elastically stiff but thermally thin casings and are: 

1 All thermally induced stresses attain their maximun 
(absolute) values either immediately when the assembly i 
subjected to a new environmental temperature, or else afte 
a long exposure. 

2 The bond stresses obtained in a given assembly, whe: 
the grain configuration and propellant material cannot b 
changed, may be reduced by selecting a case which: I: 
thinner, has a smaller modulus of elasticity, and has a smaller 
linear thermal expansion coefficient (to reduce the initia 
stress) and one more nearly equal to that of the propellant 
(to reduce the final stress). 

3 The bond stresses are virtually unaffected by change- 
in the Poisson’s ratio of the casing material. 
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ie IS well known that the lengths of hypersonic nozzles 
contoured to produce uniform test section flows are often 
prohibitive. In such a nozzle the flow along the centerline 
very rapidly expands to the desired test section conditions, 
and the remainder of the nozzle is devoted to bringing the rest 
of the flow to these conditions. The rate at which the uni- 
form flow region widens with nozzle length is, of course, given 
by the tangent of the test section Mach angle. For hyper- 
sonic Mach numbers, where the Mach angle is very low, this 
section of the nozzle constitutes a major portion of the total 
nozzle length. Some illustrative values are shown in Fig. 1 
for a family of axisymmetric nozzles? computed by WADC.* 
It is evident that elimination of a significant fraction of this 
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2 Nozzles whose initial contour downstream of the throat is 
taken as conical and in which pure source flow is assumed. 

* Cresci, R. J., “Tabulation of Coordinates for Hypersonic 
Axisymmetric Nozzles, Part I—Analysis and Coordinates for 
Test Section Mach Numbers of 8, 12, and 20,” WADC Tech. 
Note 58-300 Part I (ASTIA 204 213), Oct. 1959. 


Short Hypersonic Contour Nozzles 
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terminal portion of the hypersonic nozzle represents a sig- 
nificant saving in total nozzle length. 

If a test section flow can be tolerated in which the central 
“core” is uniform and the remainder is slightly nonuniform, 
then such a reduction in total nozzle length can be achieved. 
Such a nozzle would be determined in the following manner. 
For a given contour, choose the desired uniform core diameter 
(see Fig. 2). From this diameter on the exit characteristic, 
Point A, trace along the upstream characteristic to the nozzle 
wall, Point B. This is the nozzle contour cutoff point, 
downstream of which any deviations of the wall from the 
original contour will not be felt inside the test section uniform 
center core provided these deviations do not initially (i.e., 
immediately downstream of the contour cutoff point) produce 
a wall more concave than the original contour. If greater 
concavity is permitted somewhat further downstream of the 
cutoff point, it cannot be sufficient to produce a shock that 
would intrude upon the model flow field. If the axial location 
of the chosen core diameter on the exit characteristic is taken 
as the nozzle exit, the length of the shortened nozzles for the 
WADC contours may be readily calculated. Ratios of these 
lengths to the full nozzle lengths are presented (Fig. 3) as a 
function of Mach number for various initial wall angles. 
For comparison, Fig. 3 of Cresci’s paper is reproduced in 
Fig. 4. It gives the full nozzle lengths in terms of exit or test 
section diameters. The savings in nozzle length shown in 
Fig. 3 are noteworthy. 
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on the wall shape between the exact contour cutoff point and 
the exit of the shortened nozzle. Without solutions of the 
flow in specific nozzle shapes, it is not possible to state with 
quantitative precision the extent of these nonuniformities. 


WADC nozzles 


Fig. 1 Length™of final Mach cone as fraction of total length 
for WADC nozzles 80 
60 
ac: 
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side the chosen uniform core. This will, of course, depend ind. MACH NUMBER 


Fig. 3 Length of shortened nozzle as fraction of full length for 


However, for a test section uniform core diameter of 0.4 times 
the test section diameter, or greater, an examination of the 
Mach 12 and 20 WADC contours and consideration of ap- 
propriate Mach angles strongly suggests that the worst case 
(Mach 12, initial cone angle of 12 deg) would be a cone- 
frustum of about 3.5-deg wall angle between the exact con- 
tour cutoff point and the exit. In such a nozzle, the flow 
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angle of the outer six-tenths of the test section flow would 
vary from zero to about 3.5 deg at the wall, which is a con- 
siderable improvement over a conical nozzle flow. Further- 
more, since the major portion of the expansion takes place in 
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FULL NOZZLE LENGTH 


TEST SECTION DIAMETER 


the upstream section of the nozzle, the axial gradients in the 
nonuniform flow outside the center core will be very much less 
than at the corresponding point at the exit of a purely conical 
nozzle. The case for the short nozzle is even stronger at 
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higher Mach number, lower initial cone angle and (especially) 
larger uniform flow cores, since for these nozzles the wall angle 
at the test section would be significantly less than 3.5 deg. 
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Fig. 4 Full length of WADC contour nozzles 
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M. SEGAL! 


‘ 
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‘This paper presents a simplified method for obtaining 
the burning area of a spherically shaped grain. With this 
information the thrust-time curve can be plotted. A de- 
scriptive geometry rather than a mathematical approach 
was taken because it can be developed easily and does not 
require extensive calculations. This method will describe 
the burning characteristics not only of a spherical motor 
but also of the hemispherical or elliptical portion of a 
cylindrical motor. The application of this method is 
demonstrated for a 20-in. diameter spherical grain. 


Fens method has been developed for predicting the 
performance of spherical motors. Since this method is 
applicable to any spherical configuration, it can also be used 
to describe the burning characteristics of the hemispherical 
ends of a cylindrical motor. With minor modification this 
method can be used to predict the performance of other 
shaped grains (ellipses, cassanian ovals, etc.). 

In order to analyze the grain, two approaches are possible: 
The burning area can be calculated mathematically, or the 
burning area can be presented by means of descriptive 
geometry. 

The mathematical approach is complex, costly and time 
consuming; however, the descriptive geometry approach is 
simple, fast and inexpensive. Using this method no complex 
mathematics, expensive computer operation or high priced 
engineering talent is needed to analyze the grain. Within 
two days a draftsman can generate the information needed to 
plot the required thrust-time curve. 

The burning area at various times is obtained by projecting 
the lines which encompass a given burning area into a true 
length view. These lines can be broken up into increments, 
laid out in flat pattern and measured with a planimeter. 


Application of a Descriptive Geometry Analysis to the 
Design of a Spherical Motor 


A 20-in. OD spherical motor with a star-shaped grain, as 
shown in Fig. 1, has been selected for illustration. 


Analysis of burning process 


The first problem is to determine how the grain regresses 
during burning. Fig. 1 shows the grain pattern at a cross 
section along the equator of the sphere. It can readily be seen 
that burning proceeds in two distinct phases. 

Phase one encompasses the time from ignition to web burn- 
out; that is, to the moment when the flame first touches the 
chamber wall. During this time, since the burning areas re- 
main essentially the same (neutral burning), burning rates are 
constant and pressures and thrust can be obtained by a simple 
proportional relationship. 

In phase two (Fig. 2) the silver is burning out. Since the 
burning area is diminishing rapidly, the neutral condition 
present in phase one no longer exists. Now thrust and pres- 
sure are no longer proportional to burning area, and a. more 
involved analysis is needed. 

Attempt to determine the burning area 

Because burning of the grain occurs in two distinct phases, 
the analysis must be carried out in the two steps previously 
described. 

Phase one: Let us consider the end condition of this phase, 
the moment when the flame first touches the chamber wall. 
The total burning area at this time is analogous to the surface 


Received Oct. 29, 1959. 
1 Design Engineer, Proposal Division. 


370 


area of a pumpkin, made up of various ripples. The problem 
is to find the surface area of such a ripple. This is done b 
dividing the ripple into four equal parts (see Fig. 2), a simpli- 
fying operation which enables us to obtain the true length o! 
the three lines which encompass the desired area. Once thi 
true lengths of these lines have been obtained, they can bi 
broken up into increments and laid out into a flat pattern 
which can be measured with a planimeter. 

Let us select, as an example, a typical burnback area, such 
as burnback area “#,”’ represented by the shaded portion in 
Fig. 2. Line £,,£, isa part of a circle with the known radius 
(r:). Line £,,F;, is a cirele segment within the torus of 
radius (rz) and therefore also known. The problem then is t: 
find the true length of the remaining line E;,A;. The tru 
length of this line is shown in Fig. 3 as line £,£; in view 2 
which is a projection from view 1. The true lengths of thes 
three lines have now been obtained; if they can be spread ou! 
in a flat pattern, the area they circumscribe can be measure: 
This has been done, as can be seen in view 3. This cycle ¢ 
then be repeated for burnback areas D, C, B and A and thes 
burning areas can then be plotted into a thrust-time curve fo: 
this phase of motor firing. 

Phase two: Sliver burnout begins at the time when th: 
flame first touches the walls of the chamber. At this time, th: 
burning areas at web burnout and at the beginning of tailofi 
are equal. Pressure and, consequently, the burning rate ai 
this instance are known. However, as burning progresses, th: 
sliver diminishes rapidly. As a result the pressure within th: 


Fig. 1 Burning area of a spherical grain 


Fig.2 Quarter section of spherical grain segment showing burn- 
ing area at burnback position 
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sphere drops sharply, which in turn slows down the burning 
rate. 

Now when we trace the burning intervals during sliver burn- 
back, these distances do not represent uniform time periods. 
The burning rate is no longer constant. While this nonuni- 
form burning complicates the analysis slightly, the same 
concepts can be used to obtain the thrust-time curve during 
sliver burnout. First the burning area for each of the selected 
sliver burnback positions is determined by the method de- 
scribed in phase one. Chamber pressures and burning rates 
are then calculated by substituting these measured values in 
the mass balance equation. Therefore, the rate at which the 
sliver propellant is consumed between two burnback positions 
is the average of the burning rate at these two positions. A 
relationship has thus been established between the geometrical 
distance between burnback position and the time it takes the 
flame front to move between them. 


Practical application 


One of the requirements for the hypothetical spherical 
motor was an essentially neutral thrust-time curve. Since 
thrust and burning area are approximately proportional until 
veb burnout, our problem is to establish that at each moment 
of motor firing the burning areas are equal. 

First we assume that the burning area at web burnout 
would be equal to the surface area of the sphere. Since this 
srain must be neutral burning, we then designed the port area 
to be equal to this figure. To do this it was necessary to as- 
sume a web and ray thickness. 

Now the grain was graphically laid out and the burning 
areas were projected with the help of the Modified Gore Lay- 
out, as previously described. The layout method produces 
burning area considerably less than the assumed value but 
this must be expected since the overall purpose of the first 
grain design is to establish sliver volume, port volume, and an 
indication of burning area neutrality. 

Using this new burning area and anticipating what the new 
loading would be, a new cycle of ballistic calculations was 
made which resulted in a slightly different web thickness. 
Using this data, a final layout was made, and the desired 
thrust-time curve was plotted. 
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Fig. 3 Flat pattern layout of burning area of a spherical grain 
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Conclusion 


The approach which has been used is that of a descriptive 
geometry analysis conducted in the following way: By pro- 
jecting the lines which encompass a given burning area, the 
true lengths of these lines are determined. Breaking them up 
into increments, area segments are drawn which are then 
measured with a planimeter. 

While it is fully realized that a more sophisticated method 
based on mathematical philosophies is needed, the method 
presented will provide the preliminary design engineer with 
a fast, sufficiently accurate means for predicting the per- 
formance of spherical grains. 


Curves for Rapid Determination of 
Orbital Transfer Requirements 


PHILIP J. BONOMO! 
Westinghouse Electric Corp., Baltimore, Md. 


Nondimensional curves are presented which permit the 
determination of minimum energy requirements for orbit 
to orbit transfer in a central force field. The curves are 
general and apply to minimum energy transfers with any 
body (Earth, sun, ete.) as the primary gravitational source. 
Four examples are given to illustrate the utility of the 
curves: Transfer from a 1000-mile Earth orbit to an orbit 
in the vicinity of the moon, successive transfers from a 
200-mile Earth orbit to a 22,400-mile ‘‘stationary”’ orbit, 
and interplanetary transfer from Earth to Mars and from 
Earth to Venus (Hohmann transfer ellipses). 
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RANSFER from one Keplerian orbit to another within 
a central force field is effected by modification of the ve- 
locity vector of the transferring body. In a two-body situa- 
tion, a unique velocity, corresponding to a circular orbit, is 
associated with each radius r above the center of the primary 
gravitational source; i.e., Veirete is proportional to the square 
root of 1 over r. Variations in this velocity will modify the 
circular orbit into an elliptical orbit of major axis greater or 
less than the diameter of the circular orbit. A decrease in 
velocity will result in an interior ellipse, while an increase 
in velocity will result in an exterior ellipse. The maximum 
increase in velocity is limited, for transfer purposes, to the 
escape velocity associated with the radius of the circular orbit 
(Vescape = / 2 Vo.rcte). These facts permit transfer from one 
circular (or elliptical, if transfer is initiated at perigee or 
apogee) orbit to another and form the basis of the transfer 
requirement curves described next. 
An orbital transfer problem is specified by and the 
transfer requirements are given by k and e.? 


See Nomenclature. 
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Development of Curves | 


The velocity at any point in an elliptical orbit is given by 
[1] 
and the escape velocity associated with the radius r by 


2 
> [2] 
r r 


= 


Assume that the initial orbit is circular (or that the initial 
orbit is elliptical with transfer initiated at perigee or apogee), 
and that the velocity of the orbiting body is instantaneously 
changed; the orbiting body will then enter a transfer elliptical 
orbit. Assume further that the new total velocity Vyeq is 
some fraction k of the escape velocity associated with the 
initial radius r (k < 1.0). Then from Equations [1 and 2] 


2 
€ = k2V = (*) [3] 
r a r 


The initial radius r and the desired radius rae, will be either 
the perigee and apogee or apogee and perigee radii of the 
transfer ellipse (depending on whether rae;/r> 1 or raes/7 < 1). 
In either case, for an ellipse 

2a =r, +7. = + Paes 
Combining Equations [8 and 4] and solving for k 
Tdes 


2a 2a 


Defining \ as the desired radius to initial radius raes/r, k may 


be expressed as 
1/2 
F + 


From Equation [2] 


e=(y)~” 


where ¢€ and y are as defined below. Note that Vieq = 
k € Veso,,, Also, since raV. = for an ellipse, the velocity 
V at raes is given by 

r View 


V = Vg = 
Tdes 


[7] 


Equations [5 and 6], which are plotted in Fig. 1, together with 
Equation [7], are sufficient to determine the velocity require- 
ments for transfer from one orbit to another, as: will be il- 
lustrated. Since all quantities are nondimensional, the equa- 
tions apply for any two-body transfer situation. 


Illustrative Examples 


Four examples illustrate the application of Fig. 1. 
“4 (Continued on page 421) 
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